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DEVELOPMENT OF A COMPUTATIONAL 
MODEL FOR PREDICTING SOLAR WIND 
FLOWS -PAST NONMAGNETIC 
TERRESTRIAL PLANETS 

SUMMARY 

A summary report is provided of the work performed under 
NASA Contract No. NASW-3182. This work relates to the develop- 
ment and application of a computational model for the determi- 
nation of the detailed plasma and magnetic field properties of 
the global interaction of the solar wind with nonmagnetic 
terrestrial planetary obstacles. The theoretical method is 
based on an established single fluid, steady, dissipationless, 
magnetohydrodynamic continuum model, and is appropriate for 
the calculation of supersonic, super-Alfvenic solar wind flow 
past terrestrial ionospheres. This report provides a concise 
statement of the problems studied, a summary of all the 
important research results, and copies of the publications. 


1 . INTRODUCTION 

This is the final summary report under Contract No. NASW- 
3182 for the National Aeronautics and Space Administration. 

All of the important results from the research performed under 
this contract have been reported in the open literature, both 
in scientific journals and as technical papers at scientific 
meetings with appropriate acknowledgement of NASA support. This 
summary report provides a statement of the problem studied, a 
descriptive summary of the most important results, a reference 
list and copies of all publications resulting from the research, 
and a list of all participating scientific personnel. 
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2 . 


STATEMENT OF PROBLEM INVESTIGATED 


The problem toward which the research under this contract 
was directed was the modeling of the three-dimensional global 
interaction of the solar wind with non-magnetic terrestrial 
planet obstacles. The viewpoint on which this was to be ac- 
complished is based on a continuxim rather than particle approach. 
The theoretical method employed is based on an established 
single fluid, steady, dissipationless magnetohydrodynamic model 
that is appropriate for the calculation of supersonic, super- 
Alfvenic solar wind flow past magneto/ionopause obstacle shapes 
typical of terrestrial planets. The overall objective was the 
enablement of rational quantitative modeling studies to be per- 
formed on the global solar wind interaction problem for non- 
magnetic terrestrial planets. This accomplishment would involve 
the development and utilization of a multi-faceted computational 
procedure that embodies the theoretical model and incorporates 
advanced numerical methods so as to enable more general and 
detailed theoretical studies of the plasma and field properties 
in the interaction region than heretofore possible. Comparisons 
of predictions from the model were anticipated both with 
existing theories where possible and with observational results 
where available. Additionally, the model was to be employed in 
several collaborative efforts with other space scientists in 
which predictions from the model would both augment other theo- 
retical analysis and assist in interpreting results and com- 
parisons with observations. 

3. SUMMARY OF RESEARCH RESULTS 
3.1 Development of the Basic Computational Model 

The theoretical model employed here for predicting the 
global interaction of the solar wind with terrestrial planet 
magneto/ionospheres involves solution of the continuum partial 
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differential equations of magnetohydrodynamics of a perfect gas 
having infinite electrical conductivity and zero viscosity and 
thermal conductivity which express the conservation of mass, 
momentum, energy, and magnetic field. The computational method 
developed to solve these equations is described in detail in 
Ref. 1. In siammary, two separate finite-difference procedures 
are employed; one to treat the flow domain encompassing the bow 
shock and ionopause obstacle back to the terminator plane, and 
another to treat the flow domain downstream of the terminator. 
The various reasons for this combination are discussed in 
Ref. 1. The net result is the achievement of a computational 
efficiency not possible by any other existing means. 

One of the most important tasks of this segment of the 
study was verification of the accuracy and range of validity 
of the computational procedure. This verification was accom- 
plished through comparisons with results of other theoretical 
methods, and also by exercising the procedure on computations 
for obstacle shapes and solar wind flow conditions that span 
the entire range of interest for applications to terrestrial 
planets within the solar system. These results are reported 
in Ref. 1 and demonstrate both the accuracy and robustness of 
the computational procedures. 

Finally, one of the most significant achievements of the 
study was the incorporation of the entire computational proce- 
dure embodying the theoretical method into a modular, com- 
pletely-automated, user-friendly code. This organization allows 
any general user access to the results of the model at very 
reasonable computational cost. This feature is demonstrated in 
Ref. 1 where examples are provided of detailed maps of the 
global interaction region for all important plasma and field 
variables . 
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3.2 Extended Development and Evaluation of 
the Model with Observations 


The advanced computational model that was developed above 
was further refined and generalized. Detailed evaluations of 
the model were then made by comparisons with observational 
plasma and magnetic field data for Venus. These developments 
are detailed below. 

A new family of non-magnetic terrestrial planet obstacle 
shapes for solar wind/ionospheric interactions were determined 
which take account of gravitational variation in the scale 
height parameter. The results are detailed and reported in 
Ref. 2. An extensive catalog of solar wind flows past these 
new shapes, consisting of contour maps of all the important flow 
and magnetic field properties, were obtained and are presented 
in Ref. 3. This catalog of results has proven extremely instru- 
mental to a wide number of investigations engaged in related 
solar wind studies. 

A significant refinement included in the computational 
model in this extended development has been the allowance of 
an arbitrary oncoming direction of the interplanetary solar 
wind. This refinement is detailed in Ref. 3. This particular 
feature has been instrumental in showing that the asymmetry of 
the bow shock at the Earth previously observed and attributed 
to magnetic effects is either non-existent or at a level below 
that measurable by current spacecraft instrumentation. . These 
and extensive related results employing the model which further 
demonstrate the accuracy of both the computational procedures 
and the basic theory are reported in Refs. 4 and 5. 

Finally, a provision for the determination of time his- 
tories of plasma and field properties along an arbitrary space- 
craft trajectory has been incorporated into the computational 
model. This feature has proven invaluable in providing 
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theoretical results for direct comparisons with spacecraft 
observations. Details of this provision and results of time 
history comparisons with plasma and field data from the Pioneer- 
Venus spacecraft are provided in Ref. 2. We note in particular 
that model comparisons in Ref. 2 with magnetic field observa- 
tions for several orbits at quiet time conditions are given and 
display notably good agreement. Because determination of the 
frozen magnetic field is the last step in the computational 
model, its accurate prediction serves as a critical test of the 
correctness and accuracy of the theoretical and computational 
models . 


3.3 Collaborative Studies with the Model 

In conjunction with this investigation, a number of col- 
laborative studies have been undertaken with other space sci- 
entists in which theoretical presictions from the present model 
have been employed to interpret observational plasma and field 
results and to augment other theoretical analyses. These 
studies are described below. 

The most extensive collaborative effort has been with 
Dr. James A. Slavin and Professor R. E. Holzer at the Institute 
of Geophysics and Planetary Physics, U.C.L.A. Results from 
the current predictive model provided the theoretical basis of 
an examination of solar wind flows past all the terrestrial 
planets in the solar system. The model was applied first to 
extensive and successful modeling of the mean bow shock shapes 
and positions about these various planets. These results are 
reported in Ref. 4. Next, the model was employed to aid in 
interpreting these comparisons in light of both observations 
and previous analyses. Those results are given in Ref. 5 and 
represent the most detailed and comprehensive study completed 
to date of the global interaction phenomena of the solar wind 
with all the terrestrial planets (Mercury, Venus, Earth, and 
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Mars) from the aspect of observations of the bow shock position 
and shape. A number of significant findings and conclusions, 
several involving ongoing controversies, emerged from this study. 
These included: (1) the average bow shock position of the solar 

wind about the Earth is predicted to within 2% of its special 
location by the present model, (2) no East-West asymmetry of the 
Earth's bow shock as suggested by Romanov and others is dis- 
cernable within present observational error when appropriate 
account is taken of the average aberrated direction of the on- 
coming solar wind, and (3) the appropriate interplanetary Mach 
number to be used in such modeling studies is the sonic rather 
than magnetosonic one. All of these results both confirm and 
verify the current model's applicability to terrestrial planets. 

Further studies with Dr. Slavin and Professor Holzen have 
also been made employing the model to investigate the asymptotic 
behavior of planetary Mach cones. Mach cone angles determined 
from observational distant bow shock shapes and positions for 
solar wind flows past Venus, Earth, and Mars were compared in 
Ref. 6 with far downstream predictions from the gasdynamic 
model. The results verified that the model which is already 
known to predict good results in regions ahead of and up to the 
terminator also yields good results up to certain downstream 
distances for all the planets; i.e., -4 R^j^ at Venus, -6 R^j^ at 
Earth, and -10 R^j^ at Mars, where R^j^ denotes the particular 
obstacle nose radius. 

Discrepancies, however, appear farther downstream of these 
points presumably due to the difference between the MHD fast 
mode Mach number and the sonic Mach number inherent in the gas- 
dynamic model. The tendency to achieve better agreement further 
downstream for these various obstacles between gasdynamic theory 
and observation as the distance from the sun increases (i.e., 
Venus, Earth, Mars) is attributed to the increase in accuracy 
of the gasdynamic approximation with decreasing IMF strength. 
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These results suggest that gasdynamic theory predictions for 
the far-downstream flow about obstacles in the solar system will 
be very accurate for describing flows about the large bodies in 
the outer solar system. 

A related collaborative effort to the previous studies that 
employed the current predictive model was carried out with 
John D. Mihalov of the Space Sciences Division at Ames Research 
Center. In that effort, which formed the basis of Mr. Mihalov 's 
dissertation for Engineer's Degree from Stanford University 
(Ref. 7), the predictive model was employed in a series of com- 
parative calculations involving the detailed time histories of 
plasma and magnetic field properties in the Venusian ionosheath 
region for selected quiet time interplanetary conditions. These 
results are summarized in Ref. 8 and again verify the correct- 
ness of the present model. They also suggest some additional 
phenomena not included in the model which may be occurring 
throughout the Venusian ionosheath to retard the flow. This 
includes mass loading from photoionization and charge exchange, 
the possible presence of a viscous boundary on the ionopause, 
the leakage of ionospheric material into the solar wind flow, 
or various combinations of these effects. 

Another different series of collaborative efforts involving 
use of the predictive model was carried out primarily with 
Professor C. T. Russell and Dr. J. G. Luhmann of the Institute 
of Geophysics and Planetary Physics, U.C.L.A. as well as others 
at that institute. These studies were focused primarily on 
employing the magnetic field predictive capability of the 
present model to perform a variety of different investigations. 
In Ref. 9, a study was made of the magnetospheric source of 
energetic particles observed upstream of the Earth's bow shock. 
Calculations were performed in which those magnetosheath field 
lines predicted by the model to drape over the magnetopause 
were traced from the magnetopause to the bow shock in order to 
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locate regions at the shock that should be populated with mag- 
netospheric particles. Subsets of those field lines that 
connect to potential sites of magnetic merging on the magneto- 
pause were also traced in the event that leakage occurs 
preferentially where normal components of the field are present 
across the boundary. 

In Ref. 10, the predictive model was employed to investigate 
patterns of magnetic field merging sites on the magnetopause. 
Predictions of the magnetospheric field based on the Hedgecock 
and Thomas model and predictions of the magnetosheath field 
based on the current model were used to determine the relative 
orientations between the two fields at locations in the vicinity 
of the dayside magnetopause. Areas on the magnetopause with 
various degrees of antiparallelness between the two fields for 
various orientations of the IMF were obtained for the purpose 
of locating potential field merging sites and displayed as 
contour diagrams. The results suggest that large fractions of 
the magnetopause surface are suitable for merging for IMF's 
that are primarily southward or radial in direction. 

In Ref. 11, the predictive model was employed to investigate 
the continuing question of whether the Mars-3 spacecraft obser- 
vation of January 21, 1972 was of a Martian magnetosphere or of 
a compressed IMF in the magnetosheath. In this study, the gas- 
dynamic model was employed to generate the global flow field 
and then the magnetic field computational module was repeatedly 
employed to investigate whether an appropriate IMF could be 
determined which would produce the time history variation of 
magnetic field that was observed. Based upon the results, it 
was found that a good simulation of the observed magnetic vari- 
ation could be made purely on the basis of the magnetosheath 
field without invoking an entry into a Martian magnetosphere to 
explain the observations. 
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A final current collaborative study is underway involving 
Dr. William Knudsen of Lockheed Palo Alto Research Laboratories 
regarding Pioneer-Venus velocity potential analyzer observations 
and interpretations of the Venusian ionopause. The present 
model is being employed to examine a more accurate shape de- 
termination of the Venusian ionopause based on measured iono- 
spheric properties at axial locations from the subsolar point 
back beyond the terminator and into the wake region. Employ- 
ment of the predictive model to provide details of the plasma 
and magnetic field properties in those regions near the 
ionopause is essential for understanding the key physical 
phenomena present in the data. Presently, flow field deter- 
minations from the model have been made for a series of 
new and different ionopause models. 


9 



3.4 References 


1. Spreiter, J. R. and Stahara, S. S.: A New Predictive 

Model for Determining Solar Wind - Terrestrial Planet 
Interactions. Jour, of Geophy. Res., Vol. 85, No. A12, 
pp. 6769-6777, Dec. 1, 1980. 

2. Spreiter, J. R. and Stahara, S. S.: Solar Wind Past Venus: 

Theory and Comparisons. Jour, of Geophy. Res., Vol. 85, 

No. A13, pp. 7715-7738, Dec. 30, 1980. 

3. Stahara, S. S. , Klenke, D. , Trudinger, B. C. , and Spreiter, 

J. R. : Application of Advanced Computational Procedures 

for Modeling Solar-Wind Interactions with Venus - Theory 
and Computer Code. NASA CR 3267, May 1980. 

4. Slavin, J. A. and Holzer, R. E. : Solar Wind Flow About 

the Terrestrial Planets 1. Modeling Bow Shock Position 
and Shape. Jour, of Geophy. Res., Vol. 86, No. A13, 
pp. 11401-11418, Dec. 1, 1981. 

5. Slavin, J. A., Holzer, R. E. , Spreiter, J. R. , Stahara, 

S. S., and Chaussee, D. S. : Solar Wind Flow About the 

Terrestrial Planets 2. Comparison with Gas Dynamic Theory 
and Implications for Solar-Planetary Interactions. Jour, 
of Geophy. Res., Vol. 88, No. Al, pp. 19-35, Jan. 1, 1983. 

6. Slavin, J. A., Holzer, R. E., Spreiter, J. R. , and Stahara, 

S. S. : Planetary Mach Cones: Theory and Observation. 

Submitted to Jour, of Geophy. Res. 

7. Mihalov, J. D.: Comparison of Gas Dynamic Model for Solar 

Wind Flow Around Venus with Pioneer Venus Orbiter Data. 
Stanford University Engineers Thesis, May 1981. 

8. Mihalov, J. D. , Spreiter, J. R. , and Stahara, S. S.: 
Comparison of Gas Dynamic Model with Steady Solar Wind 
Flow Around Venus. Jour, of Geophy. Res., Vol. 87, No. A12 
pp. 10363^10371, Dec. 1, 1982. 

9. Luhmann, J. G. , Walker, R. J. , Russell, C. T. , Crooker, 

N. U., Spreiter, J. R. , Stahara, S. S. , and Williams, D. J. 
Characteristics of the Magnetospheric Source of Inter- 
planetary Energetic Particles. Institute of Geophysics 
and Planetary Physics, U.C.L.A. Publication No. 2442, 

May, 1983. Submitted to Jour, of Geophy. Res. 


10 



10 . 


Luhmann, J. G. , Walker, R. J. , Russell, C. T. , Crooker, 

N. U., Spreiter, J. R. , and Stahara, S. S.: Patterns of 

Magnetic Field. Merging Sites on the Magnetopause. 
Institute of Geophysics and Planetary Physics, U.C.L.A. 
Publication No. 2443, May 1983. Submitted to Jour, of 
Geophy. Res. 

11. Russell, C. T., Luhmann, J. G. , Spreiter, J. R. , and 

Stahara, S. S. : The Magnetic Field of Mars: Implications 

from Gas Dynamic Modeling. Institute of Geophysics and 
Planetary Physics, U.C.L.A. Publication No. 2468, July, 
1983. Submitted to Jour, of Geophy. Res. 


4 . PUBLICATIONS LIST 


Stahara, S. S., Klenke, D. , Trudinger, B. C. , and Spreiter, 

J. R. : Application of Advanced Computational Procedures for 

Modeling Solar-Wind Interactions with Venus - Theory and 
Computer Code. NASA CR 3267, May 1980. 

Spreiter, J. R. and Stahara, S. S. : A New Predictive Model for 

Determining Solar Wind - Terrestrial Planet Interactions. Jour, 
of Geophy. Res., Vol. 85, No. A12, pp. 6769-6777, Dec. 1, 1980. 

Spreiter, J. R. and Stahara, S. S. : Solar Wind Past Venus: 

Theory and Comparisons. Jour, of Geophy. Res., Vol. 85, No'. A13, 
pp. 7715-7738, Dec. 30, 1980. 

Mihalov, J. D. , Spreiter, J. R. , and Stahara, S. S. : Comparison 

of Gas Dynamic Model with Steady Solar Wind Flow Around Venus. 
Jour, of Geophy. Res., Vol. 87, No. A12, pp. 10363-10371, 

Dec. 1, 1982. 

Slavin, J. A., Holzer, R. E., Spreiter, J. R. , Stahara, S. S., 
and Chaussee, D. S.: Solar Wind Flow About the Terrestrial 

Planets 2. Comparison with Gas Dynamic Theory and Implications 
for Solar-Planetary Interactions. Jour, of Geophy. Res., 

Vol. 88, No. Al , pp. 19-35, Jan. 1, 1983. 

Slavin, J. A., Holzer, R. E., Spreiter, J. R. , and Stahara, 

S. S. : Planetary Mach Cones: Theory and Observation. Sub- 

mitted to Jour, of Geophy. Res. 


11 



Luhinann, J. G. , Walker, R. J., Russell, C. T. , Crooker, N. U. , 
Spreiter, J. R. , Stahara, S. S., and Williams, D. J. : 
Characteristics of the Magnetospheric Source of Interplanetary 
Energetic Particles. Institute of Geophysics and Planetary 
Physics, U.C.L.A. Publication No. 2442, May, 1983. Submitted 
to Jour, of Geophy. Res. 

Luhmann, J. G. , Walker, R. J. , Russell, C. T. , Crooker, N. U. , 
Spreiter, J. R. , and Stahara, S. S.: Patterns of Magnetic 

Field. Merging Sites on the Magnetopause. Institute of Geo- 
physics and Planetary Physics, U.C.L.A. Publication No. 2443, 
May, 1983. Submitted to Jour, of Geophy. Res. 

Russell, C. T. , Luhmann, J. G. , Spreiter, J. R. , and Stahara, 
S. S. : The Magnetic Field of Mars: Implications from Gas 

Dynamic Modeling. Institute of Geophysics and Planetary 
Physics, U.C.L.A. Publication No. 2468, July, 1983. Submitted 
to Jour, of Geophy. Res. 


5. PARTICIPATING SCIENTIFIC PERSONNEL 


(1) Stephen S. Stahara - Principal Investigator 

(2) John R. Spreiter - Co-Principal Investigator 


12 



APPENDIX 


COPIES OF PUBLICATIONS 



JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 85, NO. A12, PACES 8769-6777, DECBtBER 1, 1980 
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Abstract ■ A coaputatlonal nodel has been 
developed for the determination of the gasdynaaic 
anc s^gnetic field properties of the solar vind 
flov around a magnetic planet, such as the earth, 
or a nonmagnetic planet, such as Venus. The 
proceoures are based on an established single- 
fluic, steady, dissipationless , magnetogasdynamic 
model and are appropriate for the calculation of 
axiEj’mmetnc , supersonic, super-Alfv6nic solar 
vind flow past a planetary magneto/ionosphere. 
Sample results are reported for a variety of 
solar wind and planetary conditions. Some of 
these are new applications ; others are included 
to show that the new procedures proouce the same 
results as previous procedures when applied to 
the same conditions. The new methods are 
completely automated and much more efficient and 
versatile tnan those employed heretofore. 

Introduction 

Tne magnetogasdynamic model for the inter- 
action of the solar vine and a planetary magneto- 
sphere, or ionosphere if the planet is nonmag- 
netic, has been found to be of great utility in 
the prediction and Interpretation of observations 
in space. The development of this model has a 
long history stesmlng from the work of Sprelter 
and Jones [1963], Dryer and Faye-Petersen [1964, 
1966], Sprelter and Sunaers [1966], AlRsne [1967], 
Dryer and Keclctan [1967], Sprelter et al. [1968, 
1570], Alksne and Webster [1970], Sprelter and 
Alksne [1969, 1970], and Sprelter and Rizzi 
[1974]. Although the usefulness of the original 
results is veil established, there is much to be 
desired. Results have been calculated for only 
a limited set of solar vind and planetary condi- 
tions. Moreover, the original solutions bordered 
— on wnat vas-barelyvposslble at the time ,_ required 
considerable hand computation, and were extremely 
laborious. Improvements In both numerical 
methods and computer capabilities have now 
rendered such procedures thoroughly obsolete, 
both from the standpoint of efficiency of calco- 
latlor, and generality of application. 

To remedy this aituation, a new, efficient 
procedure has been developed that provides the 
numerical solution for the streamlines, magnetic 
field lines, and contours of density, speed, 
temperature, and magnetic field intensity 
similar tc those familiar from the original 
literature. Tne basic theoretical model employed 
Is that des Tibed by Spreiter, Alkane, and their 
colleagues in the work cited above. The present 
solution makes full use of recent advances in 
nonlinear computational fluid dynamics. These 
advances Include an algorlthB of beam and Warming 
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[1976] to determine the gasdynamic flov near the 
nose of the magneto/lonopause. The raiainder of 
the flov field is calculated by using a method 
developed by Kutler et al. [1973] and by Chaussee 
et al. [1975]. The magnetic field is also calcu- 
lated quite differently than it was previously, 
being baaed on the decomposition of Alkane and 
Webster [1970]. This paper presents an outline 
of the principal features of the analysis and a 
sample of the results. Complete details of the 
computational procedures are provided by Stahara 
et al. [1977]. 


Mathematical Model 


The fundamental assumption of the present 
work and all of the geophysical work cited above 
is that Che average bulk properties of solar 
wind flow around a planetary magneto/ionosphere 
can be described adequately by solutions of the 
continuum equations of magnetogasdynamics of a 
perfect gas having infinite electrical conduc- 
tivity and zero viscosity and thermal conduc- 
tivity. These equations, which express the 
conservation of mass, momentum, energy, and 
magnetic field, are given by the following 
expressions: 


3t 





0 


( 1 ) 


E B 

± ^ (pv^v^ . pd^^ - 

&T. *ik 4xG “ 6xG *ik^ “ ® 


(2) 


^ + pe + p6 + 1^)- 


SB 


pVj^(-2|- + e + £ + t) + s 




sr->T'Vk- Vj>' 

K 




( 3 ) 


(4) 


where 



\ ■ r- - \\\> 
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given by 


P 



( 6 ) 


Paper number 60A1121 . 
o:4f-o:-:'*.'60/060Aii2isoi.oo 



Sprelter and Stahara: Solar Wlnd/Flanetary Interacclonc Model 


In accordance vlch standard usage the sytabols d, 
p, V, T, e • C^T, and h « e+p/c “ C T refer to 
the density, pressure, velocity, temperature, 
internal energy, and enthalpy; and C refer to 
the specific heats at constant voluae and pres- 
sure; R “ (C “C^)w “ 8.31x10^ ergs/g *K is the 
universal gal constant, and t represents the 
eean nolecular weight nondiaensionaliaed so that 
t • 16 for atonic oxygen. For fully ioniaed 
hydrogen, n is thus 1/2. The magnetic field B 
and the Foynting vector S for the flux of electro- 
magnetic energy are in Gaussian units. The 
gravitational potential 6 and acceleration | are 
assumed to be due to massive fixed bodies so that 
their time derivatives are zero. Because of the 
omission of dissipative terms in these equations, 
surfaces of discontinuity may develop in the 
solution, across which the fluid and magnetic 
properties change abruptly, but in such a way 
chat mass, momentum, magnetic flux, and energy 
are conserved. These are approximations to 
comparatively thin layers across which similar 
changes in the fluid and magnetic properties 
occur in- the corresponding theory of a dissi- 
pative gas and correspond physically to the bow 
wave, magneto/ionosphere boundary, and possible 
otner chin regions of rapidly changing properties. 

The shape of the magneto/ionospheric obstacle, 
assumed to be axisysisetric about a line through 
the planet center extending parallel to the 
free- stream solar wind velocity vector, may be 
specified in either of two ways: by a numerical 
table of coordinates or by use of simplified 
models from the cited references, which enable 
tne calculation of approximate magneto/ ionopause 
shapes upon specification of a small number of 
solar wind and planetary magnetic field or iono- 
spheric parameters. 

For the earth or other magnetic planets having 
a cipole field this shape is obtained by rotating 
the equatorial trace of the three-dimensional 
magnetosphere surface determined by balancing 
the Newtonian pressure p > Ro cos^Ui against 
the magnetic pressure due to If times ^be tangen- 
tial component of the dipole field. Here, K and 
f are constants usually taken to be unity (see 
Spreiter [1976] for a discussion), c and V are 
the density and velocity of the solar wind, the 
subscript » refers to conditions in the incident 

-streM-upstream of the bow wave, and- is the 

angle between V and the normal to the magneto- 
pause. These considerations lead to the 
following differential equation for the geocen- 
tric distance to the magnetopause r as a func- 
tion of geomagnetic longitude t, measured in 
such a way that ^ is equal to r,/2 at the magneto- 
sphere nose and Increases to 3 t/ 2 as one moves 
arounc the flank of the magnetosphere to the 
remote tail [Beard, 1960): 


df 
c 

d; 


r 

c 


-6 

r sin C COE 
c 


-6 2 ^ 
r cos d 
m 



(7) 


*/2 < d * 3e/2 


In this equation, r^^ " ’ 

D • c (f‘E‘ /2eK£. l6 the geocentric 

e eq “ “ 


distance to the magnetosphere nose, where a^ 
refers to the radius of the earth, and is the 
average Intensity of the field at >the geomagnetic 
equator . 

For a nonmagnetic planet having a sufficiently 
dense ionosphere to withhold the solar wind , such 
as Venus appears to be most of the time, the 
ionopause shape is determined by balancing the 
Newtonian pressure and the ionospheric pressure 
approximated by 


P 


Pr exp 


r-r„ 




( 6 ) 


oerivea trom (i)-(5) by assuuing steady condi- 
tions, negligible velocity and magnetic field, 
and a constant ionospheric scale height H»kT/ag, 
where k = 1.36 x 10"^^ ergs/°K is Bcltzmann's 
constant, m is the mean molecular mass, and Pj^ 

IS the pressure at some reference radius ij^. 

These considerations lead to the following differ- 
ential equation for the coordinates tj^ of the 
ionopause as a function of angle £ from the 
ionopause nose as measured at the planetary 
center [Spreiter et al., 1970]: 



r . 
1 


sin 2 c 
2 (.' ■ 


r 


sin c) 


0 


in which h • exp [-(r^ - Ro)/H), where R^, is the 
value of rj at the ionopause nose where £ « 0. 

A family of ionopause shapes can be obtained by 
integrating (9) for different values of H/Rp. 
Appropriate values of H/R^ for application to 
Venus and Mars appear to be in the range from 
0.01 to 0.30. 

Two important parameters which characterize 
the solar wind flow in the magnetogasdynamic 
model are the free-stream values for the Mach 
number M « V/a and the Alfv6n Mach number 

“ V/A, in which a “ (y p/o)^^^ is the speed 
of sound, A is the AJ.fv6n speed A ■ 
and y is the ratio of specific heats, usually 
taken to be 5/3 in solar wind applications. 
Because is typically greater than 5 in the 

so lar w ind, an important simplification may b€_ 
introduced on the basis that for large the 
magnetic terms in the momentum and energy equa- 
tions (2) and (3) are generally sufficiently 
small compared with the gasdynamic terms that 
they may be disregarded therein. This decouples 
the solution for the fluid notion from the 
magnetic field, which may be calculated subse- 
quently by using the induction equation (1) and 
values for V already calculated. 

Although few detailed comparisons have beer, 
made, other chan those of Spreiter and Rizzi 
[1976], which mode use of their complete nagneto- 
gasdynamic numerical solution for the special 
case of aligned flow, it is generally regarded 
that the magnetic field calculated in this way 
repruiuces the broad characteristics of observa- 
tions. However, the process indicates magnetic 
field strengths chat may be coo large laDedlstely 
outside the magnetopause. Zwan and Wolf [1976] 
have made an interesting analysis of conditions 
in this region and concluded chat the strong 
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nagnetlc fields vould drive the plaeaa out along 
the field lines and produce a Chin region of 
depleted plaeaa density. Crooher and Siscoe 
[1977] have shown further that an observed 
anisotropy of Che magnetosheath pressure occurs 
naturally from a mirror instability of the 
resulting flov along the magnetic flux cubes. 
Neither of these properties nor any representa- 
ticr. of the long discussed possibility of 
macr.etic merging at the magnetopause are captured 
bv tne present model, but the ready ability to 
genera our solutions for any specified condi- 
tions should aid in further study of such pheno- 
mena ■ 


Calculation of the Casdynamic Tlov Properties 

The gasdynasic equations for supersonic flov 
past a blunt-nosed obstacle representative of a 
planetary nagneto/ionospher e are solved most 
effectively by using one metnod for the nose 
region, where the flov accelerates from subsonic 
to supersonic speeds, and another for the 
remaining region, where the flov is purely 
supersonic. In the present work the solution 
for the nose region is calculated by using a 
new axisymmetric implicit unsteady Euler equation 
solver (TK?), which determines the steady state 
sciution by a time-marching procedure. The 
remainder of the solution is calculated by 
using a shock-capturing marching procedure (SCT), " 
which spatially advances the solution downstream 
as far as required by solving the steady Euler 
equations . 

Implicit Unsteady Euler Eouation Solution 
for the Nose Region 

The partial differential equations employed in 
tne implicit code are the unsteady gasdynamic 
Euler equations derived from (l)-(3) by deleting 
the magnetic and gravitational terms and special- 
iiing for axisysietric flows. Upon introduction 
of the generalized independent variable trans- 
formetlcn t - T, C - i(T,X,R), these equations 
ma> be written as 

avj). 4 [(c^u + CjjE 4 

4 [(r^U 4 rijjE 4 G - 0 ^ (10) 

where 
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end J - - Cj-hj. is the transformation 

Jacobian. In these equations, T denotes time, 
X the axial downstream coordinate, and R the 


cylindrical radial distance; u and v are the 
velocity components in the X and R directions; 

«t “ p/Ip(T“ 1)) 4 (u^ 4 v^)/2 and ^subscripts 
denote partial derivatives with re'spect to the 
Indicated variable. 

The analysis commences by introducing a compu- 
tational mesh in polar (r,B) coordinates such 
that one family of coordinates consists of rays 
from the planetary center spaced at equal incre- 
ments of B measured from the obstacle nose and 
the other of curved lines intersecting each ray 
so as to divide the portion of it between the 
magneto/ionopause and the shock wave into a 
fixed number of equal segments. The coordi- 
nate transformation equation is then used to 
map the portion of the X,R,T physical space 
bounded by (1) the bow wave, (2) the downstream 
outflow boundary at B “ r/2, (3) the obstacle 
surface, and (4) the stagnation streamline at 
6 “ 0 into a rectangle in the (, n, t computa- 
tional space as Illustrated in Figure 1. Gener- 
ally, the transformation metrics at each time 
step are not known beforehand and must be deter- 
mined numericelly as part of the solution. 
Integration step size is established by using 
the eigenvalues of the^Jacobian matrices A and E, 
where A “ 3E/3U, B » iP/aU, and also where 
V - U/J, £ - (CxU 4 txE + Cpf/J, and 
F « (nxU 4 tij{E 4 npF)/J. 

Boundary conditions for a properly poaed 
mathematical problem are chat the flow satisfy 
the axlsymaetric Rankine-Hugoniot shock relations 
derivable from equation (10) along item 1 above, 
be entirely supersonic along item 2 above, be 
parallel to Items 3 and 4 above, and be symaietric 
about item 4 above. Initial conditions are 
determined by use of an approximating formula 
for the coordinates of the bow shock wave, 
dependent on y, and the shape of the obstacle, 
and by prescribing a Kewtonian pressure distri- 
bution on the obstacle. The remainder of the 
initial properties on the obstacle surface can 
then be determined from the conditions that the 
flov must be tangential to the surface and 
possess constant entropy. A linear variation 
for the flov properties between the bow shock 
and the obstacle is then prescribed. This 
provides Che Initial flov field, which is then 
integrated in a time asymptotic fashion until 
the steady state solution is obtained. 

The basic numerical algorittas used in the IMP 
code was developed by Beam and Warming [1976] 
and is second-order accurate, noniterative and 
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Fig. 1. Transformation from physical dooiain to 
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Fig. 2 . Coaparison of 
by present Iffipllclt nethod with other techniques 
and experinent for supersonic flow past a sphere; 
M » 4 . 926 , t » 1 . 4 . (a) Stagr»atioti line density, 
(b) Surface pressure, (c) Shock standoff 
distance. 


S is a fourth-order snoothing tens introduced 
to elli&inate nonlinear stabilities that aay arise, 
since the use of central differences in the 
spatial directions results in a neutrally stable 
algorithm. 

At the boundaries, Dodlfication of the differ- 
encing algoricho to account for the particular 
conditions described above is accomplished as 
follows. The obstacle surface flow tangency 
condition is incorporated through the use of 
Kentaer's [1970] scheme, while at the symmetry 
plane the variables are reflected according to 
whether they are odd or even. At the outflow 
boundary, where the flow is entirely supersonic, 
the dependent variables are determined by extra- 
polation from the adjacent interior points. For 
the upstream boundary, formed by the bow shock 
wave, the sharp discontinuity approach of Thomas 
et al. [1972] is used. The interior flow field 
bounded by these various boundaries is treated 
in shock-capturing fashion and therefore allows 
for the correct formation of secondary internal 
shocks should any occur. 

Shock-Capturing Karchlnt Solution 
for the Downstream Solution 


The shock-capturing technique of Kutler et al. 
[1973] and Chaussee et al. [1975] employed herein 
is based on the set of four equations for steady 
axlsysaetrlc flow obtained from <10) and (11) by 


spatially factored. In particular, the 'delta 
form' with Euler time differencing is employed. 
When applied to (10), the algorithm assumes the 
form 

(I + At4,.A®)(I + AT 6 B®)(U®^^ - C*^) 

4 h 

- - AT(d,.f” + d F° + C) (12) 

4 n 

where A and B are the Jacobian matrices, 1 is the 
identity matrix, 6 and 6 are^second-order, 
central difference'operatSrs, - O(nbT), and 

iT is the integration step sire. 

Equation (12) is solved at the interior points 
only. It requires two 4x4 block tridiagonal 
inversions at each time step o^ the J.ntegr8tlon. 
The solution proceeds as follows: 

1. Define AU « - u“. 

2. Form the right-hand side of (12), and store 

results in the array - 

3. Apply smoothing D “0 - (c/8)S/J. 

4. Define 0 • (1 + At6 b”)AU, and solve the 

matrix equation (1 + Aid A )0 » {'”^1 for u, 

^ ,%n+l 

storing the result in the u array. 

5- Solve the matrix equation 
(I + ATd^B“)AU - for AU. 

6. Obtain the values of from the 

-elation 0 • AU + U . 

7. Transfer the contents of to 0®, and 

repeat all steps until satisfactory convergence 
is attained. 
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Fig. 3. Comparison of implicit and inverse 
methods for shock shape and sonic line location 
and density distribution along bow shock and 
magnetosphere boundary for • 6 , y » 5/3 flow 
past the rotated equatorial trace of the magneto- 
pause. (a) Shock shape and sonic line location, 
(b) Density distribution. 
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•ettlng the t derivatives to zero. The fourth 
of this set of equatioss, which rcpresentE 
co&aervatloo of total eaergj pe^, can be 
integrated to obtain the relation h “e^ + p/p - 
const for the total enthalpy per unit mass. 

The coBputational nesh is defined by lines of 
constant X and (r-r^) / (R -E^) , where R^ and R^ 
are functions of X that describe the radial 
cylindrical coordinates of the oagneto/ionopause 
a^ now shock wave at the saae X as the field 
point (X,R) illustrated in Figure 1. The three 
renaining partial differential equations for 
conservation of mass and of axial and radial 
masentus are then transformed to a rectangular 
computational space by C “ X, n “ (R-R^)/(R-R^) 
to obtain 

|I + |I + 5 - 0 (13) 

oC Bn 

where 

i ~ I 

f - (F - (R, - 

G - G + -r — 


(U) 

The finite difference counterpart of (13) is 
integrated with respect to the hyperbolic coordi- 
nate f to yield values for the conservative 
variable X. Subsequent to each integration step 
the physical flow variables p, p, u, and v must 
be decoded from the components e^ of X. This 
necessitates the solution of four simultaneous, 
nonlinear equations and is facilitated by using 
the relations v - p " e^ - e^u, p - e^/u. 


22 - 



Fig. 4. bow shock and embedded shock locations 
for solar wind flow with • 5, y ■ 5/3 past the 
rotated principal meridian of the magnetopause. 



Fig. 5. Magnetopause pressure coefficients for 
the principal meridian magnetopause shapes shown 
in Figure 4. 


and h “ y/ (y-i) (p/p) to determine the following 
quadratic equation for u: 
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rfl 

2(y-l) 
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+ 



h + 

t 



0 (15) 


Two roots exist; one corresponds to subsonic flow 
and is discarded, since u is always supersonic 
in the presenr application, while the other corre- 
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Fig. 6. Variation of shock standoff distance 

with oncoming Mach number and ratio of specific 

beats for various magneto/ioncpause traces as 

determined by the present implicit procedure . 

(e) Magnetopause equatorial trace, (b) lonopause 

trace -E/R^ “ 0.01. (c) lonopause crece 

-B/R - 0.5. 
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flovB past the rotated eqxiatorlal trace of the 
magnetopauae; conblned near (bluet body) aad 
far (marching) aolutlons. 

spo&de to super Bonic flov aad givea Che desired 
solution. 

Since only the bov shock wave is treated as a 
sharp discontinuity and any others that may be 
present are 'captured' by the difference algo- 
rlttc .selection of the appropriate finite 
differ ujce scheme to advance the calculation in 
the ( direction is of prime importance. As in 
the blunt body aerodynamic analysis of Kutler et 
al. [1973] and Chaussee et al. [1973], the numer* 
ical integration of (13) is accomplised by using 
the finite difference predictor-corrector scheme 
of MacCcrmack [1969], the most efficient second- 
order algorithm for shock-capturing calculations. 

Calculatioro of the Streamlines 

The streamlines are determined by integrating 
trajectories through Che known velocity field, as 
this procedure was found to be more accurate than 
the alternative mass flov calculation. The calcu- 
lation of 8 particular streamline is Initiated 
at the bov shock, where its slope is calculated 
with an exact gasdynamic relation contained impli- 
citly in both Che blunt body (IMP) and marching 
(SCT) solution, and continued by stepwise inte- 
gration in Z using a modified third-order Euler 
predictor-corrector method. Bivariate linear 
interpolation free the flov field grid points is 
employed to obtain Che velocity components (u,v) 
required at the stepwise points along the stream- 
line trajectory. 

Calculation of the Magnetic Field 

With the flov properties known from the gasdy- 
namic solution, calculation of the magnetic field 
E proceeds by integrating tbe steady state 
counterpart of (4), commonly interpreted as indi- 
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Fig. 8. Shock shape for - 8, y " 5/3 flow 
past an ionopause shape with B/k • 0.1: combined 
near (blunt body) and far (marching) field solu- 
tions . 
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Pig. 9a. Streamline nap for » 8.0, y ■ 5/3 
flow past the rotated equatorial trace of the 
magnetopause: combined blunt body and marching 
flov field. 

eating that the field lines move with Che fluid. 
This leads to a straightforward calculation in 
which Che vector distance from each point on an 
arbitrarily selected field line to its corre- 
sponding point on an adjacent field line in the 
downstream direction is determined by numerically 
integrating /V dt over a fixed time interval 6t. 
Once the coordinates of the field lines are 
determined, the magnetic field at any point may 
be calculated from the relation 

where tt is the vector 
length of a small element of a flux cube. 

Such a procedure is valid generally, but its 
use in the present calculations is confined to 
only one component of the magnetic field. The 
other components are determined by use of the 
Alkane and Webster [1970] decomposition dependent 
upon the azlsymmetrlc properties of the gasdynamic 
solution and tbe linearity of the magnetic equa- 
tions. 

Calculation of the Contour lines 

Contours are calculated for nondlmenslonalieed 
velocity |v density p/p^, and magnetic field 



Fig. 9b. Velocity map for " 8.0, y " 5/3 flow 
past Che rotated equatorial trace of tbe magneto- 
pause: combined blunt bodv and marching flov 
field. 
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roEoonentt and by application of 

1 jfjrching procedure. After a value for the 
rec: 3 -r line is specified, the shock boundary is 
*firched for intervals which bracket the selected 
;Jter one such point is located by intex- 
,.:ior.. the reaainder of the contour is deter- 
. •: o\ 'walking' around the contour, searching 

i»;r. step for the interval and then interpo- 
.r.t tc find the point throueh which the contour 
, passes- This is rep>eated until a boundary 
IS reached. The closed contours are 
ir. a similar manner, linear interpolation 
.1 „sed throughout the process. The coordinates 
' tie contour lines are output as listings, pen 
-;:;s. or both. 

Results and Discussion 

Tc verify the correctness of the new procedures 
a-.; to eeaonstrate their capability for calcula- 
solar wind flows for a variety of conditions, 
a large number of test cases have been run, and 
i sample have been evaluated by comparison with 
t:e-. lously available theoretical and experimental 
ret.lts Figure 2 presents such a comparison for 
s.persor.ic flow of air (y“ 1.^) past a sphere at 
■ -.926. The variation of density along the 
siagnaiion streamline is provided in Figure 2a, 
sT..le tne variations of surface pressure p^ and 
sho:> standoff distance i with angular distance 
!:x the nose are given in Figures 2b and 2c. In 
iirsc and similar comparisons the results provided 
5. :ne new procedures are in essentially perfect 
agreement with those of previous numerical solu- 
tions and ctxperiment. 

Figure 3 exhibits a comparison of results 
predicted by the present method with those 
sriginally calculated by Spreitcr el al. (1966, 
T96E] for the same solar vind and magnetosphere 
conditions, using the inverse method of Inouye 
anc Lomax [1962] for the blunt nose region and 
:he method of characteristics for the downstream 
supersonic region. Figure 3a displays the loca- 
tions of the bow wave and sonic line for flow 
past the equatorial trace of the magnetopause 
for ■ 6 and y » 5/3. The density distribution 
along the magnetosphere boundary and along the 
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Fig 9c. Temperature map for " 8.0, y • 5/3 

flow past the rotated equatorial trace of the 
magnetopause ; combined blunt body and marching 
flow field. 



Fig. 9c. Density map for M •= 8.0, y “ 5/3 flow 
past the rotated equatorial trace of the magneto- 
pause: combined blunt body and marching flow 
field . 

shock wave are given in Figure 3b. Essentially 
perfect agreement is obtained between the present 
implicit technique and the inverse method . 
Similarly perfect agreement has been obtained 
with all of the results presented previously in 
the cited references by Spreiter and his 
colleagues. The present method is, however, both 
more efficient and more versatile. 

Figures 4 and 5 present new results that illus- 
trate the geometric flexibility of the present 
solution and its ability to capture embedded 
shock waves as well as the bow shock, a feature 
which the inverse method cannot duplicate. This 
case is for ■ 5 and y ■ 5/3 flow around the 
axisyiBDetric shape generated by rotating the 
principal meridian of an indented magnetopause 
about its axis and may be compared with an experi- 
mental test by Spreiter et al. [1968]. This 
particular profile shape, derived by Spreiter 
and Briggs [1962] frcm the Beard [1960] approxi- 
mation to the classical Chapman and Ferraro [1931] 
theory, contains a pronounced dent with a concave 
corner in the vicinity of a magnetic neutral 
point. Spreiter and Summers [1967] argued that 
the presence of an embedded shock wave would make 
it impossible for the magnetosheath plasma to 
follow such an indented contour and that there 
would form instead a free surface approximating, 
im the meridian plane, a tangent line across the 
indentation and capping an embedded cusp region. 
Figure 4 presents results for both the indented 
and the free surfaces, denoted by solid and dashed 
lines. Fot the indented surface a shock wave 
is located on the body at approximately £ • 80', 
just as in the experiment, while for the free 
surface the pressure coefficient displays an 
approximately constant value as anticipated. 
Finally, we note that the calculation of a 
supersonic flow with an ODbedded shock and subse- 
quent subsonic pocket provides a severe test f ct 
any blunt body procedure. The ability of the 
present code to provide convergent results for 
such a flow demonstrates the capability for 
further extension and application to more general- 
ised profiles than was heretofore possible. 

As a final illustration of the range of condi- 
tioni for which the implicit code has been tested. 
Figure 6 sumarltes results for cbe shock standoff 
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distance provided by the gasdynaslc calculations 
for three values for y for three different 
magneto/ lonopause shapes and a range of values 
for from 3 to 25. These conditions span the 
range of Interest of geometrical and solar vlnd 
conditions for which the computer programs devel- 
oped herein would be normally applied. 

In order to confirm the ability of the march- 
ing code to continue the solution to some arbi- 
trarily specified downstream location, solutions 
have been carried far downstream for a wide 
variety of cases typical of solar wind Inter- 
actions. Figure 7 exhibits the location of the 
bov wave downstream from the nose region to 
x/D ■ -11 for flows with y - 5/3 and ■ 3, 8, 
and 25 past the rotated equatorial trace of the 
magnetopause. For this calculation, starting 
renditions for the marching code (SCT) are provi- 
ded by the blunt body (l>fP) code on the line 
x/D “ 0.0, which is the usual location at which 
the two solutions are joined by the computer 
program. The marching code then determines the 
remainder of the flow field back to the speci- 
fied downstream location. Since for the shapes 
considered herein the downstream flows are quite 
smooth, the marching calculation is very effi- 
cient (less than 30s, Q3C 7600, OPT«2 compiler). 
Similar results are presented in Figure 8 for 
M “ 8 and y » 5/3 flow past an ionopause shape with 






trace o£ the magnetopause. 


b7r^" 0.1. Those results have been carried to 
x/Rp»-20 in keeping with the observation that R.^ 
for Venus and Kars is only slightly greater than 
the planetary radius, whereas for the earth, D is 
of the order of 10 earth radii. 

To illustrate the capability of the present 
procedures to determine streamlines, contour maps 
of various flow properties, and magnetic field 
lines and contours, as well as to demonstrate the 
automated plotting capability for displaying these 
results. Figures 9 and 10 have been prepared. 
Figure 9 illustrates the computer-generated 
streamline locations and contour maps of velocity 
ratio V/V , temperature ratio T/T , and density 
ratio c/o for the complete near- .and far-field 
flow about the equatorial trace of the magneto- 
sphere for " 8 and y«5/3. Based on this gasdy- 
namlc solution. Figure 10 exhibits the corre- 
sponding results for the magnetic field (B/B^)^ 
and (3/B^)„ In cne plane of magnetic symmetry for 
B^ parallel and perpendicular to V^. In addition 
to demonstrating the overall smoothness of the 
computed results these two figures illustrate the 
ability of the present techniques to provide the 
completely automated production of report quality- 
plots of both gasdynamic and magnetic field 
properties for solar wind flows past axisysmetric 
magneto-ionopause shapes. 

Concluding Remarks 

Advanced numerical techniques have been 
applied to produce an efficient operational 
computer solution for the well-established dissi- 
pationless magnetogasdynasic model for axisyssse- 
tric, supersonic, super-Alfvdnic solar wind flow 
past either magnetic or nonmagnetic terrestrial 
planets. The solution consists of the following 
assemblage of computer codes; (1) blunt body code, 
to determine the gasdynamic solution near the 
magneto/ ionopause nose, (2) marching code, to 
determine the gasdynamic solution downstream of 
the magneto/lonopause nose, (3) streamline code, 
to determine coordinates of flow field stream- 
lines, (^) magnetic field code, to determine 
frozen-in magnetic field, (5) contour code, to 
determine coordinates of contour lines of flow 
and magnetic field properties, and (6) plotting 
code, to plot selected flow and magnetic field 
results. 

Comparisons are reported which demonstrate 
the accuracy of the present techniques by compar- 
ison with previously established theoretical 
methods and with experimental data. New results 
are presented for a variety of solar wind flows 
which illustrate the flexibility and generality of 
the methods. 
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A-strett . Advanced computational procedures 
arf applied to an improved model of solar wind 
y.cK past Venus to calculate the locations of the 
ionopause and bow wave and the properties of the 
floving ionosheath plasma in the Intervening 
re£:cr.. The theoretical method is based on a 
s:ngle-fluid, steady, dissipationless, magneto- 
'-.virodynemic continuum model and is appropriate 
f:r tne calculation of axisymaetrlc supersonic, 
scper-Alfvenic solar wind flow past a nonmagnetic 
tlanet possessing a sufficiently dense ionosphere 
1 C stand off the flowing plasma above the subsolar 
pcir.t and elsewhere. Determination of time 
-isicries of plasma and magnetic field properties 
iiong an arbitrary spacecraft trajectory and 
provision for an arbitrary oncoming direction of 
int interplanetary solar wind have been Incor- 
tcraied in the model. An outline is provided of 
tne underlying theory and computational procedures, 
end sample comparisons of the results are pre- 
sented vith observations frott the Pioneer Venus 
critter . 


Introduction 

Tne magnetohydrodynamlc model of Spreiter et al. 
for solar wind flow around the Venus Iono- 
sphere, stemming from the closely related earlier 
studies for the earth and its magnetosphere 
[Spreiter, 1976; Spreiter and Jones, 1963; 

Spreiter and Rizzi, 1974; Spreiter et al. , 1966, 
196S, 1970; Dryer and Faye-Peterson, 1966; Dryer 
and Heckman, 1967], provides a theoretical basis 
fcr the understanding and interpretation of 
observ'ations from the viewpoint of a fluid rather 
than a particle description of the flow. The 
basic pattern for the associated calculations 
using this model was established by Spreiter et 
al. [1966] for the case of the earth, or any other 
planet with a dominant dipole magnetic field; and 
by Spreiter et al. [1970] for the case of Venus 
and Mars, or any other planet with no significant 
cagnetic field but a sufficiently dense ionosphere 
tc stand off the solar wind above the planetary 
surface. Tnis consists of (1) representing the 
Ccgnetopsuse or loncpause by a magnetohydrodynamlc 
tangential discontinuity, (2) solving for the 
snape of the magnetospherlc or ionospheric 
cf Steele using simplified approximations for the 
pressure of the deflected solar wind flow et the 
obstacle boundary and for the planetary magnetic 
and atmospheric prope-.les, (3) solx'lng for the 
location of the bow ue.'e and for the density, 
velocity, pressure, and temperature throughout 
the flow as a gasdynamlc problem, disregarding 
the electromagnetic forces because of the charac- 
teristically high values of the free-stream 
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Alfvdn Mach number, and finally (4) calculating 
the magnetic field as a frozen-in convected field 
asBuming chat the other fluid properties are 
provided sufficiently accurately by the gas- 
dynamic solution. 

Prior to the work reported recently by Stahara 
et al. (1977} , the utility of this model was 
severely restricted because results were avail- 
able only in the form of plots in archival 
journals for a limited set of solar wind and 
planetary conditions. Results for intermediate 
points or for other conditions had to be deter- •' 
mined by interpolation. Stahara et al. [1977] 
improved on the basic model in many significant 
ways. New and more effective algorithms were 
introduced to perform the calculations, enabling 
the determination of results for a wide range of 
specific conditions typical of solar wind/ 
terrestrial planet interactions for both magnetic 
and nonmagnetic planets. 

In the current work suzaarlzed here, that model 
has been extended and generalized in several 
Important directions. These Include (1) extension 
of the gasdynamlc capability for treating free- 
stream Mach numbers as low as M^ • 2, (2) devel- 
opment of a new family of ionopause shapes which 
includes the effect of gravitational variation in 
scale height of the ionosphere, and (3) develop- 
ment of the capability for readily determin3.ng 
the plasma properties along an arbitrary space- 
craft trajectory, simultaneously accounting for 
an arbitrary direction for the Incident solar 
wind relative to a chosen coordinate system. 

The Mathematical Model - Formulation 
of the Fluid Representation 

The fundamental assumption underlying the 
present work and that reported In all of the 
references cited above is that many properties 
of solar wind flow around a planetary magneto/ 
ionosphere can be adequately described by the 
continuum equations of magnetohydrodynamics for a 
single-component perfect gas having infinite 
electrical conductivity and zero viscosity and 
thermal conductivity. 

Governing Equations 


The magnetohydrodynamlc equations for the 
conservation of mass, momencum, energy, and 
magnetic field are as follows: 
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flux are conserved. These are approxisiacions to 
comparatively thin surfaces across which similar 
but continuous changes in the fluid and magnetic 
properties occur in the corresponding theory of a 
dissipative gas, and they correspond physically 
to the bow wave, ionosphere boundary, and possibly 
other thin regions of rapidly changing properties. 
The conservation of mass, momentum, energy, and 
maguetic flux lead to the following relations 
between the quantities on the two s'.des of any 
such discontinuity: 
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and the equation of state of a perfect gas is 
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In these equations and those to follow, the 
s\nbols £, p, v, T, e ' CyT, and h • CpT refer 
to the density, pressure, velocity, temperature, 
internal energy, and enthalpy, and Cy and Cp 
refer to tne specific heats at constant volume 
and pressure_in the usual way. The symbol 
R = (Cp-Cy)K • 8.2*10^ erg/g°K is the universal 
gas constant, and K is the mean molecular weight 
nondimens ionalized so that M ■ 16 for_atomic 
oxygen. For fully ionized hydrogen, M is thus 
1/2. The magnetic field | and the Poyncing vector 
S for the flux of electromagnetic energy are 
expressed in terms of Gaussian units. The 
gravitational potential 6 and acceleration of 
gravity g are assumed to be due to massive fixed 
bodies so that their time derivatives are zero. 

G « 6. 67*10~®cm^/g s^ is the gravitational con- 
stant. These equations apply in the region 
exterior to the ionosphere boundary, as shown in 
Figure 1, and also in a degenerate sense in the 
ionosphere. 

Because of the omission of dissipative terms 
in these equations, surfaces of dlsconclnulty may 
develop in the solution, across which the fluid 
magnetic properties change abruptly, but in such 
a way chat mass, momentum, energy, and magnetic 


Here, (n, t) denote unit vectors normal and 
tangential to the discontinuity surface, where 
qjj represents the local normal velocity of the 
discontinuity surface and v* = Vj^ - qjj is the 
fluid normal velocity component relative to the 
normal velocity qjj of the discontinuity surface. 
The square brackets ere used to indicate the 
difference between the enclosed quantities on 
the two sides of the discontinuity, as in 
(Q) ■ Q2 - Qj , wnere subscripts 1 and 2 refer to 
the conditions on the upstream and downstream 
sides, respectively, of the discontinuity. 
Correspondingly, the average of and Q 2 is 
represented by <Q> - (Q^ + Q 2 )/ 2 . 

Five classes of discontinuities are described 
by (7)-(10). Of these, only two, the fast shock 
wave and the tangential discontinuity, are of 
significance in the present analysis. The former, 
having properties that satisfy the relations 
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Fig. 1. Illustration of overall features of solar wind interaction with nonmagnetic 
terrestrial planets. 
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is appropriate for representing the bow wave; and 
the latter, having properties given by 


V* - - 0 [vj i 0 IBJ i 0 

[; ] y 0 [p + B^/8t] - 0 


( 12 ) 


is appropriate for representing the lonopause 
that separates the flowing solar wind and the 
planetary ionosphere. 

High Alfv^n Mach Kunber Approximation 


Two ieportant parameters that characterize 
steady magnetohydrodynamic flow past a stationary 
obstacle are the free-stream Mach number 

“ v^/a^; and the Alfvdn Math number 
in which a « (Yp/p)^^^ is the speed of sound and 
A ' (B^/drt)^' 2 speed of a rotational or 

A.lfvdn wave along the direction of the magnetic 
field, and y « Cp/Cy is the ratio of specific 
heats of tne plasma. 

For typical solar wind conditions at Venus, 
both l-i„ and M;^^ are high, values of the order of 
10 being representative. As a result, important 
simplifications of the equations may be intro- 
duced. Since the order of magnitude of the 
momentum flux term in (2) and the kinetic energy 
flux term in (3) is related to that of the 
magnetic terms in the same equations by the 
square of M/^_, consideration of the magnetic 
terms may be decoupled from the solution of the 
gasdynamic portions of those equations. The equa- 
tions for the fluid motion thereby reduce approxi- 
mately to those of gasdynamics , and the magnetic 
field E can be determined subsequently by solving 
the remaining equations using the values for o 
and v so determined. Furthermore, because of the 
large value for M^, the terms involving | and d 
can be disregarded in the solar wind flow because 
of the relatively small effects of gravity on the 
flow. This is not true, of course, in the iono- 
sphere, where the comparative stationarity of the 
gases results in a very nearly hydrostatic 
variation of the pressure with altitude. 

The precise specification of the range of use- 
fulness of the decoupled approximation remains 
elusive in the absence of exact solutions of the 
complete magnetohydrodynamic equations. The only 
such solutions for the present category of flows 
are those of Sprelter and Rlzzi [197t] for the 
special case in which B_ is aligned with v^. 

Tney indicate results virtually identical to those 
of the simpler decoupled theory when M^ is about 
10 or greater and notably differing results as 
K. is decreased below 5. 

Aar 

For otner magnetic field alignments one can 
only make estimates of the relative magnitudes 
of the divergences of the neglected magnetic 
terms B^B^/dr + Sfk compared tt the 

terms cv^vj. + piik PVy(v^/2 + e + p/e) 
retained in the gasdynamic analysis and then 
speculate about how large the magnetic terms can 
become before they invalidate the solution. In 
regions where Che momentum flux and kinetic 
energy terms are large in comparison with the 
pressure and enthalpy terms, the souare of the 
local Alfvdn Mach number » (pv^/2) * (B^/8r.) 


may be considered to provide a rough measure of 
the dominance of the gasdynamic terms in the 
calculation of the fluid motion. ^Since this 
quantity is related to M^ by • (o/c )(v/v )*. 
(B^/B)*M;,^, and p/p„, v/v^, and B/B_^ are indepen- 
dent of Mj>^ at large My^^, it follows that M^ 
tends to remain large for large M;^ except »^ere 
pv* « p^v^ or where B^ >> B^. These conditions 
are most likely to be encountered near the 
stagnation point, where v vanishes; but p reaches 
a maximum value comparable to p v^ there, and 
that contravenes the initial premise upon which 
the estimate is predicated. A localized region 
of small M^ does not therefore necessarily signal 
a failure of the approximation. 

As a further consideration in making such 
estimates, it should be noted that the orienta- 
tion of B^ is very Important in the determination 
of B'/B^ throughout the flow. In particular, 
near the stagnation point, B'/Bb >> 1 if B,, makes 
a substantial angle with v ; whereas B^/B^ 
vanishes at the stagnation point when B is 
parallel to y„, since B is everywhere proportional 
to pv [Sprelter and Rlzzi, 1974). For this case 
therefore ■ (c„/p)-'i^, and does not vanish 
even at the stagnation point because p/p„ is 
about 4-5 for representative conditions. 

A definitive theoretical evaluation of the 
decoupled approximation will have to await the 
development of more advanced methods for the 
solution of the complete magnetohydrodynamic 
equations. In the meanwhile, considerable in- 
sight can be gained by comparing the predictions 
with conditions actually measured in space during 
appropriately steady conditions. The small scale 
of the Venus flow field compared with that of the 
earth's magnetosphere makes the data from Pioneer 
Venus particularly valuable for this purpose; and 
the results presented here show our initial 
efforts to carry out such comparisons. This 
should be considered no more than a beginning, 
however, as we only have data from two orbits; 
many additional comparisons should be carried 
out to obtain a more comprehensive understanding 
of the relation between the theoretical and 
observational results. 

Under the conditions outlined above, (2), (3), 
(6), and (9) are thus simplified to the following: 
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These equations, together with (1) and (7), are 
just those of gtsdynamics and must be solved for 
the location of the bow wave and the gasdynamic 
properties v, p, c, and T of the flow field about 
the ionospheric obstacle. With these properties 
so determined, (4) and (10) can be used to deter- 
mine the magnetic field B. Although all of the 
solutions described here are for the steady state 
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scale heights E/R^ and (2) gravitational!)’ varying scale heights associated vlch a 
constant I/K atmosphere and characterized by H/R^. 


obtained by setting 3/3t ■ 0, and i.e., 

qj, • 0, the unsteady equations are presented In 
the foregoing because one of the computational 
methods used to determine the gasdynamic solution 
employs an unsteady procedure, Integrating In 
time until the steady state solution Is obtained 
asymptotically. 

Determination of the Ionosphere Boundary 

The determination of the Ionosphere boundary, 
or lonopause, makes use of the Idealization that 
the ionospheric plasma Is spherically symmetric, 
hydrostatically supported, Infinitely electrically 
conducting, devoid of any magnetic field, gravi- 
tationally bound to the planet, and Incapable of 
mixing with the solar vlnd plasma, as sketched 
In Figure 1. This Interior plasma is considered 
to be separated from the flovlng solar wind 
plasma by a tangential discontinuity across 
which the relations of (12) apply. Because the 
magnetic pressure B^/6ir Is much smaller than the 
ionized gas pressure p In the lonopause, the 
relation [p + B^/Sn] « 0 of (12) may be 
simplified to 


p “ (p + B^/Se) (17 

*^atm flow 

The variation of the pressure with radial 
distance from the planet center is given accord- 
ing to the hydrostatic equation by 


(18) 


where p and c are the ionized gas pressure and 
density, r is the radial distance measured from 
the center of the planet, and g“ lg| ■ 6s(rg/r)*, 
where subscript s denotes values at the surface 
of the planet. Combination of (18) and the 
perfect gas law given by (6) provides 


P " Pr exp 



(19) 


where pj^ ls_the pressure at some reference radius 
Rr and H - RT/Mg represents the local scale 
height of the ionospheric plasma. 

If H is regarded as constant, i.e., if 
variations of g and T with r are disregarded, 

(19) can be Integrated directly to obtain 


P " Pj^ exp 



( 20 ) 


Because of limited knowledge of the Ionospheric 
properties of Venus at the time, the simple 
variation of p with r given by (20) was adopted 
In the Initial study of Spreiter et al. [1970] 
and also In the more recent analysis of Stahara 
et al. [1977] involving the initial applicf.ion 
of advanced computational methods to this problem. 
Data from the Pioneer Venus spacecraft [Knudsen 
et al., 1979a, b] have provided much additional 
knowledge of conditions in the Venus ionosphere 
and, in particular, have shown chat the assump- 
tion of an Isothermal (T • const) ionosphere at 
the heights of concern here Is quite reasonable. 
With this assumption, including the variation of 


cp/dr - -pg 
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gravity with height, the following result for the 
pressure is obtained: 


P 






(21a) 


where 


« ‘ «s • 


(21b) 


Kf being the planetary radius and Kg « RT/Mgg. 
Equations (20) and (21) provide the two nodels 
enployed in this study for the ionosphere pres- 
sure variation required in the pressure balance 
relation of (17). 

On the basis of the typically hypersonic 
values for it has been assui&ed in all the 
nrevious analyses cited above that the pressure 
cf the solar wind plasma on the lonopause could 
D£ represented adequately by the Newtonian 
approximation 

2 2 

p • Kc V cos t 

^ ec at 

_ 1 _ ( 22 ) 
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where t is the angle between the outward normal to 
the ionopause and the flow direction of the undis- 
turbed solar wind incident on the bow wave and K 
is a constant usually taken as one but whose 
actual value to provide the correct gasdynatnic 
stagnation pressure at the nose of the ionopause 
IS as indicated. 

The Newtonian approximation stems from gas- 
dynamic analysis of hypersonic flows and is 
aspropriate in the range of it from 0® to somewhat 
less than 90®. In view of the close association 
between the magnetic pressure B‘/8 t and the gas 
pressure p in magnetohydrodynamics, it seems 
equally plausible to assume that the corresponding 
relation 


p + B^/St - Kp_v^cos^i). (23) 

is-an appropriate, if not_6uperior,_ relation for 
the combined magnetic and gas pressure of solar 
wind plasma on a planetary ionopause or magneto- 
pause . 

For the high Kach number flows typical of solar 
wind conditions, K approaches 0.681 for T » 5/3 
and 0.6-4 for 7 « 2. A discussion of effects of 
r.incr constituents, such as ionized helium, in the 
sclar wind and of some differences between a fluid 
enc collisior.less particle representation has been 
given by Spreiter (1976), but the most important 
effect of the introduction of the Newtonian 
aporoximatlon is that the calculation of the shape 
of the ionopause decouples from the calculation of 
the properties of the external flow. In this way 
w€ arrive at the following equation for the pres- 
sure balance at the ionopause location R^: 

Rc v^cos^v • p_/.(R. ) (24) 

« K 1 


or 


A(R^) 




g,T « const 



(25a) 


(25b) 


depending upon whether or not the gravitation 
variation is included in the scale height. It is 
convenient to choose the reference radius Rp to be 
the distance R^ from the center of the planet to 
the nose of the ionopause. Then pj{ • Pp « Kx> v^, 
and “ “ 


cos 


A(R^) 


(26) 


at all points along the ionopause. To proceed, 
cos^ij' must be expressed in terms of Rj and £, where 
Ri(£) represents the radial coordinate of the 
ionopause and £ is the angle measured at the center 
of the planet with respect to a line that extends 
directly upstream. The appropriate expression is 
(Spreiter et al., 1970) 


2 

cos U 



(R^dScosB + dR^sin£)‘ 

dR^ + (R.dS)^ 

1 1 


(27) 


Substitution of this expression into (26) results 
in the following ordinary differential equation 
for the ordinates of the ionosphere boundary: 




Fig. 3. Comparison of former and present compu- 
tational procedures for determining the gas- 
dynamic flow properties of solar wind-magneto/ 
ionopause interactions. 
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Results for various ionopause shapes obtained 
by Integrating (28) for different values of H/Rq 
using the constant scale height Bodel equation 
(25a) were provided by Sprelter et al. [1970). 
Similar results using the Isothermal model equa- 
tion (25b) In the range 0.01 £. H/R^ S. 0.5 are 
provided in Figure 2, where for comparison the 
constant scale height shapes for corresponding 
K/P.Q values are also illustrated. Data from 
Pioneer Venus indicate that the range of Interest 
for application to Venus appears to be 
0.01 s. H/Rj, <. 0.10. They also indicate that the 
theory provides a good representation of the 
ionopause shape on the dayside but chat conditions 
on the nighcslde are both more complex and variable 
than presently accounted for by the theory (see, 
for example. Brace et al. [1979], Russell et al. 
[1979b], or Brace et al. (this issue). To provide 
insight into some of the changes to be encountered 
rn the theory, some results are also included 
herein in which the nigntside ionopause Is 
arbitrarily flared in and out from che calculated 
shapes of Figures 1 and 2. 


Calculation of the Casdynamlc 
Flow Properties 


Determination of the gasdynamlc flow properties 
is the most difficult and time-consuming portion 
of the numerical solution of the foregoing equa- 
tions. As sketched in Figure 3, two different 
methods are used In the calculations: one for the 

nose region, where both subsonic and supersonic 
flows occur, and another downstream of that 
region, where a more computationally economical 
procedure can be employed, since che flow is 
everywhere supersonic. In che original analysis 
of Sprelter et al. [1970] an Inverse iteration 
method was used for che nose region, and Che 
method of characteristics was employed for che 
remaining supersonic region. Both of these 
procedures are inferior to more recently developed 
methods and have been superseded in che current 
model. The nose region is treated using a new 
axisymaetrlc implicit unsteady Euler equation 
solver specifically developed for the present 
application. That procedure determines the 
steady state solution in the nose region by an 
asymptotic time-marching procedure. The down- 
stream solution is determined by a shock capturing 
marching procedure which spatially advances the 
solution downstream as far as required by solving 
the steady Euler equations. Inasmuch as these 
methods are complex and lengthy to describe, only 
an outline of the procedures is provided here. 


The Nose Region Solution 


The nose region procedure solves che unsteady 
gasdynamlc equations (1), (6), (7), and (13)-(16) 
for axisyiaetrlc flow. The analysis commences by 
mapping the portion of che lonosheach bounded by 
the stagnation streamline, the bow wave, the 
ionopause, and a line through the terminator 
(E ■ ^/2) into a rectangle. The physical boundary 
conditions are also transferred appropriately. 
Initial flow field conditions are prescribed by 


using a linear variation of the flow properties 
between those calculated at an assumed location 
for che bow wave and those calculated at the iono— 
pause, using a Newtonian pressure distribution 
approximation on che latter boundary. Starting 
''Ich this initial flow, the unsteady equations ere 
integrated in a time-asymptotic fashion until the 
steady state solution is obtained. 

The basic numerical algorithm employed was 
developed by Beam and Warming [1976] and is second- 
order accurate, noniterative, and spatially 
factored. At che ionopause and bow shock bound- 
aries a modification of che differencing algorithm 
is introduced to account for the particular 
features of the present application. The cangency 
condition at Che ionopause is created using the 
method of Kencaer [1970], and the sharp discon- 
tinuity approach of Thomas et al. [1972] is used 
at che bow wave. Any secondary shocks which may 
occur in the nose region are automatically 
determined via the shock-capturing characteristics 
of the algorithm. 

In our earlier work [Stahara, et al., 1977) it 
was found that for certain ionopause shapes which 
have a significant amount of lateral flaring at 
the dawn-dusk terminator, for example, constant 
scale height shapes for H/Rj, i. 0.5, and/or cases 
involving low free-stream Mach numbers M„ d. 3, 
the axial component of velocity at some points on 
the terminator plane B “ ti/2 may become subsonic. 
Although this has no effect on the nose region 
solver, the downstream solution cannot be obtained 
for these cases, since the marching- region solver 
which determines the solution downstream of this 
starting plane requires supersonic axial veloc- 
ities in order to proceed. In the present effort 
this limitation has been removed by developing the 
capability for adding an additional portion of the 
flow field, located downstream of the terminator, 
to the blunt body solution as Illustrated in 
Figure 6. This effectively generalizes che 
capability of the present procedures to treat a 
wide variety of ionopause shapes. Including all 
of che shapes described by the constant scale 
height and constant temperature models found using 
(25a) and (25b), as well as to treat lower free- 
stream Mach numbers down to about M >2. 

The Downstream Region Solution 

The solver for the downstream supersonic region 
is based on che original work of Ruder et al. 

[1972, 1973] and Chaussee et al. [1975] and solves 
for axlsysmiecric flow the steady state gasdynamic 
equations in conservation law form indicated by 
(1) , (6), (7), and (13)- (16) with the time 
derivatives set to zero. Of these, (1^), which 
represents the conservation of energy, can be 
Integrated for the steady flow to yield the 
following equation for the total enthalpy: 

h^ • h + (u^ + v^)/2 “ const (29) 

where h • e + p/o • CpT is che enthalpy per unit 
mass. This relation also holds in the nose region 
in the limit as che steady state flow solution is 
approached. The computational mesh defined by 
lines of constant X and (R - Rb)/(Rs - R^), where 

and R^ are functions of X that describe the 
radial cylindrical coordinates of the bow shock 
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Fig. 4. Illustration of additional flow field segment to the obstacle nose -solution 
for determining the gasdynamic flow properties of solar wlnd-ionopause interactions. 


wave and ionopause at the same X as the field 
point (X, R) , is transformed into a rectangular 
computational grid. The gasdynamic equations are 
similarly transformed and Integrated numerically 
to obtain the solution. Since only the bow shock 
wave is treated as a sharp discontinuity and any 
others that may be present are 'captured' by the 
difference algorithm, selection of the appropriate 
finite difference scheme to advance the calcula- 
tion in the X direction is of prime importance. 
Following the analysis of Kutler et al . [1972, 
1973] and Cnaussee et al. [1973] the method used 
is the finite difference predictor-corrector 
scheme of KacCormack [1969], the most efficient 
second-order algorithm for shock-capturing 
calculations. 

Calculation of the Streamlines 


and the streamline determination is made by step- 
wise integration in X using a modified third-order 
Euler predictor-corrector method. Bivariate 
linear interpolation between the flow field grid 
points IS employed to obtain the velocity com- 
ponents (u, v) required at the stepwise points 
along the streamline. 


Calculation of the Magnetic Field 


With the flow properties known from the gas- 
dynamic calculations, the magnetic field B may be 
determined by integrating (4) and (10); or 
alternatively, the following equations may be 
derived from them: 


_2. 

Dt 


JJ^B-dS-O 


_D_ 

Dt 


I) 


- (B'V)v 

p - - - 


(32) 


The accurate determination of the gasdynamic 
streamlines is crucial, since they provide the 
basis for the frozen magnetic field calculation 
described in the following section. The stream- 
lines are determined by integrating fluid particle 

-trajectories-through -the-known velocity-field, 

since this procedure was found to be more accurate 
than the alternative mass flow calculation. The 
calculation of a particular streamline is 
initiated at the point where the streamline 
crosses the bow shock. At that point, exact 
values of the streamline slope dRg/dX are given 
by the following relation involving M^, T, and the 
angle dg between the free-stream direction and the 
tangent to the shock wave: 

dR- (2cotd_)(M^sin^6. - 1) 

— I . i 1 (30) 

2 + K (X + 1 - 2sin 6-) 

* a 

which is contained implicitly in both nose and 
marching region gasdynamic solutions. At other 
points, the local streamline slope is given by 
the ratio 


dRj/dX “ v/u 


(31) 


in which S is an arbitrary surface moving with the 
fluid and D/Dt is the material or substantial 
derivative dD/Dt » i/it + (vV). These equations 
are commonly interpreted as indicating the field 
lines move with the fluid. For the steady state 
in which’ 9/9 1 ^ “ 0 Md’vJJ » Vj^’, th^^ equations 
lead to a straightforward calculation in which the 
vector distance from each point on an arbitrarily 
selected field line to its corresponding point on 
an adjacent field line in the downstream direction 
is determined by numerically integrating J v dt 
over a fixed time interval At as indicated in 
Figure 5. Once the magnetic field lines are 



I'tor* 


lonopouSf 


Fig. 5. Illustration of quantities used for mag- 
netic field line calculation in the plane of 
magnetic symmetry. 
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determined, Che magnetic field at any point may be 
calculated from the relation 


I H 


where Lt is the vector length of a small element 
of a flux tube. 

Such a procedure is valid generally, but its 
use is confined in the present calculation to only 
the component of the magnetic field (B)j. chat lies 
in the plane of magnetic symmetry defined by the 
plane containing Che axis of symmetry of the iono- 
pause and the magnetic field lines upstream of the 
bow wave and furthermore is associated with only 
Che component of the interplanetary magnetic field 
that lies in the plane of magnetic symmetry and xs 
perpendicular to v_^. The remainder of the magnetic 
field calculation makes use of a decomposition due 
to Alksne and Webster fl970] in which the 
axisymmetrlc properties of Che gasdynamic solution 
and the linearity of the magnetic field equations 
(i) and (10) or (32) are employed to derive the 
following relationship for the magnetic field Bp 
at any point P: 


?P 





(i4) 


As illustrated in Figure 6, subscripts ", and n 
refer to contributions associated with the com- 
ponent E." of B^ parallel to y^; the component B«x 
perpendicular to v^ in the plane that contains the 
point P, the center of the planet, and the vector 
v„; and the component B„^ normal to the latter 
plane, where Cj, is a unit vector in the latter 
direction. The unit ratios (Bp/B„).t and (Bp/B^)jj 
can be calculated directly from the gasdynamic 
solution by the expressions 



P-P 


V 



Vp 


(35) 



three-dimensional magnetic field. 


where Rp is the radial cylindrical coordinate of 
the streamline through P, as indicated in Figure 6. 

In carrying out the determination of (Bp/B„)_^ 
using (33), values for are determined 

initially at the points where the streamlines and 
perpendicular-component field lines intersect. 

A generalized quadrilateral interpolation scheme 
followed by a fifth-order smoothing is then 
employed to determine the corresponding values at 
the computational grid points where .values for 
p/p„ are available for calculation of (Bj,/B„)j.. 

At the bow shock, exact formulas are used for the 
magnitude and direction of tZ: 


- 1 + cot^e(l+D^) 

- 2D X CSC 6 X cot 6 X cos(6-{) 

!(/ • e + sin ^ { [D X cot e 

X sln(6-6)]/(lA£|/|i4_i)) 

where 

- 1 - 4(Mfsin^e - 1) 

X (YM^sin^e + 1)/[(Y+1) Vsin^e] 

coto = ^(Y+l)K^/t2(«fsin^S- 1)) - ij tan6 



(36) 


(37) 


and is the angle between the free-stream direc- 
tion and the hZ vector. 

Finally, the resultant magnetic field vector 
can be expressed in terms of components relative 
to any orthogonal (X, Y, Z) coordinate system. 

If the coordinate system is oriented so chat the 
point P lies in the (X, Y) plane, the magnetic 
components are 


By/B • I(|b|/B )„ cos 6 cos o 

+ (|b|/B^)j^ cos sin Q^)]cos 

Bv/B “ K|b|/B )„ sin p cos o (38) 

+ (|b|/B ), sin p sin o ))cos o 
- ® - p n 

B,/B » (E/B ) sin o 
Z "x “ n n 


where e is the local flow angle given by 


t 


tan 


-1 


V 

u 


(39) 


and the interplanetary magnetic field angles 
and 0 ^ indicated in Figure 6 are defined by 


o 

P 




I 

a 

>- 


• tan 

[‘aJ 


(AO) 
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Fig. 7. Illustration of sun-planet (Xg, y^, rg) and solar wind (x, y, z) coordinate 
systeas and the aziauthal (fO and polar (tp) solar wind angles, both shown in a 
positive sense. 


Q 

n 





« tan 


-1 


/(Bj^ )^ + (By 


(41) 


The generallratlon of these results when the point 
P is at some arbitrary (X, T, Z) point not in the 
(X, Y) plane are provided below in the spacecraft 
trajectory section. 


Calculation of the Contour Lines 


In order to provide convenient summary Informa- 
tion of the detailed plassia and magnetic field 
properties determined by the above model, contour 
and field line maps are determined for the impor- 
tant field properties. Contours are calculated 
for nondlmensionalized velocity |y|/v_, density 
o/p„, and magnetic field components (|Bi/B„)>', 

ITetails of these calcula- 
tions are provided by Stahara et al. [1980). 

Since the temperature is a function of lyl/v^ only 
for a specified and Y, l.e., 


T/T 


1 + 


Y-lt 2 

1 

[IyI] 

2' 

l~J 

: * 

1''- j 



(42) 



Fig. 8. Illustration of solar wind (x, y, z) and (X, Y, 2) coordinate systems and 
the interplanetary magnetic field and magnetic field angles ( 0 ^, 0 ^). 
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- — Fig. 9i CaIculated_lonopause and^bow shock locations and ionosheath properties for 
■ 8.0, Y « 5/3, H/Rq - 0.10, and either parallel or perpendicular to v^. 

(a) Streamline map. (b) Velocity map, (c) Density map, (d) Temperature map. 

(e) Parallel magnetic field map. (f) Perpendicular magnetic field map. 


velocity contours may also be considered as tem- 
perature contours with only a relabeling required. 

Solar Ecllptlc/Solar Wind Coordinate 
Transformation 

To compare results of the present calculations 
with observational data, it is necessary to con- 
sider the appropriate transformations between the 
spacecraft and solar wind coordinate systems. 

Part of the data required as input to the theo- 
retical model consists of oncoming interplanetary 
values of solar wind temperature, density, and 
velocity and magnetic field vector components. 
Tnese are usually reported in a sun-planet or 


solar ecliptic reference frame, whereas the 
natural coordinate system for the theoretical 
model is one which aligns the axial direction with 
the oncoming solar wind, since the gasdynamic 
calculation is assumed to be axisymmetrlc about 
this direction. Thus the interplanetary input 
data must be transformed to the solar wind system 
to initiate the theoretical determination. Once 
the gasdynamic and magnetic field calculations in 
the solar wind system are complete, those results 
must then be transformed back to the aun-planet 
system to allow direct comparison with the space- 
craft data obtained as a function of time along 
its trajectory. 

For the sun-planet system (x^ , yg, Zg) the 
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Fig. 9. (continued) 


origin is placed at the planetary center, the 
Xg axis Is directed toward the sun, the y^ axis 
is directed opposite to the planetary orbital 
notion, and the Zg axis is orthogonally northward. 
Tne direction of the oncoming solar wind is such 
that the azimuthal angle, including aberration 
effects due to planetary motion, of the solar wind 
velocity vector in the planetary orbital plane is 
and the out-of^orbltal plane polar angle is »p. 
The positive sense of the azimuthal angle is for 
west-to-east flow, and that of the polar angle is 
for north-to-south flow, as indicated in Figure 7. 
Also iilustreted is the solar wind (x, y, *) 
coordinate system defined by (fi, Op) . This system 
is somewhat Inconvenient because the direction of 
solar wind flow le in the negative x direction; 
therefore the internal gasdynamic and magnetic 
field calculations arc performed in the (X, Y, Z) 
system illustrated in Figure 6. 

Tne coordinates and vector transformations from 


the ecliptic sun-planet (Xg, yg, Zg) system to the 
(X, Y, Z) solar wind syscem are given by 


Q ^ 


-coBflcosi -slnflcose sind 

P P P 


Qxs' 

'Jy 

m 

sinC -cosfl 0 


'^Ys 



-cosHsind sinflsind cosd 

P P P 





J 

J 


where (Qy, Qy, Q^) represents the Cartesian com- 
ponents of any vector referred to the (X, Y, Z) 


system, and (Qxg, Oj^, Qz,) represents the 
corresponding vectcr in the (xg, yg, Zg) system. 
Thus for a transformatlotj of coordinates. 


(Q,. Q„. QJ - (X, Y, Z) 

A y z 

(0^ . 0^^ . Q, ) - (*,. y^. r^) 

t ^6 6 


( 44 ) 
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(f> 

Fig. 9. (continued) 


while for a vector tranaformatlon of, aay, the 
aagnetic field. 


<Q... Qy. Q.) ■ (B,.. 2,) 


- Qv - 'Jz ^ - ®y - ^ ) 

S ' S E S S S 


(«5) 


Tne inverse t tans format Ion fron the solar wind to 
tne sun-planet systeo is given by 




cosTdCoe^ 

P 

sin^ 

-cosnsiod 

P 

-cosfl 

sinfisin^ 

sine 

P 

0 

cos^p 


( 46 ) 


Properties Along a Spacecraft Trajectory 

One of the primary alas of the present effort 
has been to develop the capability to determine 
plasma and magnetic field properties, as predicted 
by the theoretical model, at locations specified 
along an arbitrary spacecraft trajectory, and in 
such a form as to enable comparisons to be made 
directly with actual spacecraft data. To 
accomplish this, the following procedure has been 
developed and implemented in the current model. 
First, from the knows oncoming Interplanetary 
conditions provided in a sun-planet reference 
frame, the azimuthal and polar solar wind angles 
(n, «p) are employed in (43) to establish both the 
location (Xp, Yp, 2p) of the trajectory point and 
the interplanetary magnetic field components 

By^, Bj^) in Che solar wind (X, Y, 2) frame. 
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M - 3.0, y - 5/3. 


In the next step, In which the axisymnietric gas- 
dynasic and unit aagnetlc field calculations are 
carried out, it is convenient to take advantage of 
the axisytsnetric property of the gasdynaaiic flow 
and introduce still another coordinate systen 
(X, y', r') related to the (X, Y, Z) systco by a 
rotation about the X axis by an angle 6 such that 
tne y' axis passes through the point P. This 
rotation defines a new coordinate systen (x', y', 
s'), where 


x'' 


1 0 o' 


x' 

y ' 

m 

0 cose sine 


Y 

z’ 


0 -sine cos6 


Z 



. 


, . 


in which 

e « tan 

and 


Z 

_E 

'P 


(A8) 


X 

P 


Thus the (x', y') plane which contains the X axis 
and the arbitrary point P corresp onds directly to 
the plane (X, R) » (Xp, A ^ + Z|) , in which the 
axisytmietric gasdynanic flow properties are calcu- 
lated. In particular, the velocity magnitude v, 
density p, and flow angle b at the point P are 
found by bilinear interpolation through the (X, R) 
flow field grid. The vector velocity in the (X, 

Y, Z) system is then given by the transformation 


fv 

X 

''y 

m 

1 0 
0 cose 

0 

-sine 


v cosb 
v sinb 
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0 sine 

cose 


0 


(50) 


and then in the solar ecliptic system by the 
transformation indicated by 


r "s 

V 
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Calculation of the magnetic field at an 
arbitrary point is somewhat more complicated, 
since these components are dependent upon the 
orientation of the incident interplanetaty mag- 
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wind directions; also, various bow shock shapes for different Interplanetary solar 
wind conditions. 




Fig. 12. Ccnpsrlscu of observed (OPA) and theoretical time histories of lonosheath 
plasma properties for Pioneer Venus orbit 6 based on inbound and outbound inter- 
planetary solar wind conditions using a gasdynamic solution for ■ 13.3, T « 2.0. 
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Interslanetary Conditions 



—o— Inbound 13.3 J.o 6.5® -1.4® 3.99 -6.94 -2.69 



Tinr from periapsiE, min. 

(a) 




-60 -40 -20 C JO 40 60 

Tliaa troB parlapaia, mm 
(b) 

Fig. 13. Comparison of observed (OHiC) and theoretical time histories for the 
magnitude of the magnetic field for Pioneer Venus orbit 6 based on inbound and 
outbound interplanetary conditions using gasdynanlc solution for “ 13.3, If ' 
(a) Magnetic field magnitude. (b) Magnetic field components. 
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wind directions; also, various bow shock shapes for different interplanetary solar 
wind conditions. 


netic field. With the known (By , By , Bz ) 
components the corresponding components (Bx„, 
Ey^, in the rotated system are given by 


B’ 1 


1 o' o' 


\k ] 
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( 52 ) 


In this reference frame the parallel, perpendic- 
ular, and normal interplanetary components are 
identified as 
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It 




( 53 ) 
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The magnetic field angles 
system are given by 


o' and o' in the rotated 
p n 
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( 55 ) 


The magnetic angle i|i associated with the incident 
perpendicular component and the unit magnetic 
field ratios (|b|/B_,)„, (|b|/B^)^, and (B/Bjj, in 
the rotated system are next determined by bilinear 
interpolation through the flow field grid. Then 
the magnetic field components (B^, By, Bj) in the 
rotated system are calculated from 
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( 57 ) 
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Fig. 15. ComparieoTi of observed and theoretical tlae histories of ionosheath plasma 
properties for Pioneer Venus orbit 3 based on inbound and outbound interplanetary 
solar wind conditions. 
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The magnetic field components in the solar wind 
(X, y, Z) system are then determined from the 
rotational transformation 


with B^ either parallel or perpendicular to y^ . 
These and similar determinations for the following 
35 other sets of conditions have been included in 
a catalog of results evailable from Stahara et al. 
[I960): 

» (2.0, 3.0, 5.0, B.O, 12.0, 25.0) 
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and finally in the sun-planet system from 
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-cosOcosd slnO -cosnsine 
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(60) 


Results 


Using the computational procedures described 
above, a large number of solutions have been 
calculated for a'wlde variety of conditions repre- 
sentative cf those that might be experienced at 
Venus. A sample of these results is illustrated 
in Figure 9 for • 8, f - 5/3, and H/Rq • 0.10 


B/R - (0.01, 0.10, 0.25) 
o 

i/R - (0.10, 0.20, 0.25) 
o 

Also Included In that reference are a few results 
for conditions well outside the range likely to be 
encountered at Venus but that are useful because 
they illustrate the extended capability of the 
present procedures. 

Verification of the new procedures has been 
done in a variety of ways in addition to compari- 
son with previous results for sets of conditions 
for which both sets of results could be determined. 
These include consideration of a variety of 
special test cases in which the location in the 
flow field and the incident interplanetary r ignetic 
field orientation were systematicall> charged so 
as to produce both sysssetric and antlsytmietrlc 
changes in the resultant ionosheath magnetic 
field, as well as to reverse the roles of the 
perpendicular and normal components. All of these 
various permutations of the magnetic field calcu- 
lations were successfully verified. 

The final and ultimate check of the theoretical 
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Interplanetary Conditions 


Theory 


Inbound 

Outbound 


tt, y n dp 

7.38 2.0 3.3® 0.X5® -1.78 -7.86 0.08 

5.»6 2.0 3.7® 4.90® «. 38 -6.45 -0. 3i 



(a) 

Fig. 16. Coaparlson of observed (0MA.G) and theoretical tine histories for the magnetic 
field for Pioneer Venus orbit 3 based on Inbound and outbound interplanetary solar wind 
conditions using gasdynanic solutions ■ 7.38, Y « 2.0 for inbound and « 5.96, 

Y « 2.0 for outbound calculations. (a) Hagnetlc field magnitudes, (b) Magnetic field 
components. 


model and numerical procedures lies in the com- 
parison of the predicted results with data 
actually measured by a spacecraft. Attention is 
now turned, accordingly, to the results of pre- 
liminary comparisons with data obtained from 
orbits 3 and 6 of the Pioneer-Venus orblter. 

These orbits were selected in conjunction with the 
experimenters on the basis of completeness of data 
and apparent steadiness of conditions during the 
time of interest for the study. 

To perform calculations suitable for comparison 
with data from a particular orbit, it is necessary 
to specify appropriate conditions corresponding to 
those actually measured in both the ionosphere and 
in the oncoming interplanetary solar wind plasma. 

On the basis of data -from the orblter retarding 
potential analyzer (OPRA) indicating the iono- 
spheric scale height to be approximately 200 km 
for orbits 3 and 6, a value of 0.03 was selected 
for H/Rq or H/R^,. For such small values the two 
ionospheric pressure models indicated by (20) and 
(21) yield essentially the same obstacle shape, 
as can be seen from Figure 2. With regard to oo- 
cotting interplanetary conditions it is required 
that values for the solar wind bulk velocity v^, 
density p^, temperature T^, and magnetic field 
be specified. The first three are provided by 
the orbiter plasma analyzer (OPA) , and the 
magnetic field is given by the orbiter flux gate 
magnetometer (OMAC) . OPA provides either ion or 
electron data, but not both simultaneously. For 
orbits 3 and 6, ion measurements were available 
and have been employed. Data on the oncoming 
direction of the solar wind, as given in terms of 
the angles (G, »p) , define the coordinate rotations 
required to align the gasdynamic calculation in Che 


free-stream solar wind direction; while information 
of solar wind speed, density, and temperature serve 
to define the free-stream gasdynamic Mach number 
M . 

Although it is known the temperatures of Che 
ions and electrons may be substantially different 
in the solar wind plasma (a ratio le/T^ - 4 is 
often representative), they are assumed to be equal 
(Te/Ti • 1) in the derivation of the single-fluid 
magnetohydrodynamic equations upon which the 
present calculations are based. Now that the 
theory has been put into a more readily usable 
form and more extensive data on the temperature 
of both electrons and ions are becoming available 
for comparative studies, steps should be taken to 
improve understanding of the relation between 
single- and multiple-component plasma theories for 
space plasma flows, and the calculative procedures 
refined as appropriate. 

With this information the detailed gasdynamic 
and unit magnetic field calculations in the iono- 
sheach region can be carried out. To provide an 
idea of the detail obtained by the present com- 
putational procedures, Figure 10 displays the flov 
field grid for one of the gasdynamic solutions 
used in the comparisons. Values for each of the 
plasma properties v/v^, d/p_, T/T„, (B/B^).., 

(B/B_)i and (B/B^)jj are determined at each inter- 
section of the grid lines, including the bow 
shock, stagnation streamline, and lonopause 
boundary. The final output of the calculation 
consists of these quantities and the coordinates 
of the bow shock and lonopause, expressed either 
in tabular form or as plotted contours, and also 
as time histories along a specified trajectory. 

For comparisons with spacecraft data the time- 
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Fig. 16. (continued) 


history predictions are the most convenient. For 
the theoretical predictions, evo sets of results 
are usually generated on the basis of the last/ 
first interplanetary solar vind properties 
measured before/after bov shock Inbound /outbound 
crossing. 

Figure 11 displays some overall flow field 
results for orbit 6. Indicated are bow shock 
locations for the three combinations of and 
plasma specific heat ratio t, i.e., y) ■ 

(13.3, 5/3), (13.3j_2), (3.0, 5/3) for flow about 
an ionopause with E/P.^ « 0.03. Also indicated are 
two sets of points (solid lines with circles, and 
dashed lines with squares) representing the space- 
craft trajectory as viewed in two solar wind 
oriented coordinate systems. The trajectory 
indicated by the solid lines and circles is that 
based on the last measured direction (0, ^) - 
(6.5', -l.i') of the interplanetary soler wind 
Just prior to crossing the bov shock on the in- 
bound leg, while the dashed lines and squares 
denote the trajectory based on the first measured 
direction (f), ftp) » (A. 9', 7.6*) of the solar wind 
immediately after crossing the bow shock on the 
outbound leg. These results illustrate the 
extremely large dependence of spatial position of 


a trajectory point, as viewed in solar wind coor- 
dinates, on solar vind direction. For the partic- 
ular angles indicated, the shift in X coordinate 
of a trajectory point can be as great as a quarter 
of the Venusian planetary radiuj and obviously 
results in substantial differences in predicted 
flow and magnetic field properties. Although 
Spreiter and Rizri 11972] noted the importance of 
such an effect in the interpretation of the date 
from Mariner 4 in its flyby of Mars, the effects 
of angular shifts in the solar wind direction are 
usually ignored in most discussions of the inter- 
action. The present results, however, indicate 
that this purely geometrical effect can be very 
significant, even for directional shifts of less 
than 5', and must be accounted for in any 
quantltatlve comparison of theoretical and 
observational results. 

The selection of values for M_ and Y in 
Figure 11 represents an attempt to resolve the 
uncertainty in the appropriate values for these 
parameters. Because only solar vind ion temper- 
atures from the OPA were available for orbit 6, 
the initial calculations of M^ were based on the 
assumptions that Tg • Tj, which leads to M^ • 13.3 
and Y ■ 5/3, as appropriate for a gas with parti- 
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rig. 17. Conparison of observed (OMAG) and theoretical time histories of the magnetic 
field for Pioneer Venus orbit 3 based on inbound solar vind interplanetary eoncicions 
using a gasdynamic solution for • 3.0, Tf » 5/3. (a) Magnetic field magnitude. 

(b) Magnetic field components. 


cles possessing 3 degrees of freedom for their 
molecular notion. The resulting location for the 
shock wave is not in good agreement with the 
observations. To investigate whether changing 
the value for T to 2, which corresponds to a gas 
with 2 degrees of freedom as might be supposed to 
occur if the magnetic field sufficiently aligns 
the molecular motion, the calculations were 
repeated with M^ “ 13.3 and T " 2. The results, 
indicated by the dashed line, are in improved, 
but still not completely satisfactory, agreement 
with the observed shock locations. Finally, if 
it is assumed that the oncoming interplanetary 
electron temperature is not equal to the ion tem- 
perature but is substantially higher, we are led 
to consider low values of the order of 3-5 for M . 
Results for M^ » 3.0, Y • 5/3 are Indicated in 
Figure 11 by the solid line. It may be seen that 
results calculated for these values and for the 
Inbound solar wind direction are in very good 
agreement with the theoretical results. These 
comparisons make evident both the possibilities 
and the uncertainties in the use of the theory at 
the present level of understanding of the parameter 
selection process. 

Figure 12 displays the time history comparisons 
of the predicted bulk plasma density, speed, and 
temperature with OPA data. The solid lines with 
circles correspond to results based on inbound 
interplanetary conditions, while the dashed lines 
with squares correspond to outboxind conditions, 
'/hile the one data point in the ionosheath for 
each quantity is in general agreement with the 
theoretical calculation, Che lack of more detailed 
plasma data prevents a definitive conclusion. The 
OFa requires approximately 9 min to acquire suffi- 


cient data CO enable determination of the bulk 
plasma quantities. This time interval usually 
presents no problem when the spacecraft is in the 
interplanetary solar wind, but the large resolu- 
tion time effectively averages the plasma quan- 
tities in the ionosheath over such a large spatial 
range that only overall comparisons of the bulk 
plasma properties are possible with this 
instrument . 

The situation is quite different for the 
magnetic field, since the OMAG provides essentially 
Instantaneous measurements. Figure 13 displays a 
comparison of the data from orbit 6 with two sets 
of theoretical results based on the Inbound and 
outbound Interplanetary conditions. In these com- 
parisons both sets of predictions exhibit good 
agreement with Che data for the appropriate part 
of the trajectory. Thus, on the Inbound leg the 
predictions based on the inbound interplanetary 
conditions are in good agreement with the data, 
while the outbound-condition predictions are 
notably inferior, particularly with regard to 
shock crossing; and conversely. Since steadiness 
of the interplanetary conditions was a criterion 
used in the selection of orbit 6 for these com- 
parisons, the results of Figure 13 serve to 
illustrate the need for simultaneous measurements 
of interplanetary conditions anc plan’-tary flow 
field conditions in more definitive evaluations 
of the theory. 

Results fox corresponding comparisons for 
orbit 3 are given in Figures 16, 15, and 16. In 
Figure 16, bow shock locations are presented for 
five different combinations of M„, and 7. The Mach 
numbers M^ - 7.38 and 5.96 correspond to the in- 
bound and outbound interplanetary conditions. 
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Fig. 17. (continued) 


respectively, for K',^1, and as oeasured for 
the ions by the OPA; the tvo values for T ■ 5/3, 

2 are indicative of present uncertainty about the 
effective ratio of specific heats for the solar 
wind plasma. The bov shock location for “ 3 
and T « 5/3 given previously in Figure 11 for 
orbit 6 is also included, because once again it is 
closer to the observed shock crossings. Also 
provided in Figure 14 are the orbital trajectories 
as viewed in solar wind coordinates for the in- 
bound (n, tp) “ (3.3°, 0.15‘) and outbound (fi, dp) 
• (3.7‘, 4.9”) solar wind directions. 

Tne comparisons for the bulk plasma properties 
for orbit 3 shown in Figure 15 are marked again by 
the paucity of observational data. What informa- 
tion there is suggests agreement for the plasma 
speed and density but a notable discrepancy for 
the temperature Tnis is thought to be indicative 
of the manner in which the bulk properties from 
the theoretical model are interpreted in relation 
to the observational measurements; l.e., the theo- 
retical values correspond to those for a single- 
component plasma in which it is considered that 
the electrons and ions have equal temperatures, 
whereas the measurements are for the ions only in 
e multicomponent plasma in w*hich the ion and 


electron temperatures are usually quite different. 
A more penetrating analysis of the meaning of the 
temperature in the theory and its relation to that 
deduced by the experimenters from their date 
appears to be a necessary and important subject 
for future study. — 

Results for the magnetic field comparisons are 
displayed in Figure 16, which provides time 
histories of the magnitude and individual com- 
ponents based on both inbound and outbound inter- 
planetary conditions. Although the shock crossing 
are somewhat in disagreement, since the calculated 
results are for » 7.56, 6.96 and Y » 2, the 
general trends displayed by the calculated and 
observed results are notably similar, particularly 
if account is taken of the drift v’ith time along 
the trajectory of the transition from the calcu- 
lations based on inbound conditions to ones based 
on the outbound conditions. 

To demonstrate the improvement obtained in 
magnetic-field results when a gasdynamic flow 
field solution is employed which more closely 
agrees with the observational bow shock location, 
the magnetic field time histories from orbit 3 are 
compared in Figure 17 using gasdynamic results for 
“ 3.0, Y « 5/3, and the directions Cp) • 
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Fig. 18. Illustration of effect on bow shock location of variability of oncoming 
solar wind direction; ■ 8.0, Tf ■ 5/3 flow past a sphere/cylinder ionopause 
obstacle. 


(3.3°, 0.15°) measured at the Inbound shock cross- 
ing. As can be seen, there is a marked improvement 
in the agreement near the bow shock and quite good 
agreement throughout the ionosheath for both the 
magnitude and the individual magnetic field com- 
ponents. In addition to illustrating again that 
further study needs to be made of the parameter 
selection process required to relate Che calculated 
and observational results optimally, these compari- 
sons emphasize once again the need to account in 
the calculations for the actual direction of the 
interplanetary solar wind flow. 

To illustrate the significance of this point. 
Figure 18 has been prepared to show chat even a 
modest change of 16 “ in the direction of the 
incident solar wind can have substantial effects 
on the location of the ionopause and bow wave, as 
viewed in a coordinate system aligned with the 
Sun-Venus direction which is typical of most 
presentations of observational data. It would be 
of considerable interest to determine how much of 
the scatter of data for the location of the bow 
wave may be attributable to this simple factor. 

A further point of growing concern is associ- 
ated wich effects of significant changes down- 
stream of the terminator of the ionopause shape 
from that calculated herein. More detailed com- 
parisons will have to be reserved for future 
studies, but the results of Figure 19 for the 
t.neoretical characteristic lines emanating from 
tne intersection of the ionopause and the termi- 
nator define the region of influence of a small 
inward or outward taper of the ionopause profile 
downstream of the terminator. For the entire 
region displayed extending downstream to more than 
16 \enus radii, these results indicate that there 
IS nc accompanying change in the location of the 
bow wave. It may also be noted that the point 
designated by number 2 in the Mariner 5 data from 
Venus (see Spreiter and Alksne 11970] for an 
earlier comparison with theory), where the lono- 
sneath plasma changes from comparatively quiet 


to disturbed, is remarkably close to the charac- 
teristic line from the terminator. Such behavior 
could result from unsteady fluctuations of the 
ionopause surface beginning approximately at the 
terminator. The data from Pioneer Venus should be 
examined to determine whether this is a general 
property of the plasma in the vicinity of these 
characteristic lines or perhaps just the result 
of some passing transient occurring at the time 
of the Mariner 5 encounter . 


Concluding Remarks 

The present application of advanced computa- 
tional procedures to the modeling of solar wind 
flow past Venus was undertaken to Improve the 
accuracy and utility of the theoretical predic- 
tions. Starting with the steady, dissipationluss, 
magnetohydrodynamic model for axisymmetric, super- 
sonic, _ super-Alfvdnic solar wind past a nonmagnetic 
planet with a shielding ionosphere, a number of 
Important theoretical extensions have been devel- 
oped and implemented. These Include the capability 
for treating lower Interplanetary gasdynamic Mach 
numbers down to about 2 and a wider variety of 
ionopause shapes Including a new family of shapes 
which includes the effect of gravitational varia- 
tion in the scale height. Additionally, the 
capability for determining the plasma gasdynamic 
and magnetic field properties along any arbitrary 
trajectory, accounting for an arbitrary oncoming 
direction of th . solar wind, has been established. 
All of these developments have been incorporated 
into a current computatlon.il model to enable 
detailed calculations of the solar wind inter- 
action with Venus. The model has been compre- 
hensively described by Scahara et al. (1980]. 

Also Included in that account is an extensive 
catalog of results for a wide range of flow con- 
ditions and Ionopause shapes representative of 
Chose that might be anticipated at Venus or at 
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sone other nonaagnetic planet vith a shielding 
tonosphere . 

Finally, theoretical results determined from 
the model have been compared vith data from two 
cf the early orbits of the Pioneer Venus orbiter. 
Tnese comparisons have indicated the importance, 
heretofore largely neglected, of the directional 
variability of the oncoming solar wind. These 
results, taken in toto, serve to verify the basic 
theoretical model and the effectiveness of the 
nitserical solution, although it is evident that 
further study should be applied to certain 
questions relating to interpretation of plasma 
quantities like the temperature. The results 
denonstrate also the value of the present theo- 
retical procedures as a research tool capable of 
routinely provlding--at small computational cost 
an: tr. a form directly compatible vith observa- 
tions — oetails of the solar wind/planetary 
atmosphere interaction process not previously 
attainable. 

With regard to future uses as well as iaprove- 
ner.ts of the present model, the obvious need for 
a detailed study involving comparison between 
theory and observations for a large number of 
orbits cf the Pioneer Venus orbiter is clear. 

Or. the basis of the preliminary comparisons for 
orbits 3 and 6 the calculated magnetic field 
appears to be remarkably accurate for relatively 


quiettime conditions. Similar comparisons of the 
plasma properties indicate a need for an improved 
interpretation of the results from the single-fluid 
theory in terms of multicomponent measurements. 
Questions regarding the possible suppression by 
the interplanetary magnetic field of' the number 
of degrees of freedom of the plasma require 
further study and could be clarified through 
systematic comparisons with data. Additionally, 
observations from the Pioneer Venus orbiter of 
the nlghtside ionosphere of Venus have revealed 
a more complex and dynamic structure than sus- 
pected. These observations point, in particular, 
toward the need for Improvement of the simple 
model used in the present method for' the deter- 
mination of the ionosphere boundary downstream of 
the terminator. This improved determination would 
involve an iterative procedure in which a balance 
of the sum of the solar wind gasdynamic and 
magnetic pressures along the ionopause surface 
would be maintained against the ionospheric p’ as- 
sure, not necessarily assumed to be spherically 
symmetric. The present method, which balances 
the extended Newtonian gas plus magnetic pressure 
distribution against the ionospheric pressure, 
might represent the first step of such an itera- 
tion; but it is also possible that dynamic effects 
associated vith intermittently passing volumes cf 
magnetized plasma tom from more upstream regions 
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of Che lonopause might greatly impair the general 
utility of any stationary model for Che ionopause 
dovnstreaB of the terminator. This is another 
important topic awaiting further study. 
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A gas dynamic model for solar wind flow around Venus has been compared with Pioneer Venus 
orbiier plasma analyzer measuremenis from limes when the solar wind flow seemed sieadiesi. The 
comparisons were made near the lerminaior. When the observed and model bow shock locations are 
matched, the model agrees fairly well with the observed parameiers: the components of the" flow 
velocitx- and macneiic field, and the proton number densiiv and isotropic lemperaiure However, the 
Mach numbers required lo fit the observed bow shock locations are less than 2'3 those estimated from 
the measured parameters of the free-stream flow, because the bow shock locations (near the 
terminator) are farther from Venus than expected. Sometimes the measured flow speeds appear lo be 
retarded near the ionopause. 


1 Introduction 

The Pioneer Venus orbiier spacecraft [Colin, 1980], 
launched on May 20, 1978, carries plasma analyzer [Intrili- 
gator et al., 1980] and magnetometer [Russell et oL, 1980] 
experiments that produce data that permit study of the flow 
of solar wind plasma around Venus. The plasma analyzer 
experiment uses 90‘ electrostatic deflection followed by 
current measurement throughout each spacecraft rotation, in 
five different angular sectors, using separate electrometer 
tubes. This study compares ionosheath data from five orbits 
with calculated values from the gas-dynamic model of 
Spreiter and Stahara [I980fl, b]. Because of the highly 
inclined orbit, with periapsis initially well w'iihin the iono- 
sphere. the data are generally acquired in the polar regions. 

2. Flow-Field Data From Venus' Ionosheath 

Solar wind plasma proton parameters are obtained from 
the Pioneer Venus orbiter plasma analyzer data by least- 
squares fils of an isotropic temperature, convecting Maxwel- 
lian velocny distribution, to flight currents. When the space- 
craft's periapsis is on the nightside of Venus, relatively 
distinct inbound and outbound crossings of Venus' bow 
shock arc expected. For orbits w'lth periapsis on the dayside, 
on the other hand, extended time intervals within the lono- 
shcath and wjike are expected, so that it becomes more 
difficult to judge the steadiness of the free-stream plasma 
flow Consequenth’. data from nightside periapsis orbits 
only, through the first four 'eclipse seasons' to orbit 803, 
were examined for relatively constant solar wind proton 
vciocity (including the flow direction), number density, and 
temperature. For the relatively constant cases, the amount 
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of ionosheath data available, and the steadiness of the 
interplanetary magnetic field also had to be checked. Final- 
ly, five orbits were selected as those with the steadiest 
inierplaneury conditions, most suitable for detailed compar- 
ison with the gas dynamic model. Table 1 gives interplane- 
tary conditions for the orbits used here. 

In Table 1, u is the measured proton bulk speed (space- 
craft reference frame), n is the proton number density. is 
the isotropic proton temperature, B is the magnetic field 
magnitude. M, = v/[yR{Tr -r Tp)/m]' ' is used for the sonic 
Mach number, with the ratio of specific heats, y - 5/3. 
except as noted, and the mean molecular weight m = 1: 

«= u/(B/(jv7m^)’'^]; Mfj_ is the perpendicular fast magneto- 
sonic Mach number [Clemmow and Dougherty, 1969]: and 
Bp = nkTpl{B'!2no) is the ratio of proton gas pressure to 
magnetic pressure. In addition, k is Boltzmann's constant. R 
is the universal gas constant, fio is the vacuum permeability, 
and nip is the proton mass. The electron temperature. T,. is 
not directly measured during these periods, but has been 
estimated by using the measured proton and approximated 
helium number densities, together with an empincal expres- 
sion of Sutler and Sctidder [1980]. logio TrlK) - (0 185 r 
0.016)logio /ir<cm" ') T (4.74 z 0 003). Here T, is a tempera- 
ture for the 'core' of the electron velocity distribution, and n, is 
the electron number density We use 7^ ~ 7r and n, — n - 2 
n„, and no is a helium number density estimate. Estimates of 
njn vary from 0.02 (orbit 65) to 0.07 (orbit 582). as a result, 
the effect of solar wind helium on the values of M, and is 
less than -i- 13%. and thus is negligible. 

In Table 1. the column labeled M^ gives the Mach 
numbers for which the gas dynamic model shock crossing 
locations agreed w'ith those observed, as discussed below. 
The column of 'comments' in Table 1 gives a rough desenp- 
livc indication of the solar wind stream regime at the times 
used here. The solar wind flow is apparently much steadier 
for the ionosheath traversals discussed here, than for orbits 3 
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TABLE 1 . Free-Stream Solar Wind Condiiions 


Orbit 

km/s 

cm*’ 

Tp.’K' 

B, 

nT 


M. 

(est.) 

M/j, 

M. 

Pr 

Comments 

63 inbound 

360 

13.4 

63.000 

8.8 

6.9 

7.8 

5.2 

2.8 

0 38) 

Lou -amplitude 

63 outbound 

362 

J4.3 

53J200 

8.8 

7.1 

8.1 

5,4 

3.4 

0.34} 

speed maximum 

66’ 

280 

43.01 

80.0001 

10.8 

fe7.8 

£5.4 

-4.5 

2.65 

fel.O 

Trough between 











mcideraie speed 











streams 

68 inbound 

377 

21.5 

157,000 

11.2 

7.2 

6.3 

4.7 

2.55 

0.931 

Low maximum in 

68 outbound 

366 

21.8 

181,000 

10.3 

7.6 

5.9 

4.7 

2.55 

1.3 I 

speed 

j'so* T 5/3 

334 

29.0 

53,500 

10.6 

7.8 

7.2 

5.3 

2.65 

0.48 

Trough between 

r = 2 






6.5 

5.0 

3.2 


higher speeds 

SS2 inbound 

382 

-13.0 

-125.000 

8.2 

7.7 

7.0 

5.2 

3 1 

0.84] 

Trough betueen 

582 outbound 

389 

-13.0 

-100.000 

7.4 

8.7 

7.6 

5.7 

2.4 

0.82} 

high-speed 


streams 


‘Steadiest solar wind, comparison with model should be most reliable. 
tNonthermal portion of proton-velocity distribution not included in these values 


and 6. for which Spreiier and Siaharo [1980i>] compared 
flight data with the same gas dynamtc model. 

The locations of the plasma dau samples used in the 
present comparisons, in an axisymmetnc coordinate system 
with the free-stream solar wind flow direction as the axis, are 
shown tn Figure 1 . These locations may be characterized as 
lying generally in the terminator region of the ionosheath. 

3. Description of Gas Dynamic Model 

The lime dependent dissipationless magnetohydrodynam- 
ic equations for conservation of mass, momentum, energy, 
and magnetic field, ignoring heat conductivity, viscosity, 
and electrical resistivity [Spreuer et oL, 1966; Spreiter and 
Alksne. \91d-. Spreiter and Stahara, 19806], excluding gravi- 
tational terms, are used here for calculation of model iono- 
sheath flow fields for comparison with the spacecraft dau. 
These equations are 


ap 

bt 


+ — (piit) = 0 
8-rt 


( 1 ) 


b , 5 

— (pv,) -i- 

bt OXi 


pL'iVl + p6,|. - 


ha 


2mo 


(2) 


Spreiier and Alksne, 1970], and the magnetic field is subse- 
quently determined by using the velocity and density as 
known quantities. Because dW'^) are dissipation-free, dis- 
continuities in p. V . p. T. and B are permitted in the solutions 
The propemes of the discontinuities are implicit in {1H4) 
because they are expressed in conservation form, and corre- 
spond to the generalized magnetohydrodynamic Rankine- 
Hugoniot shock-jump equations, if no approximations are 
introduced. Similarly. (l)-<4) correspond to those of pas 
dynamics and the limiting case of high .^Ifvdn Mach number 
when the magnetic terms are omitted from (2) and (3 1 
[Spreiter et al., Spreiter and Alksne, \910\ Spreiier and 
Rizzi, \914\ Spreiier and Siohara, 19806] With this approxi- 
mation, (1H4) are then solved for the axisymmetric flow and 
magnetic field and the shock location. The boundary condi- 
tions applied are (1) the gas dynamic shock jump equations 
al the shock, (2) supersonic flow along a dow-nstream 
surface, onented normal to the axis of symmetry', and (3) 
flow parallel to both the lonopause surface, and to the 
sUgnation streamline [Spreiter and Siahora, 1980o]. 

The axisymmetric obstacle shape (lonopause) is obuined 
by integrating an ordinary' diflTerential equation obtained 
from the Newtonian approximation for pressure balance 
across the boundary. 


b 

fptr 


, b 
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■-(R, - R.y 

bt 


2ha) 


put 


T e + — 4- 

W J 

cos-i// = p,„ exp 
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(5) 


bB, b bB, 

— = — - uiB,) -r- = 0 (4) 

bt OX I OX, 

where the Po>-ntinp vector is Si = n/pol(i>,6‘ - and 

the ideal gas equation of state, p = pRTIm is used Here p. p. 
t,',, T. and e refer to density, pressure, velocity component, 
temperature, and internal energy per unit mass, respcctivc- 
Iv. The value m - 1/2 is used for fully ionized hydrogen, and 
B IS the magnetic field vector. 

In the momentum and energy equations above, the as- 
sumption that A/o is high is invoked, so that the magnetic 
terms in (21 and (31 are disregarded as small. The mechanical 
ponion>- of the mapneiohydrodynamic equations thereby 
reduce to those of pas dynamics [Spreiter ci al.. 1966, 


by uking p„ = p,„ [Spreiter et al., 1970). Here p„ is the 
SUgnation pressure of the solar wind at the nose of the 
ionopause. and i//is the angle between the outward normal to 
the ionopause boundary and the free-stream solar wind flov 
direction. 

The ionopause altitude at the nose of the ionosphere is 
specified, and an isothermal upper atmosphere with gravita- 
tionally varying scale height, h = RTlmg. is assumed 
[Spreiter and Siahora, 19806] Here, g = g.lRJr)' is the 
local acceleration of gravity . a subsenpt s denotes a value at 
Venus' surface. R, is the radial location of the ionopause 
surface, and r is the radial distance from the center of the 
planet. T is the lonosphenc temperature, and R,„ is a 
constant reference radius, taken here to be R, - H. with H 
the ionopause altitude at the ionopause nose Both T and m 
arc assumed to be constant in the ionospheric region of 
interest Effects of electric currents ssahin and at the bound- 
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Fif I Locations of samples of lonosheath plasma data used in this siud>’. ploiied in a cylindrical coordinaie s\-siem 
ui;h the avial coordinaie aliened wuh ihe free-sircam solar wind flow direction The ho» shock shape matches shock 
crossing locations for the outbound leg of orbit 5£!! (Af. = 2.0. wuh the lonopause determined from (.Sl 


an' of the ionosphere are assumed lo be negligible in 
deicrminmc the pressure profile of ihe ionosphere. 

Spreuer and Staharu [1980c/] have described Ihe compuia- 
iional methods used to calculate the flow field. In the nose 
region extending down.strcam lo ihe vicinity of the icrmina- 
to.'. both subsonic and supersonic flow occurs. Here, the 
fluid properties arc calculated by using the second-order 
accurate, noniterative finite difference algonthm of Beam 
and Warminy [19261. with Euler time differencing, to solve 
lime dependent. Eulcrian gas d\ namic equations for axisym- 
meinc flow, derned from equations (1) to (.2) with the 
magnetic terms deleted The computational mesh points are 


located ever)’ 4' from the stagnation point to the terminator, 
and at 17 intervals between the lonopause and the bow 
shock. Further description of the wa\' that the boundary 
conditions are imposed is given b>' Sprcnct and Siuharn 

[\mb] 

Dow'nstrcam from the region m which all subsonic flow is 
found, the time independent form of the same axisv'mmeiric. 
gas dynamic flow equations is soKed by using the shock- 
capturing. marching method of Kuilcr ei at. [ 197?]. At least 
15 steps are used proceeding downstream from the termina- 
tor to 1 /?! . with mesh points at 17 intervals, again, between 
the lonopause and the bow shock. 
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The comparison of calculated and observed values was 
done for flow fields with Mach numbers for which the model 
shock locations agreed with those observed, as well as for 
Mach numbers estimated from upstream measurements. In 
the case of multiple bow shock crossings, the innermost ones 
were used to compare locations with the model. The Mach 
numbers required for agreement with the observed shock 
locations are given in Table 1 in the column headed 
These flow-field calculations use hJR,„ = 0.03, as did 
Spreiter and Siahara [1980h], except for one case described 
below. This value corresponds to a scale height of the 
da\'side ionosphere of about 200 km. in agreement with the 
measured ion and electron temperatures of about 2000'K and 
3000°K. respectively [Miller et al., 1980]. An alternative 
method for calculating the ionopause shapes has been pre- 
sented by Knudsen et al. [1982]. The flow field using this 
shape was also calculated for one leg (orbit 582 outbound) 
and the results compared with those obtained by using the 
earlier ionopause shape 

Small variations of the ionopause nose distance, R,„, were 
introduced for the flow-field calculations for the various 
orbits, by using for ionopause nose altitude. R,„ - R,. the 
ionopause alvi’.udo obtained from the empincal expression of 
Braie et al. [1980]. The measured free-stream proton con- 
vective pressure was substituted in the empirical expression 
for the magnetic field pressure measured at the outer bound- 
ary of the ionosphere. The nose altitudes used ranged from 
296 km for orbit 66 (an upper limit) to 343 km for orbit 63 
(inbound) 

The streamlines for the case of orbit 582 outbound, with 
the ionopause shape as described by (5). are given tn Figure 
2 _ 

4, Co.mparison of Observations With Calculated 
Results 

Of the five steads orbits (Table 1). comparison with model 
calculations is first shown for orbit 129 inbound (particularly 
steady) and orbit 582 outbound (larger number of samples 
dunng the lonosheath traversal), in Table 2. In this table, A" 
and R- refer to the locations of the observations in a 
cylindrical coordinate system, in units of Venus radii (1 Ry 
= 6051 km). The axial coordinate X is along the sun-Venus 
direction and positive toward the sun. (Figure 1 . on the other 
hand, is a plot in cylindrical coordinates (Xw; R,*), with the 
upstream solar wind fiov.' direction as an A",* axis, of the 
locations of the comparisons.) The gas dynamic calculations 
used the sonic Mach numbers A/, estimated from free-stream 
measurements (Table II Proton convective pressures. 
nm^v'. of 5.4 >• lO'^ and 3.3 x 10*" dynes cm"' for the free 
stream are used to deduce ~?00 and —330 km ionopause 
altitudes at the nose for orbits 129 and 582. respectively. One 
immediaiely finds, as did Spn-iirr and Siahara [1980i] for 
iv,o other orbii^ that the calculaied bov. shock locations lie 
well wiihin ihe observed locations. Also, the calculaied 
speeds and proton number densities tend to be larger than 
the ooserved \alucs For orbit 582 (outbound), the calculat- 
ed macriciic field magnitudes tend to increase as the bow 
shock 1 ^ approached, whereas the obsersed values have an 
oppo'iie trend. 

Funner comparisons arc now given for flow fields with 
loucr M.ich numbers for which the model shock locations 
agree \'>iih those obsersed. Sprciicr and Siahara (1980/)] 
.iieueo th.it suer, losser implied .Mach numbers might arise 


from elevated electron temperatures in the free stream. 
Table 3 gives a dctaDed comparison of observed and calcu- 
lated values for orbit 129, one of the cases with the steadiest 
solar wind. 

Least-squares fitting uncertainties are given in separate 
columns of Table 3 for the proton bulk speeds (oy.,.) and the 
polar and azimuthal flow direction angles (o/poi and 0}.^,:); 
for the proton number densities and isotropic temperatures 
the corresponding uncertainties are given in the same col- 
umns. These uncertainties, which are the changes in value of 
one plasma parameter that results in an increase of x" t>y 
one, are taken as one standard deviation [see Bevinfuan, 
1969] and should represent only part of the total uncertain- 
ties. Response function uncertainties are one source of error 
not accounted for by the uncertainties in the plasma parame- 
ters given in Tables 2. 3. and 4. The 1. 2. and 3 components 
of the bulk velocity (magnitude v) are positive in the antiso- 
lar direction, the direction of planetan' motion, and the 
northward direction norma! to Venus' orbital plane, respec- 
tively; these components make up a nght-handed. orthogo- 
nal system. The velocity components have been corrected 
for tips of the spacecraft spin axis av\’ay from an ecliptic 
north-south direction, for the inclination of \ enus' orbit to 
the ecliptic plane and for the spacecraft velocity relative to 
Venus. 

The observed and calculated temperatures in the table are 
not directly comparable. The observed temperatures are 
obtained from flight data under the assumption of isotropic 
proton temperatures. The theoreiicalK calculated tempera- 
tures are obtained from the one-fiuid gas dvnamic model 
described above: they are used for the sum of the proton and 
electron temperatures. 

The resolution of the basic magnetic field measurements is 
il/16 nT [Riis.’tell et al.. 1980] For orbit 129. a hich- 
resolution (—100 samples per min) plot was supplied, and the 
average magnetic field component values were obtained b>' 
eye from the plot. For the other orbits, the magnetic-field 
data were made available as tabulations of averages over 64 
s. The three components of the magnetic field make up a 
right-handed, onhogonal system, with component 1 positive 
in the solar direction, and component 3 positive northward in 
the direction normal to Venus' orbital plane. 

For orbit 129, five lonosheath plasma samples are avail- 
able: the upstream plasma parameters appear similar to 
those for the other orbits considered for this study, but the 
proton temperature is at the low end of the range. To a 
greater extent than for most of the other cases, use here of 
expected values for electron temperatures results in a calcu- 
lated sonic Mach number too large to permit fitting the 
observed shock crossings, the calculated shock locations 
would fall inside of those observed. By matching the gas 
dynamic shock location to the observed shock locations, it is 
found that the calculaied shock for a .Mach number of 2 65 
(and a 297-km ionopause nose altitude) can both match the 
last inbound bow-shock crossing (of three crossings) ob- 
served in the magnciomeier data, and falljusi a fev. seconds 
earlier than the first of three prominent bow -shock crossings 
ohsersed in the plasmu anahzer data, on the outbound leg 
The second of these three outbound bov,’-shock crossings is 
only -70 s later than the first. (It is difficult to obtain an 
inbound shock crossing time from the plasma anaU zer data 
because a pronounced deflection of the polar floss angle is 
not obsers ed ) 
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Except for the final speed sample, near ihe outbound bow 
shock, the orbit 129 lonosheath speeds remain lower than the 
calculated \alues In particular, the speed sample closest to 
the tonstpause on the inbound lep is only about half the 
Calculated \ dlue The calculated proton number densities are 
all within the fining uncertainties of the measured values, 
given m Table The observed proton fluxes for this orbit 
are no; given in the table, however, they generallv agree with 
the calcui..ted values, within fitting uncenainiy. 

Tne loni'sheath plasma flow deflections, obsened pener- 
allv in the polar regions, are typically less than the calculated 
values m the nonhward-soulhward direction (component ?) 
and greater than the calculated values in a plane perpendicu- 
lar to the northward-southward direction (except for the 
sample cK'vesi to the lonopausc on the inbound leg). This 
ivhxe'rv .ition seemx to suggest an asymmetric obstacle. 


The calculated magnetic-field magnitudes for this orbit are 
now consistently less than the measured values (for the 
inbound leg sample closest to the lonopause. the two values 
are quite close, how-ever). Except close to the ionopause. a 
similar tendency is present for all of the other orbits. This is 
similar to the earlier results of Spinier and Swhara | IVKOh] 
(orbits 5 and 6) As an appro-ximaie descnpiion for orbit 129. 
the calculated northward-southward magnetic field compo- 
nent tends to agree with the observations; the calculated 
magnetic-field component along the sunward direction tends 
to be less than the observed values; and the remaining 
orthogonal magnetic-field component is observed to be both 
greater than (outbound lep) and less than (inbound leg) the 
calculated values. 

If y. the ratio of specific he. its. is assumed to be 2 rather 
than a cas-dvnamic .Mach number of .'.2 is found tv' fit 
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TABLE 2. Comparison of Observed and Calculated Values (M, Used in Model) 


X, Rv 

Rz, R\' 

V . km/s 

Ucaie. kffl/S 

n. cm"’ 


B . nT 

5«ic- nT 

-0.507 

2.327 

238 

Orbit 129 Inbound, M , = 
* 63 r 10 

7,2 

• 

22.3 

• 

-0.-147 

1.956 

238 

280 

42= 8 


84 

21.0 

28.1 

-0.365 

1.559 

-136 

271 

144 = 91 


46 

18.7 

21.5 

0.342 

1.250 

272 

Orbit 582 Outbound, M , = 
285 12 = 5 

7.6 

23 

26.0 

23.1 

0.385 

1.372 

271 

286 

13 = 3 


30 

25.5 

22.3 

0.418 

1.493 

266 

292 

15 = 6 


36 

25.2-24.2 

23.5 

0.447 

1.631 

254 

303 

26= 7 


46 

24.8-24.1 

26.1 

0.470 

1.764 

251 

9 

40 = 15 


• 

20 4 

• 

0.490 

1.908 

304 

9 

40 = 6 


• 

17.8-16 4 

■ 


"Outside of calculated bow shock. 


the obsen-ed shock locations. Spreiter and Siahara [1980fc] 
noted that such an assumption caused model shock locations 
near the terminator to move outward and match better their 
observed locations. The results for calculated parameters for 
this case are also given in Table 3; the changes from the -y = 
5o cases are negligible, except that the calculated tempera- 
tures increased by about 30-509c. 

The obsen'ed and calculated bulk speeds, number densi- 
ties. and magnetic field strengths for the other orbits of this 
study are compared in Table 4 The speeds are presented in 
the reference frame moving w'lth Venus. Since the magnetic 
field data v. ere supplied as tabulations of averages over 64 s, 
there are cases in which the times of these averages do not 
coincide well with the times specified for the plasma data 
For those cases, values for these averages on both sides of 
the plasma data times are given. (The plasma dau usually 
require at least 90 s for acquisition.) The values in parenthe- 


ses for orbit 582 (outbound) in the table result from a 
calculation using the lonopause shape of Knudsen e) cl. 
[1982]: they are discussed later 
In Tables 2, 3. and 4. the calculated speeds are generally 
larger than the observed values. The calculated number 
densities and magnetic field magnitudes tend to be smaller 
than the observed values The comparisons are discussed in 
the next section. Also, when the model shock locations 
match those observed, there is a slight tendency for the 
measured proton fiux (nv) to be larger than the calculated 
value, excluding orbit 582; namely, of 14 flux values from 
these earliest four orbits, four measured values are larger 
than the calculated values, including the fitting uncertainties 
of the density, and only one measured value is similarly 
smaller than the calculated value For orbit 582. onh the 
outbound portion has w-ell-defined density measurements, 
and the observed flux values nse by a factor of about 3. 






TABLE 3. 

Comparison of Measured and Model Values for Orbit 129 






V. 

<Tr,. 


Of.az- 

0/.roi. 

Ul. 

^ » .Cfcic • 


C..IC. 

lx. 


.V. /?, 


km^s 

km's 

km/s 

deg 

deg 

km.'s 

km.'s 

km's 

km's 

km/s 

km's 







Inbound 






-0..507 

2 327 

238 

5 

275 

(278) 

0.8 

1.2 

234 

271 

(274) 

30 

15 

(15) 

13 

42 

(42) 

-0 447 

1 936 

238 

6 

278 

(280) 

0.7 

1.9 

234 

275 

(276) 

28 

0 

(-1) 

10 

AA 

-0 765 

1.5.S9 

-136 

26 

286 

(289) 

5 7 

4.6 

135 

28.3 

(285) 

-0 

-24 

(-241 


38 

(4(1) 







Outbound 










272 

Cll) 




'JiJ 

-102 

- 7 s 

-80 

-Of 

0 .301 

1 441 

236 

6 

0.8 

1.1 

199 

(247) 

(-76) 

(-97) 





263 

(266) 






-49 

-70 

- 1 1 1 

0 42i. 

1 8.3.‘ 

259 

8 

0.9 

2.1 

243 

12.37) 

-62 

(-.30) 

I-ID) 

n 

n, 




B. 



^ 1 u-ii • 




B.,. 

cm ‘ 

cm‘ * 

Tn 

'K 

Feu. 'K 

nT 


B, nT 

nT 

nt 

nT 

Bi nT 

nT 







inbound 







S6 


.000 

276 000 

22.3 

19 9 

11 4 

7 0 


-9 1 


- 16 3 

6.3 = 10 

155) 


UfW.OOO) 

(19.8) 

|7 0) 


(-9 21 


1- 16 1) 

42 ; S 

(-‘1 

£208.000 

270.(Ki0 

(.’94.000) 

21 0 

18 4 
118 1) 

13 8 

8 1 
II- Oi 

-7.2 

-9 3 
(-9 2i 

-14 1 

- 1? 6 
( - 1 ' - 1 

144 - Vi 

1 'h) 

S .393 

.(KK) 

2.34.000 

(3.39.000) 

18.7 

17.0 
(16 5) 

16 5 

11 3 

(10 81 

-2 5 

-8 5 

1-8 41 

-8 3 

-9 4 

(-9.2) 







Ouihiiiind 






r s 

« 

1411 

£259 (KK) 

.3)9.000 

(413.000) 

30 7 

21.6 
(21 11 

— 7 6 

A 1 

1-4 4) 

- 1.8 8 

-116 

(-1141 

-23 1 

- r 7 
(-17.2) 





297.000 


17 9 


-.3 1 


-8 8 

1“ ^ 

-15.2 

hh ; 1-1 

lS4i 

<-2J9.(KKI 

(4.S(i.(KXn 

23 .2 

(17.5) 

mmm 

1-3 0i 


(-8 71 


1- 14 9| 


r, :,i.'l..itfc’ \ jiuc' tr, = 2 v = S'J. and m paremhoes 4/, = .’ 2 ■> = 2 
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TABLE Comparison of Observed and Calculated Values 


Orbit 

X.Ry 

Rz- 

u, km/s 

fulc< 

km/s 

n. cm"’ 

cm"^ 

B. nT 

Scale- nT 

63 inbound 

0 4525 

1.719 

£260 

271 

£61 = 23 

29 

dau gap 

15.2 

65 outbound 

-0 4645 

2.035 

311 

313 

30 = 8 

22 

12.9-16.7 

16.0 

63 outbound 

-0.206 

2.397 

314 

305 

76-11 

32 

22.0 

18.5 

66 outbound 

-0.635 

1.711 

236 

256 

52 = 6 

44 

13.1 

13.5 

66 outbound 

-0.300 

'2.207 

220 

236 

97 = 11 

72 

23.2-23.5 

18.4 

68 inbound 

0 477 

1.766 

328 

282 

-74 = 13 

45 

£25 4 

18.2 

68 inbound 

0.180 

1.618 

284 

298 

75 = 15 

36 

-22.7 

16.8 

68 inbound 

-0.116 

1.431 

272 

328 

26 = 14 

25 

14.7-13.1 

15.3 

68 outbound 

-0.318 

2.153 

302 

315 

56 r 8 

33 

19.1 

15.1 

582 inbound 

-0.744 

2.659 

<418 

326 

>7.2 

27 

11.3-13.2 

11.5 

582 inbound 

-0.673 

2 405 

-395 

327 

>6 4 

24 

13.2 

10.2 

582 inbound 

-0.399 

1.634 

262 

331 

£12 = 9 

16 

60 

5 4 

582 outbound 

0.542 

1.2.<10 

272 

310 (322) 

12 = 5 

19 (17) 

26.0 

189(18.1) 

582 outbound 

0.385 

1.372 

271 

299 (307) 

13 = 3 

21 (20) 

25.5 

16 8(16.3) 

582 outbound 

0418 

1.493 

266 

294 (500) 

15 = 6 

23 (221 

25.2-24 2 

16.2 (15.8) 

582 outbound 

0.447 

1 631 

254 

292 (297) 

26= 7 

24 (23) 

24 8-24 1 

15 9 (15.5) 

582 outbound 

0.470 

1.764 

251 

292 (297) 

40 = 15 

25 (24) 

20 4 

15.6(15.3) 

582 outbound 

0 490 

1.908 

304 

293 (298) 

40 = 6 

26 (25) 

17 8-16 4 

15.4 (15 0) 


Af. IS used in model. 


proceeding from the lonopause to the shock, and straddle the 
calculated values, which only increase by 305r as the same 
distance is traversed. 

5. Discussion and Summary 

The comparisons of measured and model parameters for 
flovi around Venus have shown generally good agreement 
when the observed and calculated shock locations are 
matched. There are two other important results from this 
study: first, that the bow shock locations are farther from 
Venus than e.\pected (considering the gas dynamic solutions 
appropriate to the estimated free-stream Mach numbers), 
and second, that there are trends with subsianlial disagree- 
ments between the observed and calculated plasma parame- 
ters in Venus' lonoshealh near the terminator. 

Because the shock locations near the terminator arc 
systematically farther from Venus than .predicted from the 
gas dynamic model, the model ionoshcath width was in- 
creased by lowenng the Mach number. This characteristic 
appeared in the earlier study (two orbits) of Spreiier and 
Siohara [1980t]. The resulting situation with regard to the 
model and measured Mach numbers has been given in Table 
1. Theis ei al. [1981] also found that a relatively low Mach 
number (—3.8) was appropriate for a mean location of 475 
Pioneer N'enus shock locations. The locations of the bow 
shock crossings, used to determine model gas dynamic Mach 
numbers for the five steady orbits studied here, tend to 
strjddle the locus of the published least-squares fit elliptical 
how shock shape [Slavin ct al.. 1980; sec also Slavin and 
lioizcr. 1981). The gas dynamic bow shock shapes that fit 
those crossings are more blunt than the fitted ellipse, and 
typically stand off a feu hundred kilometers higher at the 
nose, compared with the least-squares fit ellipse The great- 
est stand-off distances correspond to the lowest Mach num- 
bers. 

Possibly the bow shock is pushed away from Venus at the 
flanks as a result of mass addition to the ionoshcath flow, 
from the ionosphere. Mihalov and Barnes (1981) have given 
evidence for such mass addition al and forward of the 
icrmir.jtor plane Taylor ci al. 1 1980) have reported a layer of 
Piovking ions outside the lonopause. resulting from ‘solar 


wind acceleration' of thermal ions in the upper ionosphere. 
Mass addition might result in decreasing the tonosheath 
w'ave speeds The preliminary study of Slavin ei al. [1979] 
remarked on the possibility of ‘alteration of the magneto- 
sheath flow‘ in connection «.ith the solar wind- Venus inter- 
action. The enhanced generation of whistler-mode wave 
energy at Venus' bow shock [Scorf ei al.. 1980] could be 
pertinent, also. For the aligned MHD case, reduction of 
alone results in an increased shock distance at the flank 
[Spreiier and Rizzi. 1974], but the calculated magnitudes of 
these effects are insufficient to account for the differences 
observed by Pioneer Venus. 

The effect of magnetic stress on boundary regions in 
plasma flow around planets has been estimated by Zm o/i and 
fVo//I1976). This effect might cause discrepancies between 
observations and results from gas dynamic calculations; 
increased magnetic field strengths and reduced plasma den- 
sities near the boundary (body) are expected. Since the 
magnetic stress has been described as eliciting a squeezing 
out of plasma along magnetic field lines near the body, a 
transverse transport of momentum (a component normal to 
the undisturbed flow direction) would be introduced into the 
flow field. 

The model ionoshcath flow speeds tend to be larger than 
the observed values (18 out of 23 cases with matched 
observed and calculated shock locations). Pan of the dis- 
crepancy between observed and calculated ionoshcath flou 
speeds can be described as the difference between the 
tendencies for the calculated speed to increase and for the 
observed speed to decrease noticeably, both as the tone- 
pause IS approached. The difference between these trends 
might be due to transverse momentum flow, induced by 
magnetic stresses. .A decreased flow’ speed near the wake 
boundary has been inferred from Mariner 5 date also (see 
Pcrci-de-Tejuda [1979]. Russell [1979]. Spreiier and AlLsnc 
[1970]. and Rizzi [1971]). Mass additions near the bounoary 
might occur so as to slow the flow speed there. This effect 
has been analyzed for the case of Titan by Hurtle e; al. 
[1982]. In naive terms, the observed flow sometimes seems 
to act as if some type of boundary layer were present, so that 
the flow IS retarded The five plasma data samples with 
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observed flow speeds greater than model values arc all 
adjacent to bow shock crossings. Such examples indicate 
smaller speed jumps at the shock, that is, a weaker shock, 
than that given by the model. 

Usualh the calculated proton number density is less than 
that observed (13 out of 14 cases with matched shock 
locations, excluding orbit 582) but is within the- fitting 
uncertainty. For the outbound leg of orbit 582, the calculated 
number densities gradually increase by 40% from near the 
lonopause to the vicinity of the shock, with the observed and 
calculated shock locations matched, w-hereas the measured 
values increase by 3 1/3 times, from about 2/3 of the 
calculated value near the lonopause. to about 1 1/2 times the 
calculated value near the shock. The effect predicted from 
consideration of magnetic stress, a density depiction near 
the boundary' (lonopause). relative to values from a gas- 
dynamic model [Zwan and Wolf. 1976), is not unambig- 
uously detected in our study. 

The comrnn'^on of observed and calculated proton flux 
values has been described. The measured flux values seem 
to be slightly larger than the calculated values when the 
observed and calculated shock locations are matched, as 
de^cnbec'-e:.'lier Additionally, the orbit .^82 outbound com- 
panson seems to indicate redistribution of the observed 
lonosheath fiu.x aw-ay from the inner boundary, toward the 
shock, in comparison with the model Because the observed 
values used here tend to have been acquired near the polar 
regions, one may speculate that nonaxisymmetry has result- 
ed in an equatorial region with depleted fluxes. .MHD effects 
will introduce asymmetry, some of the limned observational 
evidence currently available for Venus has been summarized 
by RuiscU and Vaisberg (1982]. 

In principle, it is possible that some type of lime aliasing of 
the lonosheath plasma measurements could cause the mea- 
sured number densities to be too large, because the current 
measurements are peak values during vanous sampling 
limes. Then the measured flux values would be erroneously 
high To date, such an effect has not been established to be 
significant. 

The calculated magnetic field magnitudes tend to be 
appreciably lower than those observed, when the observed 
and calculated shock locations are matched. The closer 
agreement of the two sets near the ionopau.se in most of 
these cases is due to a drop of the measured values near the 
boundary, and so does not suggest the effect of the boundary 
layer due to magnetic stress, as described by Zwan and Wolf 
1 1976). an increase of the measured values near the boundary 
uould be expected as a result of that effect Instead, most of 
the companNOns suggest a substantial modification of the 
lonoNhcath flow compared with the gas-dynamic model, 
which has been used with amficially low Mach numbers 
(-5) in order to match the shock locations. For orbit 582 
(unbound, m.atchinp the observed and calculated shock 
locations has resulted in reversing the gradient of the calcu- 
lated magnetic field magnitudes from the case when the 
up^re.-.m .Mach number model was used (Table 2). so that 
the field magnitude decreases as the shock is approached, in 
.igreement unh the observations. 

Rii\\i lt i-i lit. II979J reported that Venus' bow shock is 
ucaker than the terrestrial bow shock, uith magnetic field 
magnitude jumps acro's the shock that are less than expeci- 
eu tri'm a gas d\ namic model with a Mach number of 6 The 
lesuli' I'f this Studs indic.ite greater lonoshcaih magnetic 


field magnitudes than for gas dynamic models vkith Mach 
numbers near 3. This is obviously a matter that warrants 
further study. 

The calculated and observed flow directions and magnetic 
field orientations have been compared, for cases with the 
calculated and observed shock locations matched, in an 
attempt to understand better the differences betw-een the 
measurements and the gas dynamic model. Except for orbit 
582, the observed plasma flow deflections in the east-west 
direction seem larger than the calculated deflections, the 
observed and calculated flow deflections in the north-south 
direction seem comparable. The observed east-west and 
north-south components of the inclination of the ionosheath 
magnetic field seem greater than the calculated values, 
except for orbit 129, The upstream magnetic field is most 
inclined out of the ecliptic plane for orbit 129. 

Knudsen ci al. [1982] have presented a new method for 
calculating a shape of Venus' lonopause that agrees well 
with observationx This method uses Pioneer Venus results 
for the dependence of lonosphenc density with solar zenith 
angle and for the altitude dependence of ionospheric electron 
and ion temperature ,A converging lonopause shape down- 
stream from N'enux is achicxed b> maintaining a iruns'e'-se 
ionosheath pressure near the terminator: the flow down- 
stream from the terminator is described by two-dimensional 
supersonic expansion around a corner of an ideal gas [Liep- 
niann and RoshLo. 19.^7). .A problem with this representation 
is that the empirical ionosphere temperature and density 
model cannot balance a free-stream pressure greater than 4 
X 10"** djnes cm"* \Knud.sen ct at.. 1982). whereas the 
observed free-stream proton pressure for only two (orbits 63 
and 582) of the five orbits used in our study was less than (-3 
X lO"** dynes cm"') this limiting value. The low jMach 
number case (Af = 3. hlR„. = 0.1) lonopause shape given by 
Knudsen et al. was used in the gas dynamic calculation for 
the orbit 582 outbound conditions (low free-stream pres- 
sure). to examine the changes of the calculated values The 
Mach number required to fit the observed shock location had 
to be reduced from 2 4 to 2.3. The results are given in 
parentheses in Table 4 Except near the lonopause. where 
the maximum changes from the case with the ionopause 
derived from (5) are a 4% increase in flow speed, a 7% 
decrease in number density and a 4 % decrease in magnetic- 
field magnitude, the changes in the calculated parameters are 
small. The changes near the Knudsen ionopause might be 
expected because that surface is contracted around the 
planet (low-er altitude ionopause), in comparison with the 
one denied from (5). 

In summary, the gas dynamic model for solar wind flow 
around Venus goes results that agree fairh well with the 
observed flow' velocities, proton number densities and iso- 
tropic temperatures, and vector magnetic fields, for five 
orbits with stead) flow, when the calculatec and obverved 
shock locations are matched The disagreements that have 
been found — principally a tendenev for the flow field to 
exhibit hicner number denviiie'- and magnetic field sircngihs 
than the model, together wuh the more distant shock loca- 
tions than expected — vuggevt vome significant differences 
between the actual and model flow fields Some of the 
differences ma\' be due to the selection of the anificialK low' 
Mach numbers for the calculations in order to match the 
observed shock loc.itions The comparisons all are in the 
vicinity of the lermin.tior The ionosheath seems to h.oc 
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slower speeds and be more ’compressed’ than the model, 
with greater proton densities and magnetic field strengths, 
when the low jMach numbers are used to match the shock 
locations The data also appear to provide some evidence for 
further retardation of the flow speed near the body (iono- 
pause). This may be indicative of the presence of a boundary 
layer, or the incorporation of ionospheric material in the 
flow, or both. Finally, it should be noted that since the 
comparisons of this study are for the vicinity of the termina- 
tor, some of the conclusions might be altered if conditions 
elsewhere were examined similarly. 
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Solar Wind Flow About the Terrestrial Planets 
2. Comparison with Gas Dynamic Theory and Implications 
for Solar-Planetary Interactions 

J. A. Slavin,' R. E. Holzer,' J. R. Spreiter,^ S. S. Stahara/ and D. S. Chaussee^ 


This study utilizes gas dynamic calculations in conjunction with observational bow shock models lo 
investigate the solar wind flow patterns about the terrestrial planets. .Average dayside bow shock 
position could be predicted for the earth by theory with an error of only ~29r, given the observed 
shape and location of the magnetopause. Accordingly, our findings confirm the validity of the single- 
fluid gas dynamic approximation for describing this major aspect of solar wind flow past the eanh. 
Modeled using gas dynamic theory, the solar wind interactions with Venus and Mars exhibit very 
significant differences. At Mars the mean inferred altitude of the solar wind-obstacle interface vanes 
from 510 km at the stagnation point to almost 1000 km near the terminator. The effective magnetic 
moment required to produce a magnetosphere of this size for average solar wind dynamic pressures 
and terrestnal-type internal current systems is 1.4 ± 0.6 x 10^* G cm^. Gas dynamic modeling of the 
January 21, 1972. .Mars 3 and July 20. 1976, Viking 1 lander particles and fields observations supports 
the conclusion that the Martian obstacle to the solar wind lies at altitudes too high for it to be 
associated with only an lonosphenc or atmospheric interaction. In contrast with .Mars, our modeling of 
the Venus observations has found that the bow wave is closer to the planet than would be expected for 
a purely ionospheric obstacle. The subsoiar width of the Venus lonosheaih in the Venera and PVO 
measurements is only 60% and 90%, respectively, of that predicted by the gas dynamic model. This 
result IS attributed to the presence of solar wind-neutral atmosphere interactions in the lower 
ionosheath that are not included in the gas dynamic code. 


Introduction 

In this article we report upon the second of a two-part 
investigation of the solar wind flow about the terrestrial 
planets. The first portion was presented in the paper by 
Slavin and Holzer [1981], which we will henceforth refer to 
as paper 1. In that study, bow wave observations at Venus, 
Earth, and Mars were brought together and modeled by 
using a standard set of methods for determining and compar- 
ing their shock shapes and positions. Here we continue by 
employing the bow shock models of paper 1 as boundary 
conditions for gas dynamic calculations from which the flow 
fields about these bodies are constructed. 

The principal objectives are threefold. First, the large data 
base available for the earth system is compared with the 
predictions of gas dynamic theory in order to test the method 
before its use at the other planets, where far fewer observa- 
tions have been made. Second, the gas dynamic calculations 
are applied to the Pioneer Venus orbiter (PVO) and Venera 9 
and 10 shock measurements with the intent of evaluating the 
degree to which flow abom V^nus may be modeled without 
including the effects of solar wind plasma-neutral atmo- 
sphere interactions above the lonopause. While the PVO 
observations have shown the Venus interaction to be primar- 
ily ionospheric [Brace et al., 1980; Elphic et al., 1980; Brace 
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et al., 1982], cometary processes (i.e., charge exchange and 
photoion pickup) in the ionosheath can exert a great influ- 
ence over the flow [Wallis, 1972; 1973], While passing 
through the inner solar system, most active comets are 
thought to deflect the solar wind at large distances through 
charge exchange and photoion-pickup-induced mass loading 
of the incident solar wind flow [Biermann et al.. 1967]. 
Clarification of the role played by the neutral atmosphere in 
the Cytherean interaction is cntical to the use of Venus as an 
analog of the solar wind flow about comets. Third, at Mars 
we use gas dynamic models to infer the averge height of the 
solar wind-obstacle interface and evaluate the controversial 
January 21, 1972, Mars 3 observations that show a possible 
encounter with an intrinsic field magnetosphere [e.g., Dol- 
ginov, 1978c, b, c; Russell, 1978c, b]. The question of 
whether Mars undergoes a Venus-type ionospheric interac- 
tion with the solar wind or possesses a significant planetary 
magnetic moment remains unresolved because of a lack of 
low-aJtitude particles and fields measurements. Until such 
time as new observations are made by an orbiter with a low 
periapsis, siTch arthe prbpose^d ESA Kapler Mars missibrir 
the Martian solar wind interaction must be studied predomi- 
nantly through modeling of the high-altitude (i.e., >1000 km) 
Soviet Mars satellite measurements. 

Gas Dynamic Theory 

The methodology and computational techniques employed 
in the gas dynamic modeling of flow about the planets are 
well documented in the literature [e.g., Spreiter et al., 1966; 
Dryer and Heckman, 1967; Spreiter and Alksne, 1970; 
Dryer, 1970; Rizzi, 1971; Spreiter and Stahara, 1980c, b]. 
For purposes of this study it is necessary for us to consider 
only the general properties of the solutions. Figure 1 displays 
bow shock position, streamlines, and charactenstic lines for 
a single-fluid gas dynamic calculation of flow about a theo- 
retical magnetopause surface [Spreiter et al., 1966], Far 
upstream the solar wind is assumed to be an ideal MHD 
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Fig. I . Gas dynamic model of flow about a theoretical magneto- 
pause for /V/, = 8. and an adiabatic exponent of 2 from Spreiier et al. 
[I966J IS displayed. Aberrated geocentnc solar ecliptic coordinates 
are used as in papier I . 

fluid. It is completely descnbed in the MHD equations by a 
sonic Mach number an Alfvenic Mach number Ma, and 
an adiabatic index y for the equation of state. Previous 
theoretical examinations of this MHD problem (e.g., 
Spretter and Rizzi, 1974; Chao and Wiskerclien, 1974] have 
shown that the bow shock jump conditions and the flow 
solutions are dependent separately upon Ms and Ma without 
any manipulation of these terms into a single Mach number 
being possible (although magnetosonic Mach number and /3 
may be used instead of M, and Ma; Tidman and Krall 
[1971]). For this reason there does not appear to be any 
theoretical support for the use of magnetosonic Mach num- 
bers in studies of the dayside magnetosheath flow and bow 
shock position, as is sometimes done on intuitive grounds 
[e.g.. Russell, I977h]. In some cases the desire to use Mach 
numbers based upon the MHD fast-mode velocity appears to 
stem from the need to simplify the problem by making the 
Mach number that determines the distant downstream Mach 
cone the only Mach number in the flow calculations [e.g.. 
Dryer and Heckman, 1967; Shen, 1971]. These dififlculties 
are understandable when it is realized that no general 
theoretical solutions to this MHD problem have been found. 
Only in the specialized case of aligned flow (i.e.. V||B) has 
the problem of ideal MHD flow about an axisymmetric, 
blunt, magnetospheric obstacle been solved [Rizit, 1971; 
Spreiier and Rizzi. 1974]. In that particular situation the 
MHD equations can be reduced to those of gas dynamics 
with a novel equation of state and solved numencally. 
However, the limit in which this MHD problem may general- 
ly be solved is that of a weak magnetic field in the sense of 
Ma* = (pV-l2)IIB'l%-n) & 10^. This limit correspx>nds to the 
traditional ‘gasdynamic approximation’ [e.g.. Spreiier et al.. 
1966. Dryer and Heckman. 1967] in which the terms of the 


MHD equations containing Ma are neglected on the grounds 
of their relative smallness. Solutions such as show n in Figure 
I are then obtained. Only in very recent times have any new 
attempts been made on the full .MHD descnption \L\on et 
a!.. 1980; Leboeuf et al., 1981; Wu et al.. 1981]. It is not yet 
clear when, or to what extent, these computational plasma 
physics codes will begin to contribute to our understanding 
of astrophysical plasma flow systems. 

While the interplanetary magnetic field plays essential 
roles m almost all of the plasma processes taking place in the 
solar system (e.g.. particle acceleration and transport), the 
magnetic energy in the solar wind is weak in companson to 
the flow energy, and the gas dynamic approximation is 
usually valid (see. for e.xample. Table 5 of paper I ). Still, the 
neglect of the IMF is widely felt to be the major limitation of 
the gas dynamic approach. It is not often appreciated that 
the more cntical assumption in these calculations may be the 
use of a single-fluid description. Not only are the solar wind 
protons and electrons observed to have significantly differ- 
ent mean thermal speeds but their distnbution functions and 
the processes influencing them are quite dissimilar [e.g.. 
Sinter and Scudder, 1980]. Under these circumstances it is 
not physically reasonable to expect any single-fluid treat- 
ment to be able to predict the individual properties of these 
two difl'erent plasma species. However, it is possible for a 
single-fluid model to well represent the common, or shared, 
properties of a multispecies system. In this respect the 
problem of flow past a planet is like that of the expansion of 
the solar wind itself. Single-fluid models can be constructed 
to match the bulk speed and density at 1 AU but not the 
individual proton and electron temperatures (e.g . Htmdhau- 
sen, 1972]. For this reason, in addition to introducing 
extended energy sources and other transport mechanisms, 
modem treatments of the solar wind often consider the 
protons and electrons as separate fluids. However, just as 
with the gas dynamic calculations in Figure 1. the magnetic 
field in the solar wind is still often regarded as frozen-in and 
dynamically unimportant once a given parcel of plasma is far 
from the initial acceleration region. 

The need for a two-fluid model becomes even more 
apparent in the measurements across the shock While the 
electrons and ions are shocked at about the same location 
relative to fluid length scales, the microphysics of the 
thermalization processes are different for the two species 
[e.g., Greenstadt and Fredricks, 1979]. The result is that the 
electron temperature increases by less than a factor of 2 
[e.g., Bame et al., 1979] as contrasted with ajump 2-4 times 
as great for the protons [e.g., Montgomery et al., 1970]. 
Accordingly, the present gas dynamic models do no more 
than approximate the prediction of the individual species 
temperatures [e.g., Mihalov et al.. 1980. 1982]. With respect 
to the validity of the gas dynamic approach it is important to 
note that the model successfully predicts all of the shared 
properties: shock location [e.g., Wolfe et al.. 1966; this 
paper], velocity jump [e.g., Burlaga and Ogilvie, 1968; 
Mihalov et al., 1980). density jump [e.g.. Howe and Binsack. 
1972], and total enthalpy [Montgomery et al.. 1970]. Hence, 
there is ample justification for the use of gas dynamic 
modeling as a tool in studying solar-planetary interactions 
once Its limitations are recognized. 

As discussed above, sonic and magnetosonic Mach num- 
bers are the most commonly assumed ordering parameters 
for studies of bow shock location. Mach numbers based 
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upon the Alfven speed are also used sometimes for unstated 
reasons (e.g., Harvey et al.. 1981]. ObservaiionaJly, the- 
determination of the appropriate Mach numbers for this 
problem is not straightforward for a number of reasons. The 
high averge values and small spread in and [e.g., 
Fairfield, 1971 ; Formisano, 1979] and the dominance of solar 
wind dynamic pressure m setting bow wave position [Fair- 
field. 1971] make the variations in shock position and shape 
with changing Mach number difficult to detect. The lack of a 
.Mach-number-independent means of determining the adia- 
batic e.xponent [Dryer, 1971] also introduces ambiguities into 
the application of gas dynamic theory to the problem. In this 
study the gas dynamic approximation is found to produce an 
accurate representation of the terrestrial bow shock/magne- 
topause system with a sonic Mach number 

M, = Ml-08 7p -H 1.14 r,)/1.16 (I) 

based upon total plasma pressure [including 4% He*'*' with 
ThJTh = 3.5, see Table 5 of paper 1] and an adiabatic 
exponent of 2. Not only does the gas dynamic model that 
uses M, computed in this way reproduce the average obser- 
vations but A/j is also found to order the shock shapes and 
stand-off distances determined with the various earth-orbit- 
ing satellites in paper 1. 

The effects of the choice of a -y value appear in the model 
jump conditions across the shock and in the resulting width , 
of the magnetosheath. Figure 2 shows both the y = Sfi and 2 
bow wave positions from Spreiter et al. [1966]. In this study 
we found that for a Mach number defined as in (1) an 
adiabatic exponent of 2 produces good agreement between 
gas dynamic theory and observation, although an optimum 
value might be very slightly lower. This result is in essential 
agreement with the VELA studies of Gosling et al. [1967] 
and Montgomery et al. [1970], the IMP investigation of 
Fairfield [1971], and the ISEE work of Zhuang and Russell 
[1981]. While a more detailed examination of this issue lies 
beyond the scope of our present investigation, further study 
of the problem is planned for the near future. 

The variation of shock shape with obstacle shape is also 
addressed in Figure 2. Blunter bodies require thicker transi- 
tion regions (i.e., more distant bow shocks) to accomplish 
the diversion of the incident flow than do more slender 
obstacles. It was, in fact, the one-to-one relationship be- 
tween shock and body shape that formed the physical basis 
of the original inverse method for the calculation of hyper- 
sonic flow about obstacles [Van Dyke, 1958]. Shock shape 
would be iteratively varied until the resulting stagnation 
streamline coincided with the desired body shape. As shown 
in Figure 22 of paper 1 . the bow waves of Venus, Earth, and 
Mars all possess different shapes in good qualitative agree- 
ment with their observed obstacle shapes and the aerody- 
namic relationships established by Van Dyke. For this 
reason, studies of flow about the planets require theoretical 
models with obstacle shape as an input variable. Fortunate- 
ly, such codes now exist as a result of the work of Spreiter et 
al. [1970], Chausee et al. [1978], Spreiter and Stahara 
[1980a. b], and Stahara et al. [1980] on the solar wind 
interaction with Venus. Lacking, in situ observations of 
ionopause shape, they performed pressure balance calcula- 
tions with the ratio of ionospheric scale height to subsolar 
obstacle radius H/Rq as a vanable to obtain ionopause shape. 
The larger the assumed scale height the greater the bluntness 
of their tangential discontinuity obstacle. The increased 



Fig. 2. Bow shock positions associated with theoretical earth 
(i.e.. H/Ro = 0.20) and Venus (i.e., H/Rn = 0 01) obstacles for ,Vf, = 
8 and y = 2 and 5/3 are compared. Note the increased width m the 
region between the shock and obstacle for blunter obstacles and 
larger adiabatic exponents. 

bluntness is due to the large decline in external pressure 
between the nose and terminator coupled with the iono- 
sphere's exponential pressure dependence upon scale 
height. The resulting flow models were then cataloged by 
A/„ y, and H/Ro [Stahara et al.. 1980]. Hence gas dynamic 
models are now available for a family of blunt obstacles 
parameterized by HIRq. 

Modeling Results 

Earth 

In paper 1 the location of the terrestnal bow shock was 
modeled independently with the observations of IMP 3, IMP 
4, Heos 1, and Prognoz 1 and 2. The results of that study are 
displayed in Figure 3, where they have been scaled, assum- 
ing the usual sixth root dipole-external pressure relation 
(e.g., equation (2)), to a common solar wind dynamic pres- 
sure of 1.16 = 2.1 X 10"® dynes/cm*. As shown, the 

differences remaining among the four surfaces, while small, 
appear well ordered by the average observed upstream sonic 
Mach numbers. Both nose distance and eccentricity of the 
model surfaces vary monitonically with M^. Such an order- 
ing is not present when the models are compared with the 
average Alfvenic and magnetosonic Mach numbers listed in 
Table 5 of paper 1. Given the high average Alfvenic Mach 
number in the solar wind this result is in qualitative agree- 
ment with expectations based upon the gas dynamic approxi- 
mation discussed earlier. 

Figure 4 compares the average bow wave from paper 1 
with the gas dynamic models of flow past a theoretical 
magnetopause surface by Spreiter et al. [1966]. Because it is 
our intent to predict the flow pattern behind the shock at 
Venus and Mars, we take this same approach here for the 
earth. In Figure 4 the forward ponions of the theoretical and 
observed bow waves have been scaled to coincide. Magneto- 
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Fig. 3 The average shock surfaces from five eanh orbiter 
missions obtained m paper 1 are displayed after being scaled to a 
common mean solar wind dynamic pressure. 


EARTH BOW SHOCK AND MAGNETOPAUSE 
(ALL SCALED TO Pg^= 2.1x10'® dynes/cm^) 



Fig. 4 The average bow wave from paper I and the magneto- 
pause from Hotter and Slavtn (1978) are compared with the Spreiler 
el al. {1966) gas dynamic calculations. 


pause position is then the predicted quantity. The observed 
magnetopause model with which the results are compared 
is the average model of Holier and Slavin [1978) scaled to 
P,w = 2.1 X 10"“ dynesfcm’. This model surface was chosen 
because the study producing it considered the effects of solar 
wind dynamic pressure, and the observations upon which it 
is based were made near the center of the interval spanned 
by this investigation. In this way any effect of the solar cycle 
variation in geomagnetic activity on magnetopause position 
should be minimized. The shape of Holzer and Slavin 
magnetopause model is close to that of Fairfield [1971] but 
slightly less blunt than some of the others, such as Howe and 
Binsack [1972]. The gas dynamic bow shock/magnetopause 
positions displayed correspond to a sonic Mach number of 8 
and adiabatic indices of 2 and 5/3. A sonic Mach number of 
7.2, the mean in Figure 3. would have been preferred, but 
funding constraints limited us to using previously generated 
models in most cases. The error introduced by the use of 
iV/j = 8 will be much less than the spread among the surface 
models in Figure 3. The spread in Figure 3 is larger than 
would be predicted by gas dynamic theory: this may be a 
result of the fitting uncertainties outlined in paper I. Howev- 
er, the more important factors in Figure 4 are yand obstacle 
shape. A value of y = 5/3 is seen to predict a magnetosheath 
that IS too narrow, whereas y = 2 appears to be in near 
agreement with observation. The greatest discrepancy be- 
tween the theoretical model and observational reality lies in 
their respective magnetopause shapes. This is a general 
failing of the various theoretical magnetopause calculations 
[e.g., Fairfield, \^1\\ Walker, 1979] believed to be caused by 
inaccuracies in their specification of the external pressure 
field and/or the internal magnetosphenc current systems. 
Fluid calculations that use theoretical obstacle shape then 
produce a bow wave that is too slender and lies too close to 
the planet. We avoid this difficulty by selecting an obstacle 
surface in Figure 5 that is very near the average shape of the 
magnetopause as determined by actual observation. As 
shown, a HIRq = 0.35 obstacle is nearly midway between the 
slightly more slender magnetopause models of Fairfield and 
Holzer and Slavin and the blunter result of Howe and 
Binsack. Fanher downstream the flat slope of the HIRa = 
0.35 obstacle could become a source of error for compari- 
sons with the observations, but we are concerned here with 
only the near portions of the bow wave. 



Fig. 5. A comparison of the empincal Fairfield [\91\ \ and Howe 
and Binsack [1972] magnetopause models with ihe gas dynamic 
HIR(, = 0.35 obstacle demonstrates the appropnateness of this 
theoretical surface for flow calculations. 




Slavin et al.: Solar Wind Flow About the Planets 


In Figure 6 this HIRq = 0.35 obstacle has been used to 
calculate the flow pattern about the earth's magnetosphere. 
In order to compare the assumed obstacle shape with the 
Holzer and Slavin magnetopause, the radius of the obstacle 
has been chosen so that they share the same stagnation 
points. The agreement over the forward region of the magne- 
topause influencing the shock for the interval of interest (see 
the characteristic lines of the flow in Figure 1) is excellent. 
As expected from the discussion of the preceeding figure, 
the assumed obstacle is slightly blunter along the flanks than 
the Holzer and Slavin observational model but less so than 
Howe and Binsack. The resulting theoretical bow shock 
surface is in close agreement with paper I’s observational 
model as shown in the figure. The discrepancies near the 
nose and at the terminator amount to only —2%. Error bars 
normal to the empincal shock and magnetopause surfaces 
provide a measure of the rms deviations of the observed 
boundary crossings about the model. As indicated, the 
differences between theory and observation appear much 
smaller than the rms spread in the data themselves and only 
somewhat larger than the fitting uncertainties discussed in 
paper 1. If the small overestimate in stand-off distance is 
real, then the optimum value of the adiabatic exponent may 
be slightly less than 2, as suggested by Fairfield. However, it 
is also possible that the small differences between the 
predicted and observed bow waves may be an MHD effect 
arising from the omission of the IMF in the gas dynamic 
calculations. For the special case of aligned flow, Spreiter 
and Rizzi [1974] have found that at Alfvenic Mach numbers 
greater than 5 (cf., the mean value for the observations in 
Figure 3 is 9.0; see Table 5 of paper 1) the shock surface is 
only slightly blunter than would be expected from gas 

EARTH BOW SHOCK AND MAGNETOPAUSE 

(ALL SCALED TO Pgw=2.1xlO-® dynes/cm2) 



Fig. 6. The average observational boundary surfaces from Fig- 
ure 5 are compared above with gas dynamic flow past a HIRq = 0.35 
obstacle for y = 2 and Af, = 8 upstream conditions. The error bars 
provide a measure of the rms spread in the spacecraft boundary 
crossings modeled in paper I. 



Fig. 7. The aligned MHD flow model of Spreiter and Rizzi 
[1974] IS shown above. Future studies of low conditions, using a 
larger data base than employed by this investigation, are planned to 
test these theoretical predictions. 

dynamic calculations as detailed in Figure 7. Unfortunately, 
a much larger data set than the one employed by this study 
would be required to venfy the results of their MHD model, 
given the modest rate of occurrence of low-A/^ conditions. 

Venus 

The Cytherean bow wave has now been surveyed from 
orbit by Venera 9, Venera 10. and Pioneer Venus, as 
discussed in paper 1 [see also Verigin et al., 1978; Slavin et 
al., \9S0', Smirnov et al., 19S0\ Russell and Vaisberg, 1982]. 
Since the initial measurements of the Manner 5 flyby and the 
Venera 4 descent probe in 1967 [Bridge et al., 1967; Dol- 
ginov et al., 1969; Cringauz et al.. 1968], a number of studies 
have compared bow wave position with the predictions of 
gas dynamics [Gnngauz et al., 1970; Spreiter et al., 1970. 
Rizzi, 1971; Vaisberg and Bogdanov, \91A\ Ness et al., 1974; 
Bridge et al., 1974; Wallis, 1974; Vaisberg et al., 1976a,- 
Cringauz et al., 1976; Russell, 1977a, b; Dolgtnov et al., 
\91%, Romanov et al., 1978; Verigin etal., \9H,Wolfe et al., 
1979; Slavin et al., 1979a, b; Mihalov et al., 1980; Slavin et 
al., \9%Q\ Smirnov et al., 19S0: Spreiter and Stahara, 1980b; 
Taylor et al., \9%Q,Theis et al., 1981]. General cnticisms that 
can be leveled at many of these works may be listed under 
three main categones; they generally do not uke into 
account measured ionopause and shock shape; the adiabatic 
exponent is commonly chosen to be 5/3 as opposed to the 
preferred (i.e., from the far more numerous earth observa- 
tions) value of 2; and finally the possibility of nonideal gas 
dynamic processes such as charge exchange and photoion 
pickup influencing the flow is often ignored. With these 
points in mind the approach adopted here has been to match 
the observed shape of the bow wave with a gas dynamic 
model and then compare the predicted obstacle surface with 
existing ionopause measurements in the same way as Figure 
4 does for the earth’s magnetopause. An additional advan- 
tage to this method is that it predicts the solar wind flow field 



Slavin et al . Solar Wind Flow About the Planets 



Fig. 8. The Venera 9 and 10. Pioneer Venus, and Mars 2. 3. and 
5 bow shock crossings from paper 1 have all been mapped to zero- 
aberraied solar zenith angle by means of their average surface 
models in paper I. The Viking I lander bow shock crossing has also 
been included by using the mean boundary shape derived from the 
orbiter measurements. In the upper left-hand comer the average 
height of the Venus subsolar lonopause has been marked to separate 
the attached and detached bow shock regimes. 


from Venera and Pioneer Venus. The solutions are a func- 
tion of the shape of the obstacle (e.g.. W//?o = 0.20 corre- 
sponds to MP = classical Chapman- Ferraro magnetopause), 
the adiabatic exponent (5/3 or 2). the sonic Mach number, 
g and the obstacle radius, which has been chosen so as to 
^ make the noses of the theoretical and observed bow shocks 

C/) 

c coincide. As in the terrestrial case a sonic Mach number of 8 
^ was used, but again the error introduced should not be large 
« in this hypersonic regime. The net effect of using 8 at Venus, 
g as opposed to the = 6-7 expected to be more tvpical of 
<£ 0.7 AD (e.g., see Table 2 of paper 1). is to slightly underesti- 

H mate the width of the ionosheath. For both the Venera and 
E PVO observations it was found that an HIRo = 0.20 obstacle 
2 shape provided the best fit to the measured bow shock 
shape. In the subsolar region the observed distance between 
the shock and lonopause from the Venera and PVO measure- 
ments is only 60% and 90%, respectively, of that predicted 
by gas dynamics with y = 2 For Venera, in fact, the 
effective obstacle surface passes beneath the surface of the 
planet. Further, the flanks of the gas dynamic obstacle 
needed to explain the shape of the shock flare out more than 
IS observed for the Venus ionopause [Brace et al., 1980], 
Based upon the ionosphenc scale heights measured by PVO 
[Nanle et al.. 1980], only an H/Ro = 0.01 ionopause shape 
would have been expected. In the section devoted to solar- 
planetary interactions these findings will be interpreted in 
terms of significant cometary interactions in the low-altitude 


behind the bow shock. Unlike at the eanh, highly eccentnc 
orbiters at Venus must pass through the low-altitude dayside 
regions at speeds sufficiently great to alias the plasma 
measurements [Mihalov et al.. 1980. 1982]. For this reason it 
IS highly desirable to be able to infer the flow field near the 
ionopause with a method that is independent of the low- 
altitude observations. 

The Venera and PVO bow wave models are considered 
separately, owing to the very different altitudes at which 
they place the shock, as shown in Figure 8. The possible 
reasons for the difference in bow shock position between 
these missions, which we have discussed at length in paper I 
and Slavin et al. [19796. 1980], are still the subject of 
controversy. Smirnov et al. [1980, 1981] and Romanov et al. 
[1978] favor an asymmetnc shock cross section at Venus 
with the nonaxially symmetric IMF as the underlying cause. 
The near-polar-orbiting Pioneer Venus spacecraft then ob- 
serves a more distant bow wave because the shock resides at 
higher altitudes over the poles than near the ecliptic plane 
where Venera 9/10 encountered it. However, as detailed in 
paper 1 and our earlier studies, we have found the trace of 
the shock surface in the aberrated terminator plane to be 
approximately circular at both the earth and Venus, indepen- 
dent of IMF orientation. Further, no anomalously large 
aberrations of the symmetry axis of the bow wave away from 
the solar wind flow direction have been observed. For this 
reason we favor a temporal cause for the discrepancy 
between Venera and PVO. with significant changes in the 
low-altitude solar wind interaction taking place between 
solar cycle minimum when the Venera measurements were 
taken and the solar maximum epoch of the PVO pnmary 
mission. 

In Figures 9 and 10 a number of theoretical bow wave 
solutions are compared with the observed shock positions 



X- (Ry) 

Fig. 9. The average Venera 9 and 10 bow shock model from 
paper I is compared with several gas dynamic models. They are 
labeled by calculation type (i.e., GD = gas dynamic), obstacle type 
(i.e., IP = ionopause; MP = magnetopause). HtRn value if IP, and 
the adiabatic exponent. Sonic Mach number appears at the top of 
the figure. 




Slavin et al.. Solar Wind Flow About the Planets 


ionosheath. It will be suggested that at least for some 
modeling applications these processes may be regarded as 
decreasing the adiabatic exponent behind the shock from the 
value of 2 appropriate to planets without large neutral 
populations extending into the solar wind. The term ‘lono- 
sheath’ may then imply the presence of plasma-neutral 
atmosphere interactions as opposed to simply providing an 
alternate name for a magnetosheath surrounding a tangential 
discontinuity supported by an ionosphere. Finally, the evi- 
dence from all sources for such an lonospheric-cometary 
model of the Venus solar wind interaction will be reviewed. 

Mars 

In Figure 1 1 the Mars bow wave model surface from paper 
1 is compared with a number of gas dynamic flow solutions 
[e.g., Spreiter and Stahara, 1980a. b\ Stahara ei al., 1980) in 
the same manner that was employed at the earth and Venus. 
The best fit to the shape of the bow shock, indicated with a 
solid line, was produced with = 7.4, y = 2. and H!R^ = 
0.03. A sonic Mach number of 7.4 is about 20% lower than 
the expected mean at 1.5 AU. However, hypersonic flow in 
the near field fi.e., .r' > -1 is relatively insensitive to 
changing Mach number beyond Ms = 5. The use of 7.4 
because of its availability [e.g., Spreiter and Stahara, 1980f>] 
should not be a significant source of error in the model 
calculation. In marked contrast with the situation at Venus, 
the best fit model in Figure 1 1 indicates that the average 
Martian obstacle to the solar wind ranges from 510 (±20%) 
at the subsolar point to nearly 1000 km near the terminator. 
The assumption of y = 5/3, instead of y = 2. would place the 
solar wind-obstacle interface at still greater altitudes. How- 
ever, the small role of the neutral atmosphere in the solar 

pvo Venus 



Fig. 10. The average Pioneer Venus bow shock model from paper 
I IS compared with a number of gas dynamic models, as m Figure 9. 



Fig. 1 1 The average Mars bow shock model from paper I is 
compared with a number of gas dynamic models, as in the two 
preceding figures. 


wind interaction with Mars at altitudes above 3(X) km [Clou- 
tier et at., 1969) argues in favor of using the adiabatic 
exponent that best characterizes the flow about the earth as 
opposed to Venus. As treated by Slavm and Holier [1982], 
the height of the Martian obstacle in Figure 1 1 appears to be 
far too great for this planet to be undergoing a Venus-type 
interaction with the solar wind. In a later section we will 
examine additional Mars orbiter observations with the aid of 
the theoretical flow models developed in the course of this 
study. 

• This modeling of the Mars bow wave with a complete gas 
dynamic flow field represents an improvement over many of 
the earlier studies [e.g.. Russell, 1977c). In most of those 
works the shock crossings were modeled empincally to 
determine the subsolar height of the bow wave. The Martian 
magnetosheath was then assumed to have the same relative 
thickness as that of the earth in order to infer the altitude of 
the stagnation point on the obstacle. However, the width of 
the magneiosheath is a function of obstacle shape, with 
thicker magnetosheath regions being necessary to redirect 
flow about blunter bodies [Van Dyke, 1958: Spreiter and 
Stahara, 1980i). The model in Figure 11 shows an obstacle 
that is more slender than that of the earth. For this reason, 
previous investigations assuming the same relative magneto- 
sheath thickness as the terrestnal case underestimate the 
height of the solar wind-obstacle interface by about 20%. By 
comparison with the other Mars studies that have used gas 
dynamic flow calculations [Dryer and Heckman, 1967; 
Spreiter and Rizzi, 1972; Gringauz et al., 1973; Bogdanov 
and Vaisberg, 1975; Breus and Verigin, 1976; Gringauz ei 
at., 1976; Dolginov et al., 1976: Vaisberg et at., 1976i), we 
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have enjoyed the advantage of having a much larger number 
of flow models with which to work. The result has been the 
good fit to the shape of the bow wave in Figure 1 1 and hence 
the near field flow pattern In addition, we have based our 
modeling on only the adequately sampled portion of the 
Mars shock forward of ,r'> - I In this way the spatial 
biases inherent in the sparse downstream observations (e.g. 
see paper 1) have been avoided. 

Discussion 

Earth 

This study represents the first major comparison between 
gas dynamic theory and the position of the teirestnal bow 
shock since the Fairfield [1971] work conducted a decade 
ago. A detailed discussion of the data set and modeling 
techniques utilized by Fairfield is contained in paper 1. For 
the purposes of this investigation the main advantages of our 
observational shock model over that of Fairfield lie in the 
slightly larger fraction of a solar cycle sampled, 1965-1972 as 
contrasted with 1963-1968, and the greater availability of 
upstream parameters for the individual bow wave encoun- 
ters. The shape and location of these two average bow shock 
models are, in fact, quite similar (see Figure 19 of paper 1). 

In Fairfield the observational model of the shock was 
compared with the gas dynamic calculations of Spreiter and 
Jones [1963], which assumed 7 = 2 and = 8.7. Just as in 
this study, the theoretical magnetopause and bow shock 
shapes were determined to be slightly less blunt than the 
observations indicated. However, the Spreiter and Stahara 
codes now available allow for the choice of obstacle shape as 
an input parameter. Figure 6 displayed an excellent fit to the 
empirical model with HIRn = 0.35, Af, = 8, and y = 2. 

Accordingly, we concluded that our investigation has 
confirmed the good agreement between single-fluid gas 
dynamic theory and average dayside bow wave location 
obtained by Fairfield and many of the earlier works [e.g.. 
Gosling et al., 1967], It is also important to note that both the 
Fairfield study and this investigation have found that the 
single-fluid gas dynamic models perform well when a sonic 
Mach number based upon total ion and electron plasma 
pressure and an adiabatic exponent close to 2 are used. 
Finally, we have found that the average shock position and 
shape models obtained from the individual missions of paper 
1 are well ordered by A/, and not by magnetosonic or 
Alfvenic Mach numbers. All of these results support the 
validity of the gas dynamic approximation and the general 
applicability of fluid theory to the macroscopic descnptiori of 
solar wind flow past planetary size tangential discontinuity 
surfaces. Future investigations with a much larger data set 
are planned to examine quantitatively the gas dynamic 
predictions concerning the variation in shock position/shape 
with AJ,; to study the limitations of the gas dynamic calcula- 
tions farther downstream, where the position of the shock 
should be under the influence of the magnetosonic fast wave 
Mach cone; and. last, to document the expected gradual 
failure of these gas dynamic models as IMF strength in- 
creases and A/a decreases. 

Venus 

The problem of flow into and about a partially absorbing 
body in an astrophysical application was considered long ago 
in reference to the moon [Tozer and Wilson, 1967: Sonett 
and Colburn. 1968]. While there was little hope that the lunar 


surface material electncal conductivity could be high enough 
to provide for anything like the nearly total exclusion of the 
incident solar wind possible for a planetary ionosphere [e.g.. 
Cloutier and Daniell, 1973; Zeleny and Vaisberg. 1981], it 
was thought that sufficient deflection might occur to cause 
the formation of a bow shock. Under these circumstances 
the bow wave would remain strong in the sense that 7 and 
upstream Mach number would determine the jump condi- 
tions, but it would move closer to the moon as the total 
amount of absorption increased. A sufficiently high rate of 
absorption could then result in the bow shock becoming 
‘attached’ to the body with only ‘limb’ shocks remaining. It 
IS now known that the moon absorbs almost all of the 
incident solar wind and thus prevents the development of a 
global shock wave [Schubert and Lichtenstein. 1974; 
Spreiter et al., 1970] Russell [1977a] raised the possibility of 
a partial absorption of the solar wind at Venus in his 
interpretation of the Manner 10 shock observations [Bridge 
et al., 1974; Ness et al., 1974]. In Russell's controversial 
analysis [/Vexj, 1977; Russell, 1977f>], evidence that the 
Venus bow wave had been attached to the ionosphere during 
the time of the Manner 10 flyby was presented. He estimated 
that the fraction of the incident flow being absorbed at that 
time was 33%. In Figure 8 of this study, both the Venera 9 
and 10 and Pioneer Venus bow shock encounters from paper 
1 were mapped to the nose position along the x' axis. The 
nominal top of the ionosphere is marked in the upper left- 
hand comer of the figure to separate the attached and 
detached regimes. None of the Pioneer Venus events, and 
only a statistically insignificant number of the Venera shock 
encounters, fall into the attached category'. .Accordingly, it is 
concluded that the Cytherean bow wave must nearly always 
be detached. However, dunng the solar minimum epoch of 
the Venera orbiter and the Mariner 10 flyby, the shock can 
indeed lie quite close to the planet. 

The low-altitude flow field in the Venus lonosheath is 
examined further in Figures 12 and 13. The top panels of 
each figure display the exact streamline patterns from the 
best fit gas dynamic models in Figures 9 and 10. In both 
cases, although more so for the Venera model, a ponion of 
the gas dynamic streamlines intersect the observed lono- 
pause boundary. If there were no neutral atmosphere ex- 
tending up into the lonosheath. then one might conclude that 
the particles on the streamlines crossing the ionopause were 
simply being lost to the ionosphere/atmosphere in a lunar- 
type interaction. However, such a conclusion would contra- 
dict the theoretical calculations indicating that the Cytherean 
ionosphere is a sufficiently good electrical conductor to 
exclude almost all of the solar wind [e.g.. Cloutier and 
Darnell, 1973; Zeleny and Vaisberg. 1981]. In fact, the 
neutral atmosphere does permeate the lonosheath [Cravens 
et al., 1980; Nagy et al.. 1981]. These neutrals will interact 
with the solar wind plasma through charge exchange and 
photoion pickup to become mass/momentum/energy terms 
in the flow field that are not taken into account by the 
conventional gas dynamic calculations used in this investiga- 
tion. As demonstrated in previous theoretical studies [Bier- 
mann et al., 1967; Cloutier et al., 1969; Wallis. 1973. Hard 
and Siscoe in Bridge et al., 1976], the net result is to divert 
the incoming flow about the volume occupied by the neutral 
population. Accordingly, these interactions at Venus can 
result in a greater deflection of lonosheaih flow at higher 
altitudes than would be expected from the thermal and 
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magnetic pressure gradients of a neutral atom free magneto- 
sheath [e.g., Wallis, 1982], Charge exchange interactions, 
will transfer momentum and energy from the solar wind to 
the hot neutral atoms they create while the pickup of cold 
ions created by both photoionization and charge exchange 
will add mass to the flow field and transfer flow energy to 
thermal energy. Although the detailed high time and spatial 
resolution particle distnbution measurements required to 
model these processes will not be made for some time, it 
may be that one effect of these processes is to introduce 
additional degrees of freedom to the flow system. They could 
therefore lower the value of the adiabatic exponent (i.e., y = 
{n T 2)/n. for a perfect gas) appropriate to the idealized 
equation of state in our gas dynamic model calculations. For 
these reasons it is probable that plasma-neutral atmosphere 
interactions, as opposed to a lunar-type absorption of a 
portion of the solar wind flux, are the main reasons for the 
small width of the Venus ionosheath relative to Eanhiike 
magnetosheaths. The bottom panel of Figures 12 and 13 
present a qualitative flow pattern that is consistent with the 
high-altitude gas dynamic results (i.e., top panels), the need 
for a general exclusion of the solar wind by the ionosphere, 
and the anticipated effects of plasma-neutral interactions 
discussed above. The solar wind plasma on the streamlines 
closest to the stagnation line penetrate the deepest into the 


Gasdynamic Model of Flow About Venus at Solar Minimum 
(Ms=8.y=2, H/Ro=Q20. and 1 Roe'0-91 Ry) 



Ionospheric - Cometary Model Venus Solar Wind Interaction 
Region Region Region 
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Fig. 12. In the top panel the gas dynamic flow field derived from 
the position and shape of the Venus bow wave measured by Venera 
are plotted in relation to their radio occuliation ionopause [Ivanov- 
Kholodny et al., 1979) A qualitative flow model, including the 
expected effects of charge exchange and photoion pickup, is dis- 
played in the bottom panel 


Gasdynamic Model of Flow About Venus at Solar Maximum 
(Ms=8, y = 2. H/R(J=0.20. and lRoe= 101 R^) 



Ionospheric - Cometary Model Venus Solar Wind Interaction 
Region Region Region 



Fig. 13. The same as the preceding figure but for the Pioneer 
Venus modeling results. 


neutral atmosphere external to the ionopause and take the 
greatest amount of time to flow about the obstacle. These 
streamlines interact most strongly with the neutrals and are 
designated region 1 in the bottom panels. To what extent and 
in what manner the ionosheath plasma (i.e., solar wind + 
ionized neutrals -i- ionosphenc plasma!?)) in this region 
flows downstream can only be determined by further obser- 
vation and modeling of ionosheath phenomena (e.g.. the 
recently discovered 'plasma clouds,’ Brace et al. [1982]; 
Bussell et al. [1982]). Region 2 is filled with the pomon of the 
flow that interacts more weakly with the neutral atmosphere 
than region 1 because of its overall greater altitude. These 
streamlines have been drawn to follow the contour of the 
observed ionopause surface. Those streamlines residing at 
sufficiently great heights to be essentially unaffected, like the 
upstream flow, make up region 3. They correspond exactly 
to the gas dynamic results of the upper panels. The plasma 
depletion, or magnetic bamer. region [Elphic et al., 1980] 
has somewhat arbitranly been drawn as being a constant 100 
km thick [Zwan and Wolf, 1976]. This region has not yet 
been well mapped observationally or intensively modeled 
with theory. Its actual extent and degree of depletion will 
vary in three dimensions. For these reasons, no attempt has 
been made to extrapolate the two-dimensional streamlines of 
Figures 12 and 13 into the depletion layer. 

Figure 14 displays these same results, but from a different 
perspective, by showing the upstream mapping of the low- 
allitude flow field in the preceding two figures. The Venera 9 



Slavin et al.; Solar Wind Flow About the Planets 



Fig. 14 The results of the preceding two figures are displayed 
from a different perspective, which shows the fractions of the 
incident solar wind interacting directly wuh the planet and the global 
two-dimensional streamline pattern. 

and 10 solar minimum flow field places 10% of the incident 
streamlines in region 1 and 10% in region 2. At solar 
maximum the PVO mode! shows a weaker interaction, with 
the neutrals producing a region I and region 2 with 1% and 
9%, respectively, of the incoming streamlines. Hence, this 
study concludes that 10-20% of the solar wind streamlines 
incident over the dayside of Venus interact strongly with the 
planet. This result places an absolute upper limit on the 
mean electric potential drop across the ionosphere of ap- 
proximately 0.15 x 60 kV = 9 kV pole to pole. However, the 
shielding effect of the high electrical conductivity of the 
ionosphere may reduce the actual drop to order 10 V or less 
[Cloutier and Daniel!, 1973; Zeleny and Vaisberg, 1981]. 
Downstream, the region 1 and 2 streamlines will still be 
appreciably affected by the accretion of new ions and will 
make up the boundary layer region that has been previously 
observed (e.g., the Penumbra region, Gringauz et al. [1976], 
Vaisberg et al. [1976n]; mantle region, Spenner et al. [1980]). 

As regards the factor of 2 variation in the fraction of the 
incident solar wind contained in regions 1 and 2, we propose 
that this phenomenon occurs directly as a result of the solar 
cycle modulation of tonopause height first discovered by 
Wolff et al. [1979], In their model, lonopause altitude varies 
with the phase of the solar cycle. It is lowest at solar 
minimum because of the solar wind dynamic pressure peak 
and solar EUV emission trough that occur during the low 
portion of the cycle. Hence, near solar minimum the number 


of atmospheric neutrals exposed to the solar wind will be 
maximized in relation to the rest of the cycle. Theoretical 
support for a factor of 2 increase in the strength of the 
comelary component of the Venus interaction caused by the 
relatively small change in lonopause altitude in Figures 12 
and 13 is found in the calculations by Gombosi et al. [I98IJ. 
In their study a gas dynamic flow field was assumed, and the 
total number of charge exchange reactions in the lonosheath 
was computed by using Pioneer Venus neutral atmosphere 
models. They found that 1-16% of the incident protons 
underwent charge exchange as a function of the assumed 
height of the ionopause. Hodges and Tinsley [1981] have 
recently produced the same order of magnitude results as 
Gombosi et al. [1981] in a modeling study of the Venus 
exosphere. In the Gombosi et al. [1981] work a change in 
ionopause height of the amount shown in the figures pro- 
duces a variation of only about 30% in the total number of 
charge exchange interactions. By comparison, we found a 
100% growth in the solar wind flux in regions 1 and 2. 
However, since their model does not take into account 
photoion pickup, the modification of the flow field with mass 
addition, or the actual shape of the ionopause. the factor of 3 
difference between our model of the obser\'ations and their 
theoretical calculations does not appear to be serious. In 
addition, McElroy et al. [1982] have estimated that iono- 
sheath electron impact ionization will be a comparable or 
greater source of oxygen ions than photoionization. Hope- 
fully, the Pioneer Venus orbiter will continue to be funded 
until its reentry in the early 1990’s. In this event it should be 
possible to reach some definite conclusion about the exis- 
tence and causes of the solar cycle lonosheath variations we 
have proposed. 

On the basis of our modeling results it appears necessary 
to briefly review and evaluate the current state of theory and 
observation concerning the solar wind interaction with Ve- 
nus. Pnor to the Pioneer Venus mission, three factors were 
stressed as being of possibly great imponance to the physics 
of the interaction: cometary processes involving the neutral 
atmosphere [Wallis, 1972, 1973; Bridge et al., 1976], iono- 
spheric currents induced by the solar wind [Cloutier and 
Darnell, 1973; Cloutier et al., 1982], and the-effects of any 
small intrinsic magnetic fields [Russell, 1976a, b]. Low- 
altitude nightside PVO measurements [Russell et al., 1980] 
have now confirmed the Venera 4 results of Dolginov et al. 
[1969], indicating no significant planetary magnetic field. 
The suong plasma-neutral atmosphere interaction with a 
large reduction in solar wind Mach number upstream of the 
bow wave, proposed by Wallis, has also failed to matenalize 
[Mihalov et al., 1980, 1982; Slavin et al., I980J. However, 
the more modest comelary interaction within the iono- 
sheath, proposed by this study, appears consistent with the 
theoretical models of Hard and Siscoe (reported in Bridge et 
al. [1976]) as well as with the recent Gombosi et al. [1980, 
1981] and L. M. Zeleny and O. M. Vaisberg (unpublished 
manuscript, 1981) calculations. Observationally, the pres- 
ence of a suprathermal ion layer just exterior to the iono- 
pause [Taylor et at., 1980] and discovery of O' ions being 
carried off by the solar wind in the downstream lonosheath 
[Mihalov et al., 1980; Mihalov and Barnes. 1981; Intriliga- 
tor, 1982] along with clouds of apparently ionospheric plas- 
ma [Brace et al., 1982; Russell et al.. 1982] indicate the 
presence of strong interactions in the low-altitude iono- 
sheath. However, the exact role of the neutral atmosphere in 
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these processes has not yet been determined. It has not yet 
been demonstrated that sufficient mass is being picked up to 
account for the observed slowing of the solar wind near the 
flanks of the ionopause (Mihalov et al. [1982]; H. Perez-de- 
Tejada, unpublished manuscnpt, 1982). Finally, while many 
aspects of the interaction remain unresolved, the incident 
solar wind is indeed diverted about the planet by large 
induction currents within the ionosphere [Elphic et al., 1980] 
as predicted by the Cloutier and Darnell models. Pioneer 
Venus has gone on to present the basic unipolar induction 
concept with a number of new phenomena to be incorporat- 
ed, such- as ionospheric flux ropes [Elphic et al., 1980], a 
Zwan-Wolf plasma depletion layer exterior to the ionopause 
[Zwan and Wolf, 1976; Elphic et at., 1980; Vaisberg et al., 
1980], occasional intervals of strong ionospheric magnetic 
fields [Elphic et al., 1980; Luhmann et al., 1980], and a 
boundary layer along the flanks of the ionopause [Perez-de- 
Tejada, 1980; Spenner et al., 1980; Taylor et al., 1980]. 
Definitive observations regarding many of these questions 
will probably require a large number of high time resolution 
measurements from a low-altitude orbit skimming the iono- 
sheath-ionosphere interface such as that being considered 
for VOIR. 

Mars 

The solar wind interaction with Mars has been the subject 
of controversy for over 15 years. As considered in the 
reviews by Michel [1971], Hill and Michel [1975], Bauer 
[1976], Ness [1979], Russell [\919], Siscoe and Slavin [1979], 
and Gnngauz [1981], the very limited panicles and fields 
measurements gathered by the Manner 4 distant flyby and 
the Mars 2,3, and 5 orbiters have not provided a definitive 
picture of the solar wind interaction with this planet. While 
the bow shock has been sampled on numerous occasions, 
none of these satellites probed the region below 1000 km 
where the solar wind interaction must take place. Funher, 
the total number of observations were small, and the instru- 
ments were often turned off during the periapsis portions of 
their orbits [Gringauz, 1976]. Given the nature of the existing 
data, there has naturally developed a question over whether 
Mars is magnetized at all or stands off the solar wind through 
a Venus-type solar wind interaction. 

For this reason it is important that the average flow 
pattern about Mars that is displayed in Figure 11 does not 
appear at ail like that of Venus in Figures 9 and 10. Given the 
same adiabatic exponents, the Mars gas dynamic obstacle 
lies well above the ionosphere proper while the Venusian 
obstacle resides at or below the altitude of the ionosphere. 
Figure 8 funher stresses this difference between Mars and 
Venus by comparing the nose positions of the Martian bow 
shock along with the Venera and PVO observations at Venus 
from paper 1. While all three distributions overlap to some 
degree, the Mars bow wave is indeed the most distant overall 
and is far removed from the attached regime. Slavin and 
Holier [1982] have conducted a study of the Mars solar 
wind-obstacle interface in Figure 1 1 and the ability of the 
ionosphere to support it. Using PVO Venus interaction 
models, they showed that the Martian ionosphere could not 
stand off the average solar wind at these heights, even if it 
were always magnetized through the capture of IMF field 
lines, as happens only sometimes at Venus. Their results 
supported the genertil findings reached by some of the earlier 
studies, which did not consider induced magnetic fields 



Fig. 15. The inbound and outbound January 21, 1972, xMars 3 
bow shock and magnetopause crossings have been modeled sepa- 
rately with the same gas dynamic solution found to well represent 
the average Mars observations in Figure 12, 

[Gnngauz, 1976; Intriligator and Smith. 1979], Slavin and 
Holzer further concluded that if Mars were in fact lacking a 
significant intrinsic magnetic field, then the solar wind would 
often push the ionosphere to altitudes below 300 km, where 
plasma-neutral atmosphere interactions would absorb and 
deflect the flow [Cloutier et al., 1969], Such an interaction 
with the tenuous atmosphere of Mars would produce a bow 
wave location lower than that of Venus, Figure 8 shows this 
not to be the case. Further, no such low-altitude ‘iono- 
pauses’ have been observed m the extensive Viking radio 
occultation measurements [Lindal et al.. 1979J. 

Using the scale heights measured by the Viking landers 
[Hanson et al., 1977], Slavin and Holzer went on to show 
that lonosphenc plasma pressure will make a negligible 
contribution to standing off the solar wind at the altitudes of 
the gas dynamic Mars obstacle in Figure 11. With this 
knowledge the effective dipole magnetic moment of Mars 
may be computed by using the classical relationship between 
magnetic moment and stand-off distance for a dipole in the 
solar wind under the assumption of lerrestnal-type magneto- 
spheric current systems [e.g., Spreiter et al., 1966; Holier 
and Slavin, 1978] 

M = {2-rr (2) 

where klf- is a constant equal to 0.59, and 7?, is the 
magnetopause stand-off distance, 3900 km. The result ob- 
tained by Slavin and Holzer was 1.4 = 0.6 x 10^- G cm^, 
which is sufficient to produce a small magnetosphere at 
Mars. 

Having briefly reviewed the analysis of Slavin and Holzer, 
this investigation makes further application of gas dynamic 
modeling to the Mars interaction problem by examining 
some individual events from the published Mars orbiter 
records. In particular, the January 21, 1972. observations of 
Mars 3 are considered in Figure 15. Those measurements 
have been interpreted by the magnetometer [Dolginov et al., 
1973] and Faraday cup [Gnngauz et al., 1974] expenmenters 
as depicting a brief encounter with an intrinsic field magneto- 
sphere between approximately 21;30 and 22:10 (Moscow 
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Fig. 16 Both the study by Bogdanov and Vaisberg (19751 and 
paper 1 noted some unusually distant bow shock crossings by Mars 
; and 3 on December 14-15, 1971. These crossings are modeled with 
a low-Mach-number gas dynamic calculation by using an obstacle 
shape close to that determined in Figure 11. 


time). Wa///s [1975] and 7?wjse// (1978a, b] have disputed this 
analysis on numerous grounds, discussed in Slavm and 
Holier [1982] and many of the review papers cited earlier. 
As shown, we have modeled the Dolginov/Gringauz magne- 
topause crossings and the Vaisberg and Smirnov shock 
encounters by using the average gas dynamic model from 
Figure 11. The agreement is excellent between the model 
and the observations when the inbound and outbound passes 
are considered separately, but good results are also apparent 
when only average magnetopause/bow shock positions are 
used. The first pair of boundary encounters are well matched 
with a subsolar obstacle radius of 1.78 ^ms J ^ms = 
3390 km) while a smaller /?ob = l '^6 /?ms required for the 
outbound encounters. Such a situation might occur if, for 
example, a pressure pulse in the solar wind (e.g.. an inter- 
planetary shock wave) arrived while Mars 3 was within the 
magnetosphere. The height of the magetopause would then 
be depressed until the spacecraft again found itself in the 
magnetosheath. This scenario would also explain the lower 
altitude of the exit from the magnetosphere as compared 
with the entry identified by Dolginov/Gringauz. 

Solar wind ion speeds were obtained from the plasma 
analyzer measurements, but no reliable densities and tem- 
peratures [Dolginov, 1978a]. Bogdanov [1975; see Dolginov, 
1978a] has assembled solar wind dynamic pressure values 
for many of the Mars 2 and 3 passes near the planet by using 
the in situ ion speed measurements and densities extrapolat- 
ed from IMP 6 at 1 AU. For the January 21, 1972, events the 
Mars 2 and IMP 6 parameters are = 500 km/s and n = 


1.2 cm"^ to yield P,w = 0.6 x lO”* dynes/cm^ While solar 
wind dynamic pressures determined in this manner are 
subject to significant errors, it is desirable to examine the 
planetary magnetic moments that may be derived from the 
obstacle sizes in Figure 15 by using (2). The effective dipole 
moments arrived at are 3.3 x 10^*, 1.8 x 10^. and 2.5 x 10" 
G cm^ for the inbound, outbound, and average stand-off 
distances. These values are similar to. but slightly greater 
than, the 1.8 x 10^ G cm^ Dolginov [1978a. b] determined 
by inversion of the magnetic field vectors near penapsis and 
the 1.4 ± 0.6 X 10^- G cm^ obtained from gas dynamic 
modeling of bow shock shape/position by Slavin and Holier 
[1982]. Hence, given the uncertainties in the solar wind 
dynamic pressure on January 21, 1972. our gas dynamic 
analysis of these Mars 3 observations suppon the conclusion 
that Mars possesses a significant intrinsic magnetic field and 
does not undergo a Venus-type ionosphenc interaction with 
the solar wind. 

In paper 1 a pair of shock crossings by Mars 2 and 3, 
separated by only about 7 hours on December 14-15, 1971. 
were seen to lie far outside of the average bow wave 
position. They were noted earlier by Bogdanov and Vais- 
berg [1975]. who calculated a subsolar obstacle radius of 1.8 
/?MS for the event. In the absence of upstream parameters. 
Vaisberg [1976] suggested that these crossings might be 
more the result of low solar wind Mach numbers producing a 
thicker-than-average magnetosheath than low dynamic pres- 



Fig. 17. The Viking 1 lander bow shock crossing identified by 
Cragin et al. [1982] is plotted in aberrated coordinates and compared 
with the Mars 2. 3, and 5 bow wave encounters considered in paper 
I. 
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sure causing an expansion of an intrinsic field magneto- 
sphere. We examine this question in Figure 16 with a low- 
Mach-number gas dynamic model based on an obstacle 
shape similar to those utilized in Figures 11 and 15. The 
A/5 = 4 flow solution yields a moderately good fit to the 
crossings but with a diflference in flare angle that suggests a 
sightly lower Mach number as being optimum. This point 
was not pursued further because lack of funds prevented 
production of lower-Mach-number models. However, the 
important result is that the assumption of an unusually low 
Mach number still does not negate the need for a solar wind 
stand-off distance far in excess of the height at which an 
lonopause might form. Hence, not only low-Mach-number 
conditions but also a large-diameter obstacle, such as an 
expanded intrinsic field magnetosphere, appear necessary to 
explain the December 14-15, 1971, Mars orbiter bow wave 
observations. 

Recently, the high-altitude (i.e., 16,000 to 350 km) Viking 
lander retarding-potential analyzer measurements have been 
examined by Cragin et al. [1982]. As detailed in their paper, 
a distinct peak in the rms deviations of Viking 1 lander RPA 
electron fluxes was observed near the altitude of the mean 
Mars bow shock determined by the Mars 2, 3, and 5 orbiiers. 
They suggest that this signature is associated with wave 
particle interactions in the Martian foreshock. In Figure 17 
we have aberrated and plotted the Cragin-et-aJ.-suggested 
shock encounter at 1700-km altitude; this is relative to the 
Soviet orbiter crossings modeled in paper 1. The agreement 
IS good with the Viking 1 lander datum, only somewhat 
lower in altitude than the average as seen in Figure 8. 
Further, there is no evidence for the larger solar maximum/ 
minimum variation very evident in the Venus shock observa- 
tions. The Viking 2 lander followed a different trajectory, 
which did not take it outside of the magnetosheath. 

Figure 18 reproduces the Viking 1 electron current rms 
deviation and mean value as a function of altitude from 
Cragin et al. [1982] for comparison with our gas dynamic 
Mars flow model in Figure 1 1. As shown by cross hatching at 
the top and bottom of the figure, the gas dynamic model at a 
solar zenith angle of 40° places the obstacle altitude at 490 
km if the bow shock is assumed to lie near the sudden jump 
in a near 1850 km. Given the very limited nature of the 
measurements, it' may be argued that the shock could have 
lain somewhere else between 3000 and 700 km. However, it 
is interesting to note that the theoretical obstacle in Figure 18 



Fig. 18. The gas dynamic flow model describing the average 
Mars observauons in Figure 12 is used 10 model the Viking I lander 
RPA mean negauve current and coefficient of variation measure- 
ments of Cragin el al. [1982]. 



Psw (10"® dynes/cm^ ) 

Fig. 19. Seven Mars 3 and 5 shock crossings for which reliable 
solar wind dynamic pressure (i e.. 1.16 n^mpY-) values exist have 
been graphed against Aju. 

falls 200-300 km below the rapid drop in electron flux near 
700-800 km, which Cragin et aJ. suggest might indicate entry 
into a plasma depletion layer similar to that at Venus. Such a 
thickness is similar to that predicted for Mars by Zwan and 
Wolf [1976] and observed by Elphic et al. [1980] at Venus. 
In agreement with the results of Intnligator and Smith [1979] 
and Slavin and Holier [1982], Cragin et al. found that the 
ionosphere of Mars could not have stood off the solar wind 
on this occasion for any reasonable set of ionosphenc 
plasma conditions. 

In Figure 19 the dependence of the height of the Mars bow 
shock upon solar wind dynamic pressure is investigated. For 
seven of the Mars shock crossings in paper 1 there exists 
either in situ dynamic pressure measurements [i.e.. Mars 5; 
Dolginov, 1976] or good agreement between the 1 and 1.5 
AU solar wind speed observations during the interval for 
which density was extrapolated to Mars orbit [Dolginov, 
1978a]. The average shape of the shock from paper 1 was 
then used to map each of these four Mars 5 and three Mars 3 
shock encounters down to the nose position with the results 
plotted in Figure 19. The formal errors associated with these 
mappings caused by the uncertainties in the mean surface 
(see Table 7 of paper l)"and the locations of the 'crossing 
average about 0.05 Rms- For the sake of comparison with an 
ionospheric interaction we have also displayed the distance 
to the nose of the Venus bow wave as a function of upstream 
dynamic pressure determined with the Pioneer Venus obser- 
vations by Slavin et al. [1980]. The exponential variation of 
ionospheric plasma pressure with altitude makes Venus a 
relatively incompressible obstacle, as can be seen in the 
small slope of the dashed line in the figure. The much 
stronger pressure dependence of a simple magnetic dipole of 
magnitude 1.4 x 10“ G cm^ is plotted with a long-short 
dashed line. While the few Mars data points do not contra- 
dict this simple dipole relationship at low there is only a 
slight additional decrease in shock distance as pressure 
increases above 2 x 10"® dynes/cm’. There are also two 
crossings at relatively low pressures. To emphasize the 
potentially strong role of magnetic reconnection in determin- 
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mg the configuration of a planetary magnetosphere with a 
high ionosphenc electrical conductivity [Rassbach et al.. 
1974], we have speculatively labeled these two points as 
erosion events. Best fits to all seven of the points, with the 
lowest excluded, are indicated. In both cases the pressure 
dependence is weaker than that of a simple dipole but 
stronger than at Venus. Although the presently available 
observations do not warrant additional modeling, we note 
for future investigations that ionosphenc induction currents 
can be expected to stiffen the Martian dipole field lines to 
compression in a fashion analogous to that predicted for the 
small magnetosphere of Mercury [Hood and Schubert, 1979; 
Suess and Goldstein, 1979]. The results in Figure 19, sug- 
gesting a Venus-type dependence at higher dynamic pres- 
sures and an Earthlike response at lower pressures, may be 
quite reasonable for a small Martian magnetosphere. 

Conclusions 

1. The finding by Fairfield [1971] that the average posi- 
tion of the dayside bow shock/magnetopause boundaries is 
well predicted by single-fluid gas dynamic theory has been 
confirmed and extended. The new gas dynamic code of 
Spretier and Stahara [1980a] is found to improve the agree- 
ment between theory and experiment by allowing for the 
specification of an obstacle shape based upon magnetopause 
observations as opposed to using less realistic theoretical 
models. The average locations and shapes of the forward 
bow shock determined by four different satellite missions 
have been found to be well ordered by solar wind dynamic 
pressure and sonic Mach number. A single-fluid adiabatic 
exponent of 2 produced better gas dynamic fits to the 
observed width of the magnetosheath than a value of 5/3. On 
the basis of these results it is concluded that the gas dynamic 
approximation to the fluid descnption of solar wind impinge- 
ment on the terrestrial magnetosphere is generally valid for 
average upstream conditions. 

2. Application of gas dynamic modeling to Venus has 
shown that the ionosheath is significantly narrower than 
would be predicted for a simple terrestnal-type magneto- 
sheath of comparable shape. Further, the ionosheath is 
thinner at solar minimum in the Venera observations than at 
solar maximum as measured by Pioneer Venus. The most 
probable causes of both the overall thinness of the lono- 
sheath and its apparent solar cycle modulation are the 
ionized neutral pickup and solar wind-exosphere charge 

- exchange processes taking place in the io no sheat h. Th e 
lower ionopause [Wolff et al., 1979] during cycle minimum 
will enhance the number of atmospheric neutrals exposed to 
the solar wind by allowing the incoming flow to penetrate 
closer to the planet than is the case at solar maximum. From 
the Venera and PVO gas dynamic flow models it is estimated 
that 20% and 10%, respectively, of the incident solar wind 
streamlines are significantly modified by plasma-neutral 
particle interactions. 

3. Gas dynamic modeling of the Mars interaction, unlike 
the situation at Venus, places the solar wind-obstacle inter- 
face well above the heights at which this planet could stand 
off the solar wind in the absence of an intrinsic magnetic 
field. It was found that the bow shock crossings observed by 
Mars 3 on January 21, 1972, were indeed consistent with this 
spacecraft having encountered a modest intrinsic field Mar- 
tian magnetosphere near penapsis, as maintained by Dol- 
gtnov et al. [1973] and Gnngauz et al. [1974], Because its 


magnetic moment is small, the ionosphere and neutral atmo- 
sphere may be expected to sometimes panicipate directly in 
the solar wind interaction through induction and. more 
rarely, cometary processes when the solar wind dynamic 
pressure is high or the dayside magnetosphere is significant- 
ly eroded by reconnection. Possible evidence for such 
effects was found in the observed dependence of bow wave 
altitude upon upstream dynamic pressure. 
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ABSTRACT 


This study uses observations by a number of spacecraft to investigate 
the asjrmptotic behavior of planetary bow shocks. Toward this end a 
single standard method has been used to model distant bow shock position 
and shape. Mach cone angles of 13.9 +. 2°, 11.4 + 3°, and 8.1 ^ 4^ at 
Venus, Earth, and Mars, respectively, were determined from the 
observational shock models. These cone angles and their decrease with 
growing distance from the sun are generally consistent with downstream 
bow shock position being limited by the MHD fast mode Mach number. 
Gasdynamlc solutions for solar wind flow about Venus, Earth, and Mars 
Were computed up to 50 Rqq (i.e. obstacle radii) behind the planets and 
compared with observed bow shock location. In each case the position of 
the shock was well predicted up to a certain distance downstream: -4 
at Venus, -6 Rq^ at Earth, and -10 RQg at Mars. Beyond this point the 
observed shock position lies farther from the aberrated sun-planet line 
than the gasdynamlc model with the discrepancy greatest-at Venus and 
least at Mars. The better agreement between gasdynamlc theory and 
observation with growing distance from the sun is attributed to an 
increase in the accurancy of the gasdynamlc approximation with 
decreasing IMF strength. These results suggest that gasdynamlc theory 
will be most accurate in describing flow about large bodies in the outer 
solar system. 



Introduction 


Previous studies (Fairfield, 1971; Slavin et al., 1983) have shown 
that the position and shape of planetary bow shocks forward of about one 
obstacle radius behind the planet are well predicted by single fluid 
gasdynamic theory (Spreiter et al., 1966; Sprelter and Stahara, 1980). 
In these theoretical models the presence of the interplanetary magnetic 
field is n^lected on the grounds of its smallness in the sense 

/(BV&ir) <^100. The flow field is then completely determined by 
upstream sonic Mach number. Mg, adiabatic exponent, 7 *, and the obstacle 
shape/position. 

While this approach is mathematically valid as an approximation to 
the numerically Intractable MHD equations, the requirement that the 
ambient magnetic field be neligible is not met everywhere. For example, 
investigations of the regions just exterior to the Venus ionopause 
(Elphic et al., 1980), Earth magnetopause (Crooker et al , 1979), and 

Saturn magnetopause (Smith et al., 1980) have discovered a thin layer of 
- compressed magnetic field adjacent to these obstacles from which the 
solar wind plasma has been lost via field aligned flows. These layers 
are not present in gasdynamic theory, but were predicted by one 
dimensional (Lees, 1964) and quasi-two dimensional (Zwan and Wolf, 1976) 
MHD studies of the stagnation region. 

Another place where the gasdynamic approximation may be a poor 
representation is near the distant downstream bow shock. In this region 
the characteristic lines start to intersect the shock at small angles, 
as shown in Figure 1, indicating that the shock is approaching its "Mach 
cone" limit (e.g. Landau and Lifshitz, 1959). Gasdynamic models may not 
be able to accurately predict shock location in this region because the 
gasdynamic approximation, M^IO, becomes invalid as the shock weakens. 
Mach number in the shock frame decreases as less deflection of the flow 
is required and the wave begins to propagate more across than against 
the solar wind flow direction. The bow shock continues to weaken until 
the Mach number in the shock frame goes to unity and it decays into a 
compressional wave. Hence, at some intermediate distance downstream the 
gasdynamic approximation becomes invalid and the agreement between 
theory and observation may deteriorate. 

This study investigates the asymptotic behavior of planetary bow 
shocks and the ability of gasdynamic theory to describe it. Spacecraft 
observations at Venus, Earth, and Mars are used to model the shapes and 
positions of their distant bow waves utilizing a standard method. The 
measured planetary Mach cone angles are compared with the mean sonic and 
MHD fast wave Mach numbers at 0.7, 1.0, and 1.5 AU to assess the 
downstream influence of the interplanetary magnetic field. Finally, 
numerical gasdynamic flow solutions extending up to 50 obstacle radii 
behind the planets are obtained and tested against the observed location 
of the downstream bow shock. 

Modeling Bow Shock Position 

The bow shocks of Venus and Mars have been modeled using the three 
parameter second order method of Slavin et al . (1980; Slavin and Holzer, 
1981). The terrestrial shock was not modeled because the techniques 
applied by Fairfield (1971) have already been shown to be quite similar 
to those employed here (Slavin and Holzer, 1981). The fitting 



parameters are the surface eccentricity, £. , the semi-la tus rectum, L, 
and the location of the focus along the x' axis, x^. The shock itself 
is assumed to possess cylindrical symmetry about the solar wind 
direction which is taken to_^be qoposite the aberrated planet centered 
solar ecliptic x axis (i.e. -x'). All modeling is performed 

in the aberrated coordinates using either- in situ solar wind speed or a 
mean value of 430 km/sec. 

The near planet bow shock has been found previously to be highly 
symmetric about the x' direction in agreement with the predictions of 
gasdynamic theory (Slavin et al., 1980; Slavin and Holzer, 1981; 
Tatrallyay et al., 1983). Farther downstream the asymmetric nature of 

the MHD fast wave Mach cone may become significant. However, the effect 
should be small and there are too few distant shock obervations to model 
this region as a function of the upstream IMF orientation. 
Uncertainties introduced by our assumption of cylindrical symmetry for 
the downstream shock will be discussed in a later section. 

The equation for the model surface is 

r = L/(i +£cos6) (1) 

where r is radial distance from the focus (i.e. x'=x^) and 8 is measured 
about the focus from the positive x' axis direction. When the 
eccentricity in (1) is greater than unity, as is the situation in 
studies of the distant shock, the resultant curve is a hyperbola. For 
any hyperbola asymptotes exist which bound its position and correspond 
to the Mach cone for shock waves. They depend only upon the 
eccentricity and yield an expression for experimentally determined Mach 
cone angle 

C<Qgg= tan-l((£.2-l)l/2) (2) 

This expression can be differentiated to obtain the uncertainty in Mach 
cone angle , as a function of the error in the model eccentricity 

(«/£)/ ( 3 ) 

Figure 2 displays typical orbits for Venera 9/10 and the Pioneer 
Venus Orbiter, PVO, in relationship to an extrapolated near planet model 
surface (Slavin et al., 1980). While PVO provides some excellent 
cbiTerage of the forward shoch when periapsis is on the nightside, it is 
a very poor source of information on the location of the distant bow 
wave because of its tendency to parallel the shock (see discussion of 
these problems in Slavin and Holzer, 1981; Tatrallyay et al., 1983). 
Small variations in the flare angle of the shock appear in the PVO 
observations as boundary crossings at larger or smaller values of the x' 
coordinate. For this reason we have not used the PVO observations in 
modeling the downstream Venus shock. The lower inclination and less 
eccentric nature of their orbits make the Venera 9 and 10 shock 
crossings (Smirnov et al., 1981) better suited for modeling the distant 
Venus shock. Near planet shock crossings by Pioneer Venus were also not 
included because of the possibly solar cycle related variation in the 
altitude of the forward Cytherean shock discovered by Slavin et al . 
(1979; 1980). With the exception of the Mariner 5 and Venera 4/6 
downstream crossings, the Venus bow shock model produced by this study 
is based on solar minimum measurements. 



In l/r,cos(6) coordinates centered on the geometric focus second 
order curves become straight lines (see equation 1). Linear least 
square fitting techniques .may thus be applied. Figure 3 displays 
Mariner 5 (Bridge et al., 1967; Russell, 1977), Venera 4 and 6 (Gringauz 
et al., 1970), Mariner 10 (Bridge et al., 1974; Russell, 1977), and 
Venera 9 and 10 (Smirnov et al., 1980) bow shock crossings in these 
"conic" coordinates along with a least square best fit. The focus 
location on the x' axis was then systematically varied until the 
root-mean-sqxiare deviation of the crossings measured normal to the best 
fit surface was minimized. Figure 4 plots the Venus crossings in the 
conic coordinates again, but with the focus now centered on the optimum 
Xjj*»+0.45 R^. The rms deviation in x,(y'^+z'^)^' ^ space is decreased by 
a factor of 2.6 in going from to x =0.45 R^. Figure 5 displays the 
shock crossings and best fit model surface in the more usual aberrated 
cylindrical coordinates. The uncertainties in the model eccentricity 
and seml-latus rectum are approximately tlX. This is somewhat less than 
was found for the near planet Venera 9/10 and Pioneer Venus models by 
Slavin and Holzer (1981) because the inclusion the downstream 
observations permits less variation in the fitting parameters. The Mach 
cone determined from our shock model and equations (2) and (3) is 
13.9°^2°. This result is intermediate between the 10.8° and 15.4° cone 
angles implied by the PVO distant shock model of Tatrallyay et al. 
(1983) and Venera 9,10 model by Smirnov et al. (1980). The next section 
will compare this result to the sonic and MHD Mach cones expected on the 
basis of average solar wind parameters. 

In passing it is noted that the more extensive investigation of the 
distant Pioneer Venus distant shock crossings by Tatrallyay supports the 
preliminary finding of Slavin et al. (1980) that no decrease in solar 
wind Mach number took place bewteen the Venera 9,10 and PVO intervals of 
observation. If anything, Mach number appears to have been higher 
during the 1979-81 period. A large decrease in solar wind Mach number 
would have provided an alternate explanation for the change in forward 
bow shock position betweem 1975-6 and 1978-79 attributed to the effect 
of solar cycle changes in solar EUV and solar wind dynamic pressure by 
Slavin et al. (1980; Slavin and Holzer, 1983). 

Observations of the terrestrial bow shock downstream of x'=-60 Rg 
have been made by only a few missions: Explorer 33 (Howe and Binsack, 
1972; Mihalov, 1974), Pioneer 7 (Villante, 1976), and Pioneer 8 
(Bavassano et al., 1971). The shock crossings in the better sampled 
x'>-60Rg region have been modeled by Fairfield (1971) with a second 
order method comparable to the procedure employed here. As discussed in 
Slavin and Holzer (1981), the Fairfield best fit to the aberrated shock 
data was a hyperbola with £ = 1.02, L = 22.3 Rg, and x = +3.4 Rg. The 
corresponding Mach cone angle is 11.4°+3° (assuming i 1* fitting errors) 
or 2.5° less than observed for Venus. While the Mach cones inferred 
from the more distant Pioneer and Explorer observations will be 
discussed in the next two sections, we have adopted the Fairfield 
aberrated shock surface for all comparisons with our models of the Venus 
and Mars bow shocks. 

The Martian bow shock has been observed by the Mars 2, Mars 3 and 
Mars 5 orblters and the Mariner 4 fly-by (see reviews by Russell, 1979; 
Gringauz, 1980; Slavin and Holzer, 1982). Sample trajectories for the 
three Soviet orbiters are displayed in Figure 6 along with the locations 
at which they crossed the Martian bow shock (O.L. Vaisberg and V.N. 



Smirnov, private communication, 1979; see Slavin and Holzer, 1981). 
While the Mars 2 and 3 orbits parallel the downstream shock surface, the 
more perpendicular crossings of Mars 5 at about x'= -5 Rqj allow us to 
model the distant shock. Mariner 4 provides an additional shock 
encounter at approximately twice that distance (Smith, 1969). 

Figure 7 displays our best fit to the Mars shock observations. The 
shape and location of the model surface is well determined with ■ 
1.01 (^IZ), L“ 1.68 Rj^ (:tlZ), and x^«^.7Rj^. The Mach cone angle is 
8.1°+4° which is 3.3° less than the Earth vSue and 5.8° degrees below 
the Venus result. No other Martian shock model using observations from 
all four of these missions has been published. However, comparison with 
the Bogdanov and Vaisberg (1975) Mars 2/3 and Mariner 4 model 
demonstrates the importance of the Mars 5 crossings. In the absence of 
the intermediate distance Mars 5 shock data, their model eccentricity 
was approximately 7Z larger and far more uncertain than determined here. 

Mach Cone Angles 


Table 1 compares the Mach cone angles at Venus, Earth, and Mars 
determined in the preceeding section with conditions in the 
interplanetary medium. In particular, the radial scalings for solar 
wind sonic and Alfvenic Mach number discussed in Slavin and Holzer 
(1981) have been adopted and the 0.7, 1.0, and 1.5 AU values listed. 
The limiting Mach cone angle in gasdynamic theory is that associated 
with the sonic wave speed 


<P^g»sln”l(l/Mg) (4) 

Following Table 1, as mean sonic Mach number increases from 6.6 to 7.9 
the cone angle decreases from 8.7° at 0.7 AD to 7.3° near 1.5 AU. 

While no general MHD solutions for flow about a blunt obstacle have 
been found (Spreiter and Rlzzi, 1974), Mach cone type arguments using 
the MHD fast wave speed may be applied to predict asymptotic bow shock 
position in MHD theory (Michel, 1965). Unlike the situation for 
compressional disturbances in an unmagnetized isotropic plasma, fast 
mode wave speed is a function of propagation direction relative to the 
ambient magnetic field. The fast mode propagates faster perpendicular 
to B than parallel to it. Mach number based on the perpendicular fast 
mode speed is 

< 5 ) 

SO that 

<7<j_-sin“'^((M^VMg2)^/2/M^Mg (6) 

Since the perpendicular fast mode speed is significantly larger than the 
sonic speed, the values are 2-3° greater than as shown in Table 1. 
The parallel fast mode propagation speed is the greater of the sonic and 
Alfven speeds. Hence, the parallel MHD Mach number is simply equal to 
Mg due to the average M*>Mg conditions in the solar wind. At oblique 
propogation angles the rast mode speed lies in between the 0° and 90° 
values listed in the table. 



The angle the local shock propagation direction (i.e. the shock 
normal) makes with the magnetic field is function of location on the 
shock surface. Both the mean spiral configuration of the IMF and the 
three dimensional curved nature of the bow shock contribute to the lack 
of symmetry. Over the poles and on the dusk side of the shock surface, 
the tendency is toward quasi-perpendicular geometries with limiting 
angles approaching . Quasi-parallel conditions are more probable on 
the dawn side where the cone angle should be closer to o( n These 

asymmetries in the MHD Mach cone were first quantitatively modeled by 
Michel (1965). It was his conclusion that the variations in MHD Mach 
cone angle with resj>ect to the x' axis are about 1° overall, but 
generally less than 0.5° for the regions near the ecliptic where most of 
the spacecraft observations used in this study were made. For this 
reason it does not appear that any significant additional errors have 
been introduced by our earlier assumption of axial symmetry given the 
2-4° uncertainties in the observational determinations of Mach cone 
angle. In addition, the the larger Mach numbers in the outer solar 
-- system can be expected to produce distant shock shapes that are nearly 
parabolic given the Mars fitting results in this study. From equation 
(3) it is clear that a very large number of observations will be needed 
to accurately measure the smaller Mach cone angles at Jupiter and 
Saturn. 

The observed Mach cone angles in Table 1 decrease with distance 

from the sun as predicted by both gasdynamic and MHD theory^ but the 
variation is larger than expected. The observed Mach cone angle lies 
between and Vy only in the case of Mars although the error bars on 

for the Earth also fall in that range. The measured cone angle at 

Venus lies well above the sonic Mach cone angle , but the lower bar is 
only slightly greater than the perpendicular fast mode prediction. 
Accordingly, the observed locations of the distant bow shocks appear to 
be more consistent with the larger MHD Mach cone than the smaller 
gasdynamic value. The Mach numbers used in calculating the predicted 

cone angles are based upon long term averages of interplanetary 
conditions (Slavin and Holzer, 1981). Hence they may not be the optimum 
values for the various portions of the solar cycle when the Venera and 
Mars orbiters bow shock observations were made. However, no long term 
variation in solar wind Mach number has been found in the solar wind 
records (Slavin et al., 1979) and the Earth observations used by 
Fairfield covered approximately half a solar cycle. 

As mentioned earlier , bow shock encounters downstream of x '“-60 R 
involving Explorer 33, Pioneer 7, and Pioneer 8 have taken place. Mach 
cone angles have been estimated in these cases by modeling of the shock 
surface and examination of the shock normals inferred from jump 
conditions. The distant bow shock model of Howe and Binsack (1972) 

based upon only the Explorer 33 and 35 observations has an asymptote of 
17.7°. This value is greater than all of the Mach angles determined by 
this study and would imply a mean solar wind Mach number of only 3.3. 
However, the modeling method used in that study was never fully 
described and inspection of the fit shows a marked tendency to 
underestimate the altitude of the forward shock The result is an 

enhanced flare angle which is not present in the Fairfield model. Shock 
normals determined for Pioneer 7 and 8 implied Mach cone angles of 20° 

(Villante, 1976) and 22° (Bavassano et al . , 1971), respectively. The 

Mach numbers implied by the Pioneer 7/8 cone angles are 2. 7-2. 9. These 



values are much smaller than the average MHD and sonic Mach numbers 
listed in Table 1. The reason for the large Mach cone angles inferred 
from the shock normal analyses is unclear, but could be associated with 
either unusually low solar wind Mach numbers during the Pioneer 
encounters or waves on the shock boundary. 

Gasdynamic Models 

The Spretier and Stahara (1980; Stahara et al., 1980) gasdynamic 
code has been modified to produce solutions as far as 50 Rq^ downstream. 
Figures 8 , 9, and 10 dlisplay the calculated bow shock locations for the 
Venus, earth, and Mars body shapes determined by the Slavin et al. 
(1983) study and sonic Mach numbers of h and 8 . These two Mach numbers 
were chosen because they correspond to the largest and smallest Mach 
cone angles in Table 1 (i.e. 7.9 from from Venus and 4.2 from 

at Mars). The ability of the gasdynamic solutions to predict forward 
shock position has already been demonstrated by Slavin et al. (1983). 
However, the planetary Mach cones inferred from the shock observations 
in the preceeding section indicate that downstream the bow shocks are 
limited by larger Mach cones than predicted by gasdynamic theory. Below 
we investigate where these departures start and their magnitude. 

The gasdynamic models of flow past Venus are displayed in Figure 8 . 
The near shock is well represented by the Mg =8 theoretical shock forward 
of approximately x'=-4 Farther downstream the bow shock continues 

to flare outward until its slope nearly matches that of the Mg=4 
gasdynamic shock. Figure 9 performs this same comparison with the earth 
shock model of Fairfield (1971). The results are similar, but with the 
disagreement between theory and observation being less severe and 
starting farther behind the planet near x '=-6 Rq^. Finally, Figure 10 
presents the Mars observations. In this case the Mg =8 theoretical model 
does not begin to seriously diverge from the observational model until 
x '=-10 Rq£ is reached, and then only by a modest amount. 

These results are generally consistent with both the expected 
limitations of the gasdynamic approximation and the uncertainties in the 
experimentally determined Mach cones. In all three cases the planetary 
bow shocks approach Mach cones that are greater than predicted by 
gasdynamic theory, but by only a small amount at Mars where^< ^025 <C{j. 
The trend as a function of distance from the sun is also much as would 
be expected from Table 1. Gasdynamic theory is more accurate and may be 
used farther db^str earn at greater distances from the sun where the IMF 
is weaker and the Alfvenic Mach number higher. The implications are 
that gasdynamic theory will be least useful at Mercury, but highly 
accurate for dscribing the flow of solar wind about Jupiter, Saturn, and 
the outer planets. 

Conclusions 


The distant bow shocks of Venus, Earth, and Mars have been modeled 
using a standard method. Mach cone angles determined from the best fit 
shock models are generally consistent with the MHD fast mode speeds 
measured in the solar wind as a function of distance from the sun. The 
finding by earlier studies that gasdynamic theory can accurately predict 
the location and . shape of the forward bow shock has been extended and 
quantified. At Venus the bow shock position and orientation is poorly 



represented by gasdynamic theory much beyond x'=-4 Rqj* The 
disagreement Is smaller at the Earth and begins farther downstream near 
x'=-6 Rog* Finally, at Mars, the gasdynamic bow wave lies quite close to 
the observed shock surface with small discrepancies becoming evident 
only downstream of x'=~10 RQg* The better agreement between gasdynamic 
theory and observation with increasing distance from the sun is 
attributed to the decrease in IMF strength with distance making the 
gasdynamic approximation more accurate. 



TABLE 1. Predicted and Observed Mach Cones 
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Figure Captions 


1. Gasdynamlc models of flow about the earth for an adiabatic exponent 
of 2 and sonic Mach numbers of 4 and 8 are displayed. Obstacle shape in 
the Spreiter and Stahara (1980) numerical codes is parameterized in 
terms a quantity H/R . The coordinates used are aberrated geo^centric 
solar ecliptic (i *e. 1 V*| ) in units of obstacle radii, Rqj* 

2. Sample Venera 9, 10 and Pioneer Venus orbits are shown in Venus solar 
ecliptic coordinates relative to the bow shock. 

3. Mariner 5, Venera 4, Venera 6, Mariner 10, and Venera 9, 10 bow shock 
crossings are displayed In planet-centered "conic" cooridnates along 
with a linear best fit. 

4. The same as the preceedlng figure, but with the conic focus moved to 
x^=+0.45 Ry to improve the fit. 

5. The Venus bow shock crossings and best fit are displayed in aberrated 
planet centered solar ecliptic coordinates. 

6. Sample Mars 2, 3, and 5 orbital trajectories are plotted in relation 
to the observed crossings of the bow shock. 

7. Mariner 4, Mars 2, Mars 3, and Mars 5 bow shock crossings are plotted 

in planet centered solar ecliptic coordinates. The best 3 parameter 

second order least square fit to the boundary is also displayed. 

8. The observed location of the Venus bow shock is compared with two 

gasdynamlc models of flow about the planet. 

9. The terrestrial bow shock model of Fairfield (1971) is comapared with 
two models of flow about the earth. 

10. The observed location of the Mars bow shock is compared with two 

models of flow about the planet. 
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Abstract 


The earth's bow shock is frequently cited as an example of an astronhysical 
shock where particle acceleration is observed. However, because energetic particles 
observed upstream of the bow shock may be accelerated within the magnetosphere, 
it is important to understand the properties of the magnetospheric source. A 
first order picture of the spatial distribution of magnetospheric particles in 
the magnetosheath and upstream is obtained by mapping those magnetic field lines 
which drape over the magnetopause through the bow shock. Subsets of these field 
lines that connect to potential sites of magnetic merging on the magnetopause are 
also traced in the event that leakage occurs preferential ly where normal components 
of the field are present across that boundary. The results can be used to determine 
whether the so-called "diffuse" particles observed upstream are accelerated locally 
or within the magnetosphere. 



Introduction 


The energetic particles that are observed in interplanetary space are attributed 
to a variety of sources, both nonlocal and local. The primary nonlocal sources are 
the as yet uncertain sources of the galactic cosmic rays, while the more numerous 
local sources include solar flares, corotating interaction regions, interplanetary 
shocks, end planetary bow shocks. In recent years, much attention has been focussed 
on the ^50 keV ions detected upstream of the earth's bow shock because of their 
impl ications. for shock acceleration of galactic cosmic rays in the interstellar 
medium (cf. Axford, 1981). While the question of the contribution of magnetospheric 
particles to this pooulation has been addressed by a number of authors (i.e.. Scholar 
et al . , 1981), it is still a relatively unsettled issue, ‘Yet, if planetary bow 
shocks are to be used as astrophysical laboratories, it is important to understand 
the nature of the magnetospheric source. 

Ions and electron beams with energies a few times that of solar wind particles 
have been traced back to the quasiperpendicular bow shock (cf. Gosling et al., 1981). 
The origin of these few keV particles is most likely related to the reflection of 
a small fraction of the incident solar wind from the bow shock as first suggested by 
Sonnerup (1969). These beams are probably also responsible for the generation by 
a beam-plasma instability of the HMD waves that are observed in the same upstream 
region (Barnes^ 1970: Gary, 1981; Russell and. Hoppe, 1983). A more energetic, -■ 
(>50 keV) less directed population, the so-called "diffuse" particles which ere 
observed downstream of these beams (cf. Ipavich et al., 1981; Anderson, 1981; 

Kitchell et al . , 1983) has been explained as beam particles which are stochastically 
accelerated by the Fermi process in the regions where the upstream magnetic field 
fluctuations produced by the beam are converted against the shock (cf. Jokipii, 1967; 
Lee, 1982). Indeed, both West and Buck (1976) and Crooker et al . (1981) found an 
energetic ion population in the magnetosheath that appeared to be converted, together 



with fluctuating magnetic fields, along plasma stream lines from the region of the 
bow shock where the diffuse particles are observed upstream. 

On the other hand, very energetic (MeV) Jovian electrons observed in the inter- 
planetary medium have been attributed to a magnetospheric source (McDonald et al., 
1975), and Zwickl et al . (1981) found heavy ions of unmistakably magnetospheric 
origin in the upstream particles near Jupiter. Still, it may be argued that the 
region of space near the planets is contaminated by cold planetary ions which are 
then accelerated at the planetary bow shock, as proposed for the earth's upstream 
diffuse particles. 

These things considered, it is relevant to review what has been observed near 
the terrestrial magnetopause through which the particles must exit the magnetosphere. 
Williams et al. (1981) have examined the behavior of energetic particles at the 
magnetopause boundary. These authors demonstrated that bite-outs occur in the 
trapped pitch angle distributions at large pitch angles, which is what one expects 
if some of the drifting trapped population is lost to the magnetosheath. Also, 

West and Buck (1976), Asbridge et al . (1978) and Bieber and Stone (1981) have all 
reported observations of layers of energetic particles in the magnetosheath near the 
magnetopause. Scholer et al. (1981) showed that the flux in the mapnetosheath is 
comparable to the 90° trapped flux which is feeding the leakage. At least some of 
these particles observed near the magnetopause must end up in interplanetary space. 

The distinction between magnetospheric particles and those accelerated upstream 
of the terrestrial bow shock- was attempted experimentally by Scholer et al. (1981) 
who argued that magnetospheric ions are more energetic and are accompanied by 
electrons. In another observational study, Bieber and Stone (1981) pointed out that 
upstream magnetospheric electron events usually occur in conjunction with geomagnetic 
activity on the ground. Yet it seems that some of the confusion surrounding this 
m.atter of magnetospheric sources can be cleared up most effectively by a model that 
describes where those particles which are observed to leal: out of the magnetopause go. 



A global picture of the magnetospheric particle source can be expected to 
exhibit certain characteristics. The particles are presumably accelerated within 
the magnetosphere either in the routine manner of the permanently trapped outer zone 
populations or in transient events in the magnetotail. Proceeding outward, there 
are two mechanisms by which these particles can exit the magnetopause based on the 
current understanding of its structure. One possibility is that a turbulent . 
boundary layer (cf. Russell and Greenstadt, 1983) can scatter an adiabatical ly 
drifting trapped magnetospheric particle onto magnetosheath field lines. Provided 
that the boundary layer is large in areal extent, this type of leakage would occur 
over practically the entire magnetopause, with perhaps a preference for dawn and 
dusk where the sunward drifting electrons and protons, respectively, cross the 
terminator. Alternatively, the particles may leak out preferentially at areas on 
the magnetopause where merging of the magnetospheric and magnetosheath fields occurs 
(Crooker, 1979; Luhmann et al., 1983) because there normal components of the magnetic 
field allow adiabatic motion from one region to the other (Speiser et al., 1981; 

Daly, 1983). 

Next, one must trace the 'particle motion through the magnetosheeth from the 
magnetopause to the bow shock. Here there is a question of whether the particle 
motion is adiabatic or diffusive. If along its path a particle encounters magnetic 
fl uctuations. of substantial amplitude which satisfy the condition for gyroresonance 
when Doppler shifted to the particle's reference frame (cf. Jokipii, 1967), it will 
be scattered by the fluctuations. In this case the magnetospheric particles will 
"diffuse" outwards towards the shock. Moreover, if the magnetic fluctuations have 
some net motion, say in the direction of the convecting magnetosheath plasma, the 
isotropized particles will exhibit a similar bulk motion. If the scattering is 
strong, few magnetospheric particles will be able to reach the upstream region.' 



Rather, they will be carried antisunwards along the flanks of the magnetopause with 
the magnetosheath plasma in which the magnetic fluctuations are embedded. However, 
if the magnetic fluctuations are in a frequency range that does not affect the 
particles of interest, or if the magnetic field is fairly undisturbed, the gyro- 
centers of the energetic particles will approximately follow magnetosheath field 
lines through the bow shock to interplanetary space.- While displacement of the 
gyrocenters from the field lines will be caused by the motional (VXB) electric field 
drift associated with the convection of the magnetized magnetosheath and solar 
wind plasma, this can probably be neglected for >50 keV- ions (cf. Anderson, 1981; 
Mitchell et al . , 1983). Thus, field line tracing from the magnetopause through the 
magnetosheath will, under conditions of nearly adiabatic particle motion, give 
information (to within a gyroradius) about the volume that should be populated with 
magnetospheric particles. 

While some combination of diffusive and adiabatic behavior is probably a more 
realistic description of the particle motion in the magnetosheath, several studies 
(i.e.. Palmer, 1981) suggest that the mean free path for scattering in the magneto- 
sheath can be long. In this paper emphasis is placed on the characteristics of 
the magnetospheric source that would be expected for adiabatic behavior of the 
particles within the megnetosheath, although some discussion of the effects of 
magnetosheath turbulence is included in the final section. Toward this end, magne- 
tosheaih magnetic field lines that either pass within ^'•■1 ('vl particle gyroradius) 
of the magnetopause (for modeling widespread leakage) or at potential sites of 
magnetopause reconnection where the magnetosheath and magnetosphere fields are 
entiperallel, are traced outward to interplanetary space. The magnetosheath field 
model used here is derived from the gasdynamic treatment of the solar wind flow 
around the magnetosphere (Spreiter and Stahara, 1964, 1982), In this model, the 



magnetic field configuration is calculated on the assumption that the interplanetary 
magnetic field is frozen in the diverted solar wind plasma. Thus, the "beams" of 
magnetospheric particles are found to have distinctive shapes and locations that 
depend on the interplanetary field orientation. 

Description of the Model 

The gasdynamic magnetosheath magnetic field model has been described elsewhere 
by its originators (Spreiter and Stahara, 1964, 1982). As mentioned above, this 
model assumes that the interplanetary magnetic field is frozen in the flow. 

Although few comparisons of this model with observations have been carried out for 
the earth (cf. Fairfield, 1976), several detailed analyses have been done for the 
Venus magnetosheath (Spreiter and Stahara, 1980). The fact that the Venus ionospheri 
obstacle is much less compressible than the earth's magnetosphere, and is scaled 
much smaller, is an advantage for comparisons with the steady state model. The 
time scale for variations in the interplanetary field is typically short compared 
to the time (several hours) it takes a spacecraft to travel through the dayside 
magnetosheath at the earth. But the interplanetary field can be quite steady during 
the time interval of the Pioneer Venus orbiter traversal of the Venus magnetosheath 
(^-.5 hr). The magnetcsheath magnetic field model based on gasdynamics agrees quite 
well with the observed magnetic field at Venus. Although here the object of interest 
is the less well-behaved earth, it is considered that an idea of the nagnetosheath 
field line geometries can be obtained from an examination of constant scale models 
with constant interplanetary fields. A sonic Mach number of 6 and a magnetopause 
shape derived from the Hedgecock and Thomas magnetosphere model (cf. IJalker, 1976) 
were used in the gasdynamic calculation to produce the magnetosheath fields used here 

To determine the locations of magnetic merging on the magnetopause, the magneto- 
spheric field at the" magnetopause was presumed to be given by the Hedgecock and 
Thomas (cf. Walker, 1976) model. As described elsewhere (Luhmann et al . , 1983) the 



locations where this internal field was within '^lO® of antiparallel to the maoneto- 
sheath field were then defined as urging sites. 

Routes for magnetospheric particle leakage are approximated by magnetosheath 
field lines which either draoe against. the dayside magnetopause, for the case of 
widespread leakage, or by field lines which originate at the magnetopause at the 
aforementioned merging sites. 

Figures 1 and 2 illustrate the results of the modeling in the same format for 
the cases of widespread leakage and localized reconnection site leakage, respectively. 
Th'ese figures display projected views of the three-dimensional magnetosheath field 
line configuration for different interplanetary field orientations. The underlying 
magnetosheath draping is shown as fine lines. These field lines start at 49 points 
in interplanetary space on a square grid oriented perpendicular to the interplanetary 

t 

field. The heavy lines in Figure 1 are projected views of field lines which pass 
within ^-1 earth radius of the magnetopause at the terminator plane at equi spaced 
intervals of -vlO®.' These produce a three-dimensional picture of the volume of 
flux tubes which would be populated by magnetospheric particles if the particles 
leak out over a large area of -the dayside magnetopause. The shapes of the volumes 
are notably different for each interplanetary field configuration, appearing as 
sheets for interplanetary fields perpendicular to the flow, and as a cylindrical 
,tube for radial magnetic field. Subsets of these .volumes are connected to merging 
sites, as defined above, on the magnetopause. The heavy lines in Figure 2 identify 
groups of field lines which originate near those sites, which are shown on the 
magnetopause from the viewpoint of the sun in Figure 3. Figure 2 illustrates the 
importance of southward interplanetary field in obtaining large amounts of magneto- 
spheric particle leakage if merging at the magnetopause is a necessary factor. 

Finally, Figure 4 shows the cross-sections of several of the volumes of flux tubes 
from Figures 1 and 2 at various planes in interplanetary space. This diagram gives 



a feeling for the sizes of the volumes of magnetospheric particles that a spacecraft 
might traverse while located beyond the bow shock. 

Discussion and Conclusions 

In the introductory section the effects of scattering of the particles by 
magnetic fluctuations and of drifts produced by the motional electric field 
(E = VXB where V = plasma velocity, B = magnetic field) were mentioned. Anderson 
(1981) has given a complete description of trajectory dispersion by the VXB field 
in interplanetary space where V and B are uniform. However, the drift effects in 
the magnetosheath are not so simply described because of the complicated magneto- 
sheath velocity field and the magnetic field, which varies with the interplanetary 
field orientation. Figure 5 shows several examples of particle trajectories 
originatin-g near the subsolar magnetopause, drifting in the VX'B and F fields from 
the gasdynamic magnetosheath model. As discussed by Anderson (1981) for the upstream 
region, the particle gyrocenter trajectories are dispersed antisunward from magnetic 
field lines according to their parallel (to B) velocities. Particles with small 
parallel velocities are swept into the flow, in the limit of zero parallel velocity 
their gyrocenters following streamlines. At the opposite extreme, particles with 
large parallel velocities have gyrocenter paths that lie practically along magnetic 
field lines. Most observations are of particles with intermediate behavior. Thus, 
some spatial dispersion of particles of different energies in the magnetosheath 
and upstream regions is to be expected. The anticipated spatial gradients, with 
low energies mostly antisunward, should be considered in interpreting delays in 
the onsets of fluxes at successively higher energies as observed from spacecraft 
moving from dawn toward noon. This spreading of the magnetospheric beam will distort 
the spatial distributions suggested by Figures 1-2 for particles with low parallel 
velocities. However, for energies >50 keV and pitch angles >45® field-aligned 
gyrocenter motion is probably a fair approximation. 



On the subject of scattering, it has been established that magnetic fluctuations 
O-pa^i" 

in the magnetosheath, ^from the turbulent boundary layer near the magnetopause 

(Hones, 1983), are confined to streamlines of magnetosheath plasma flow that are 
connected to the quasiparallel bow shock (Greenstadt et al., 1983, Russell and 
Luhmann, 1983). Figure 5 illustrates how sections of the trajectories of the 
magnetospheric particles intersect the regions (shaded) of fluctuating magnetic 
fields behind the quasiparallel bow shock on their way through the nagnetosheath. 

The spectrum of the magnetosheath field fluctuations, which has been studied by 
Greenstadt et al . (1983) is similar in shape, and at least equivalent in power, to 
the spectrum of the waves upstream from the bow shock. The latter has been shown 
by Lee (1982) to often justify the use of the diffusion equation in describing the 
transport of upstream energetic particles. The upstream patterns for magnetospheric 
particle leakage shown in Figures 1 and 2 will under such circumstances be smeared 
out by the diffusing regions in the magnetosheath. However, some properties of the 
leakage patterns, such as the local time sector of the transmitted population, may 
be roughly preserved. An important point in connection with this geometry is that 
the upstream leakage frequently occurs through the region of the quasiparallel bow 
shock, where the so-called "diffuse" particle populations, consisting of nearly 
isotropic energetic ions, are observed (Ipavich et al . , 1981, Paschmann et al . , 1981, 
Bonifazi end Moreno, 1931). - The present results suggest that caution__should be 
exercised when interpreting these diffuse populations as ions locally accelerated 
in the upstream region. Other properties of the observed energetic diffuse ions 
that a magnetospheric source would explain are the apparent limiting fluxes, which 
ere similar to the trapped flux near the magnetopause, the time constants for flux 
buildup after a sudden change in the interplanetary field, which may be the leakage 
end travel times to the upstream observing point, and the deoendence of the 



energetic upstream population (West and Buck, 1976), since the trapped population 
at the magnetopause as well as the rate of leakage through reconnection sites 
(i.e., the areal extent of merging regions) will be greater under the conditions 
that produce high Kp. The observed frequent absence of energetic upstream electrons 
together with the diffuse ions can be attributed to leakage efficiencies and 
propagation differences for the two species, partly due to their very different 
gyroradi i . 

One obvious test that could be performed to determine the effeciency of upstream 
acceleration without contamination by the magnetospheric source is to measure 
energetic particles upstream of the bow shocks of planets with v;eak intrinsic fields 
like Venus and Mars, However, although magnetic field measurements are available 
for these planets, energetic particles were not measured on the spacecraft missions 
to them. In any case, the present study suggests that in interpreting observations 
of energetic particles near planetary bow shocks, one must consider not only whether 
the local magnetic field lines connect to the bow shock, but also whether they 
connect to the magnetopause. 
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Figure Captions 


Figure 1. Projections (in the GSE coordinate system) of magnetosheath field lines 
(fine lines) and field lines passing within -^-l Re of the magnetopause in 
the terminator plane (heavy lines) for various interplanetary field orientations 

(a) northward or southward (Z-directed), (b) dawn or dusk (Y-directed) , (c) 
Parker Spiral (45° to X and Y axes, parallel to X-Y plane, (d) radial 
(X-di rected) . • The gyrocenters of energetic particles which leak out of the 
magnetosphere over a large area of the magnetopause will follow paths roughly 
.within the volumes defined by the heavy lines. 

Figure 2. Same as Figure 1, but here the heavy lines originate near potential 
merging sites on the magnetopause (see Fig. 3), and the interplanetary field 
is: (a) southward, (b) northward, (c) Y-directed, (d) radial, (e) and (f) 

Parker Spiral toward and away, respectively. 

Figure 3. Projection of the magnetopause surface (in the GSE Y-Z plane) showing 
where the magnetospheric and magnetosheath model fields are within 10° of 
antiparallel at the magnetopause. Reconnection or merging is presumed to 
occur between the two fields at- these sites. Merging produces normal components 
of the field on the boundary, allowing the adiabatic motion of particles from 
the magnetosphere to the magnetosheath. The field orientations are (a) radial, 

(b) southward, (c) Y-directed, (d) Parker Spiral away, (e) northward, (f) 

Parker spiral toward. 



Figure 4. Cross-sections of some of the volumes defined by the heavy lines in 
Figures 1 and 2 at various planes parallel to the GSE coordinate axes. The 
light shading shows the areas v/ithin which widespread leakage (see Fig. 1) 
would be observed, while the dark shading shows these areas for localized,- 
reconnection site leakage (see Fig. 2). The X, Y, or Z positions and inter- 
planetary field orientations are indicated in the upper right. 

Figure 5. Examples of particle trajectories in the model maqnetosheath magnetic 
and electric (VXD) field. On the left, the interplanetary field is perpendic- 
ular to the flow. The light lines show the streamline and background field 
Tine geometry in the plane of symmetry. The particles launched near the 
subsolar magnetopause with 45° pitch angle, have energies of 1 keV, 5 keV 
and 50 keV from left to right. The shading shows where magnetosheath turbulence, 
which can scatter the particles, is expected behind the quasiparallel shock. 

On the right a similar diagram illustrates 1 keV and 5 keV particles for a 
Parker Spiral interplanetary field. 
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Abstract 


Models of the magnetospheric and nagnetosheath magnetic fields are used to 
determine the relative orientations of the two near the dayside magnetopause for 
the purpose of locating potential merging sites. Areas on the magnetopause with 
various degrees of antiparallelness for different interplanetary field orientations 
are displayed as contour diagrams. For southward and 6SE-Y interplanetary field,, 
the patterns obtained are consistent with those envisioned by Crooker in an earlier 
analysis which used simplified representations for the magnetic field geometry. 
Here, the application of realistic models shows the locations of areas where any 
antiperellel component occurs. Merging sites for radial interplanetary fields are 
also illustrated. The results suggest that the geometrical configuration of the 
field is suitable for merging over a large fraction of the magnetopause for inter- 
planetary fields that are either primarily southward, (ISE-Y, or redial (GSE-X) 
in direction. 



Introduction 


There is considerable interest in the question of where reconnection or magnetic 
field merging occurs in the magnetosphere (cf. Russell, 1976). Several years ago 
Crooker (1979) presented a qualitative model of the sites of magnetic field merging 
on the magnetopause. This analysis was based on a conceptual picture of the magne- 
topause and a superposed overlying uniform field, perpendicular to the earth-sun 

>- , 

line, which was taken to represent the magnetosheath field at its inner boundary. 

The latter was rotated to mimic the effect of different interplanetary field 
directions. Reconnection'or merging lines v^ere defined as the locus of points for 
which these approximate magnetospheric and magnetosheath fields were antiparallel 
when projected onto a plane (the GSE Y-Z plane). The major point arising from this 
analysis was the role of the cusp in locating potential merging sites. However, the 
author also qualified the results by -pointing out that merging can occur where only 
components of the tv/o- fields are antiparallel (cf. Cov;ley, 1976), and that the 
draping geometry of the megnetosheath field over the magnetopause is not necessarily 
well represented by a uniform field in the GSE Y-Z plane projected onto the magneto- 
pause. The present study extends the ideas put forth in this earlier effort by 
employing realistic models of the magnetospheric and magnetosheath fields to locate 
antiparallel fields at their magnetopause interface. The new aspects of this 
analysis include the location of fields with only components that are antiparallel, 
and the treatment of radial (GSE-X) interplanetary fields which were not considered"^ 
previously. 

Description of the Model . 

The Hedgecock and Thomas model (cf. Walker, 1976) was selected to represent the 
magnetospheric field-at the dayside magnetosheath because it has a fairly realistic- 
cusp geometry. This model also has north-south asymmetries caused by the uneven 
distribution of the data that went into its construction, rather than by true 



asymmetry. It was therefore assumed that the more heavily sampled northern section 
v;as the more accurate, and a mirror image of the model north of the equator was 
used to represent the southern magnetosphere. Possibly real dawn-dusk asymmetries 
in this northern section- were retained. 'The dipole axis of the model was assumed 
to coincide v/ith the 6SE-Z axis for this study. Although the seasonal variation of 
the dipole tilt with respect to this axis, will modify the magnetospheric field at 
the magnetopause, the complication introduced by the use of a tilted model precluded 
the study of this effect. 

The magnetosheath field was modeled using the assumption that the interplanetary 
field is frozen into the medium in the gasdynamic treatment of supersonic flow 
around an axisymmetric obstacle (cf. Spreiter and Stahara, 1982). The shape of the 
obstacle was presumed to follow the shape of the Hedgecock and Thomas model magne- 
topause in the noon-midnight meridian. For the purpose of the present analysis, 
the aberration of the solar wind flow was neglected. A free-stream sonic Mach 
number of 6.0 was used in the gasdynamic code, with which magnetosheath fields at 
a distance of approximately .5 earth radii from the surface of the obstacle were 
computed for a variety of interplanetary field orientations. (The. field determined 
with the gasdynamic code becomes inaccurate near the stagnation streamline that 
bathes the obstacle surface.) 

The angle between the two modeled fields was found over the surface of the 
dayside magnetopause at the grid points shown in Figure 1. Contour diagrams of the 
cosine of this angle were then constructed to display the regions of different 
degrees of antiperal lei ness as viev^ed from the Sun. 

Results 

Figure 2 contains the contour diagrams described above for various interplanetary 
magnetic field orientations (ie. magnetosheath models). The shaded areas, representing 
fields that are within 10° of being antiparallel, are roughly consistent with 



Crooker's (1979) merging lines for interplanetary field orientations perpendicular 
to the solar wind velocity. New sites arise when there is a substantial 6SE-X or 
radial component of the interplanetary field. As expected, the southward directed 
interplanetary field produces the largest area of nearly antiparallel fields on the 
magnetopause, while the northward interplanetary field produces the smallest area. 
However, the regions for which there is some component of the magnetosphere and 
magnetosheath fields antiparallel cover a large region of, the magnetopause for most 
interplanetary field orientations. 

Discussion and Conclusions 

For several reasons, the patterns shown in Figure 2 must be considered with 
some measure of caution when comparing with observations at the magnetopause. First, 
the interpi anetary field is typically variable on the time scale of plasma convection 
through the dayside magnetosheath. Except under circumstances of exceotional ly 
steady interplanetary field orientation, the actual pattern of antiparallel fields 
projected on the megnetcoause must be predicted from the field observed in the solar 
wind at earlier times. (The convection time is several minutes from the subsolar 
shock to the nose of the magnetopause.) Second, currents near the magnetopause and 
MHD effects such as that described by Zwan and Wolf (1979) will affect the field 
geometry at that boundary. Third, the possible dependence of the merging rate on 
the local plasma velocity has 'hot ieen folded into these patterns. If slow flow 
is necessary for merging then those regions near the subsolar stagnation point will be 
favored. Similarly, the possible role of the field magnitudes was neglected here. 
Finally, reconnection itself may cause the reconnected field lines to be pulled 
over the magnetopause, disrupting any pattern imposed purely by the magnetospheric 
and magnetosheath fields. (The boundary layer may be a manifestation of the latter 
process.) Yet, the patterns shown in Figure 2 may provide a picture of the initial 
conditions set up by a sudden change in interplanetary field orientation. 



Recent statistical studies of the- spatial distribution of flux transfer events 
near the magnetopause (Berchem and Russell, 1983) show a double-banded distribution 
which flanks the equator. These authors have argued that their results are con- 
sistent with the preferred occurrence of reconnection at the equator where quasi- 
steady merging has been predicted for southward interplanetary field (Dungey, 1961, 
Sonnerup, 1976). However, their distributions are also consistent v/ith a super- 
position of the patterns for southward field and for GSE-Y fields as shown in 
Figure 2. In fact in their study, some events occurred in association with practi- 
cally GSE-Y directed fields. Also, because the highest geomagnetic latitudes were 
not covered in their sample, the merging poleward of the cusps predicted for north- 
ward directed fields (see Fig. 2) could not have been observed. Still, separation 
of their data into events associated with steady, primarily southward interplanetary 
fields and GSE-Y fields may produce the distinctive patterns shown in Figure 2 for 
these tv.'o cases. One would expect the GSE-Y fields to produce events in opposing 
corners along a diagonal through the subsolar point, while southward fields would 
produce a broad latitudinal band including the equator. 

In conclusion, the present analysis represents the extent to which one can 
practically go in using static magnetic field models to understand the global picture 
of reconnection sites on the magnetopause. Time-dependent, 3-dimensional MHD models 
of the solar wind interaction with the magnetosphere, such as that developed by 
Wu et al. (1980) and Redder et al. (1981) v/ill, of course, be the ultimate tool for 
the investigation of 'merging sites. 



Figure Captions 


Figure 1. A view, in the GSE X-2 plane, of the grid points on the Hedgecock- 
Thotnas model magnetopause at which the dot-product of the magnetospheric 
field model and the gasdynamic magnetosheath field model v/as computed in 
this analysis. Also shown are the magnetospheric model field lines in the 
noon meridian. 

Figure 2, Contours on the magnetopause (viewed from the sun, i.e., the GSE Y-2 
plane projection) of equal value of the cosine of the angle between the 
magnetospheric and magnetosheath model fields. Only contours with negative 
values, implying some antiparallel component, are shown. Values at the 
contours, starting with the contour filled with shading, are -.98, -.95, -.9, 
-.8, -.7, -.6, -.5, -.4, -.3, -.2, -.1, 0. The interplanetary field 
orientations are as follows: (a) radial (+X-directed) , (b) southward, (c) 
Y-directed (tov/ard dusk), (d) Parker Spiral toward, (e) northward, (f) 

Parker Spiral away, (g) equal Y and 2 components, no X-component, (h) equal 
X, Y, and 2 components, (i) equal and opposite Y and 2 components, no X- 


component. 
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ABSTRACT 


On January 21, 1972 the Mars-3 spacecraft observed a variation in 
the magnetic field during its periapsis passage over the dayside of Mars 
that was suggestive of entry into a Martian magnetosphere. We have 
obtained the original data and trajectory (Sh. Sh. Dolginov, personal 
corrmuni cation, 1983), and have attempted to simulate the observed 
variation of the magnetic field using a gasdynamic simulation. In the 
gasdynamic model we first generate a flow field and then use this flow 
field to carry the interplanetary magnetic field through the Martian 
magnetosheath. The independence of the flow field and magnetic field 
calculation allows us to converge rapidly on an IMF orientation which 
would result in a magnetic variation similar to that observed by Mars-3. 
There appears to be no need to invoke an entry into a Martian magnetosphere 
to explain these observations. 



Of the magnetic fields of the four terrestrial planets, the magnetic 
field of Mars is least understood, even though Mars has been visited by 
many spacecraft. The reason for this gap in our understanding is that 
of the many U.S. spacecraft to Mars, only the first, Mariner 4, in 1965, 
carried a magnetometer, and on this mission the flyby distance was so 
great that only brief encounters with the planetary bow shock and 
magnetosheath were observed. At least three Soviet Martian orbiters 
have carried magnetometers, but none of these spacecraft haVe provided 
low altitude (<1000 km) data or measurements directly behind the planet. 
Nevertheless, some of these measurements have been interoreted as indicative 
of the existence of an intrinsic planetary magnetic field. The Mars 5 
measurements behind, but not directly behind, the planet have been inter- 
preted as indicating a tail (Dolginov et al, 1976). The Mars 2 measure- 
ments have been interpreted as indicating a dipole movement with its axis 
in the equatorial plane (Smirnov et al , 1978). The Mars 3 measurements 
have been interpreted to indicate, on one occasion, an entry into a Martian 
magnetosphere (Dolginov; et al , 1973j Gringauz et al , 1974). However, 
for this interpretation to be correct, the dipole moment : would have to 
be aligned with the spin axis. 

These interpretations have been criticized on several counts. The size 
of the magnetic moment consistent with the reported Mars 5 tail encounters 
is smaller than that reported from the dayside Mars 3 entry (Russell, 1978a). 
Further, the putative tail entries are suspect because of the direction of 
the field in the region inferred to be "tail". The Mars 2 inferences are 
inconsistent with the Mars 3 inferences, but are consistent with the mapneto- 
tail being due to draped interplanetary magnetic field lines (Russell, 1981). 



We note that Venus has no discernible intrinsic field, yet still has a 
magnetotail formed from the interplanetary magnetic field (Russell et al, 
1981a). Finally, the Mars 3 measurements have been interpreted as passage 
through a magnetosheath, rather than a magnetosphere (Wallis, 1975; 

Russell, 1978b). Wallis (1975) suggested that the region that Dolginov 
et al (1973) had identified as magnetosheath could equally well be 
interpreted as the foreshock region. Russell (1978b) suggested that the 
putative magnetospheri c field had the draping pattern of magnetosheath 
field lines. At that time there was no convenient means to test the hypo- 
thesis. However, code development since that time now permits a relatively 
simple test of this hypothesis, and it is the purpose of this note to 
perform such a test. 



GAS DYNAMIC MODELING 


Numerical models of the solar wind interaction with planetary obstacl 
have now been in use for nearly two decades (cf. Spreiter et al, 1966). 
Code developments and technological change since that time have decreased 
the cost and increased the speed and flexibility of the models (Spreiter 
and Stahara, 1980a, b), and provided the capability of intercomparing 
observations and models along actual spacecraft trajectories. A feature 
of the gasdynamic model makes the model ideally suited for 'the task at 
hand. The flow solution is not affected by the magnetic field. The 
magnetic field is obtained by stretching the original upstream magnetic 
field with the derived flow field, but neither the magnetic tension nor 
the magnetic pressure react back on the flow. The flow solution is 
affected only by the upstream Mach number, the ratio of specific heats, 
and the shape of the obstacle. 

We use this magnetic field independence in our investigations below 
in the following manner. First, we calculated the flow field around the 
obstacle for a given Mach number, , and shape. We choose a magnetosonic 
liach number of 7.7 as being appropriate for average solar wind conditions 
at Mars (Russell et al 1982b). For an obstacle shape we use a Venus 
shape with an ionospheric scale height of 0.025. The solution is quite 
insensitive to this latter quantity as long as it is smaller than about 
0.1. The choice of )T, the ratio of specific heats^is not a simple matter, 
Fairfield (1971) and Zhuang and Russell (1981) both deduced that a value 
^ = 2 gives the best fit to the terrestial magnetosheath thickness. 

Slavin et al (1983) have found good agreement for Mars using If = 2 , and 
Tatrallyay et al (1983) have found that "jf = 1.9 orovides the best fit for 
Venus. However, Tatrallyay et al also find an Alfvenic Mach number depend 



ence for the value of ^ and that the best fit value of ^ is 1.6 above' an 
Alfvenic Mach nunter of 7. Thus, we might expect that the best fit "for 
Mars might not be the same as the terrestial and Venus best fits. 

On the other hand, the magnetic field strength at Mars, on the day 

being considered herein^ is large and hence one might expect that the Mach- 

number was lower than usual. In the absence of knowledge of 

the plasma density and bulk velocity, we will use the average Mach number 

7.7 and the previously determined best- fit ^ -value of 2. 

After the flow field solution is obtained, we can calculate the 
magnetic field along the trajectory of the spacecraft as if the magnetic 
field lines were dye lines in a fluid flow. This calculation is rapid and 
can be repeated at small cost in 'computing time. Thus, it is feasible 
to hunt for a suitable IMF direction, flow observation angle, or obstacle 
size. 

Previous work has used the gasdynamic simulation to determine whether 
the bow shock location of Mars implies an obstacle of the size of the 
planet plus ionosphere or somewhat larger (Russell 1977; 1970’, Slavin et al , 
1983). Once that obstacle size is determined, one has the further problem 
of whether the ionosphere is strong enough to stand off the solar wind. 

In this study, we use the gasdynamic simulation in a different way. We 
wish to examine whether we can replicate the Mars 3 magnetic field observa- 
tions, which were claimed to indicate a magnetospheri c entry, with a . 
trajectory passing solely through magnetosheath. This argument has been 
made in the past (Russell, T978b), but not as quantitatively as is now 
possible. 



CALCULATIONS 


The magnetic field observations that we are attempting to model are 
shown in Figure 1. Panel (a) shows the time series and panel (b) shows the 
field vectors and trajectories projected into the Y-Z plane. The field 
values and the trajectory have been supplied in tabular form by Sh. Sh. 

Dolginov (personal communication, 1983) to whom we are very grateful. The 
trajectory is very similar to that originally derived by Russell (1978b; 1979). 
However, use of investigator-supplied values removes the uncertainty 
associated with the original interpolation process and a possible source of 
error. Table 1 lists the position of Mars-3 at selected times. The magnetic 
field values also closely resembled those published by Gringauz et al . (1974) 
and Russell (1978b; 1979) as well as those in the early publications by 
Dolginov and coworkers. However, there is ambiguity and inconsistency in the 
various descriptions of the coordinate system used to display the Mars-3 magneti 
field data. We use herein the latest definition (Dolginov, 1978) in which 
the direction of the Y-component is in the direction of planetary motion. As 
we show below this definition of coordinates provides internal consistency 
in the behavior of the data. We emphasize that we must judge consistency in 
terms of the variations in the data, since the zero levels of the Mars-3 
magnetometer are not known. Our modeling provides an estimate of these 
uncertainties and they prove to be large. 

If we assume that the Mars-3 magnetometer's zero levels remained unchanged 
during its periapsis passage on 1/21/72, we can use the field variations to 
initialize our simulation effort. It is both our hypothesis and Wallis' (1975) 
that an inbound shock crossing occurred at 2129 + 2 (across an unfortunately 
timed data gap of 4 minutes duration) and an outbound crossing occurred at 2259. 



The latter time was also given as the outbound shock crossing by Dolainov et al 
(1975). However, there is little change in the vector field-direction there. 
This fact suggests -that the outbound shock is. a quasi-paral lei shock, i.e., 
the IMF is parallel to the shock normal at this point. The inbound shock has 
a more rapid increase, as if it were a quasiperpendicular shock. 

Table 2 lists the median maqnetic field vectors from 2115 to 2145 Moscow 
time across the feature we interpret to be the inbound shock crossing. The 
Mars-3 magnetometer returned 8 values of the magnetic field in rapid succession 
every 2 minutes only some of which were successfully recovered. The right- 
hand column labelledN lists the number of values available at each recording 
interval. Assuming that the noise on the measurements, both due to the ambient 
medium and due to the spacecraft, were randomly positive and negative, we used 
the median' values in each telemetry interval as our best estimate. Medians 
have been used because mean values are sensitive to occasional very inaccurate 
values whereas medians are not. 

For our first estimate of the IMF field, we take the jump in field across 
the inbound shock divided by three and assume this to be the IMF value. As 
shown in Table 2, the upstream values are "noisy". However, since we are 
merely attempting to determine a starting value we will simply use the value 
at 2127 as the upstream value and the value at 2131 as the downstream value. 
One-third of this difference is (-1.5, -4.5, 0.9) which provides our first 
estimate of the upstream field corrected for zero level errors. Figure 2 
panel (a) shows the time series derived using this vector &s input and an 
ionopause altitude of 430 km. The inbound shock crossing occurs about 8 
minutes early and the outbound shock about 22 minutes early. If we rotate 
the flow 9°. in the direction of aberration expected from planetary motion, we 
obtain the time series shown in Figure 2 panel (b). A 9° aberration angle is 



only about 5° larger than that expected due to Kars' orbital motion transverse 
to the solar wind flow and well within the known amplitude of solar wind 
directional variations. We note that other flow directions will also 
reproduce the observed shock locations for example a total aberration angle 
of 7° and a flow from the south of -4° will do as well. 

The model time series of Figure 2 resemble the variations seen in 
Figure 1. However, the model variations are smaller than the observations, 
especially in the component. Further, the jump in the magnetic field at 
the outbound shock seems larger in the models than in the data. Hence we 
need to make the IMF larger, more negative along the Z-di recti on and more 
aligned with the outbound shock normal. Figure 3 shows the results of using 
an input value -of (-4, -4.5, -3) nT and an aberration of 7 degrees with the 
flow coming from 4° below the ecliptic plane. The panel on the right shows 
the model and on the left the observations with the baselines corrected so 
that the average IMF readings agree. The qualitative behavior of each of the 
corresponding traces agrees. The magnitude of the variations are equal but 
the traces differ in the exact timing of the peaks and .the rapid variations 
in direction. There are many possible reasons for such differences. There 
are questions of the appropriate choice of y and ^MS for this day and of the 
accuracy of the gasdynamic simulation of an MHD problem. We may not have 
made the best choice of the IMF orientation. Further, it is certain that 
the IMF did not remain steady during the Mars flyby, as assumed in the 
simulation. 

In order to determine whether our solution is, in fact, a reasonable 
one, Figures 4, 5 and 6 compare individual vectors from the observations and 
model. Figure 4 shows the solar ecliptic plane projection. The field vectors 
to the right and left of Mars agree quite well. Only during the subsolar 



passage is there a major disagreement. There is certainly no hint in these 
data that a magnetosphere has been entered. Figure 5 shows the vectors 
projected in the dawn-dusk plane. There is very good agreement in this 
projection except near the inbound shock crossing, where there is some 
apparent temporal variation. Finally, in Figure 6, we show the two solar- 
cylindrical projections. Again, the agreement is very good, except for 
obvious temporal fluctuations and in a region near closest approach. This 
latter distortion is a pulling of the field toward the anti solar direction 
as if there were some drag on the flow in this region. Such a drag could be 
provided by mass-loading. However, calculations suggest that mass-loading 
is weak at Mars (Russell et al . , 1983). 

Earlier, we speculated that the inbound shock was a quasi perpendicul ar 
shock and the outbound 'shock was quasiparallel. This is, of course, implicit 
in the solutions we have derived, but it is of some interest to check this 
assumption. We can use the coplanarity assumption to derive the shock normal 
from the model field values. These normals are (.144, -.624, -.442) and 
(.618, .705, .347). Since the TMF value input to the model was (-4, -4.5, 
-3), these correspond to values of 77" and 12°, respectively, for the 

inbound and outbound shocks. These values are those of a quasiperpendicular 
shock and a quasiparallel shock, as we originally assumed. 

DISCUSSION AND CONCLUSIONS • 

The model values do not- perfectly replicate the data. Part of the 
differences are certainly due to temporal variations in the interplanetary 
magnetic field. Some of the difference may be attributable to non-gasdynamic 
effects, such as magnetic forces and mass loading. In view of the existence 
of so many reasons why there should be di fferences the observed agreement- 
is heartening. We have found a magnetosheath magnetic field along the 
Mars-3 trajectory which resembles the Mars-3 observations. Thus, the Mars-3 



data do not provide unambiguous evidence for an intrinsic Martian magnetic 
field. Our best fit obstacle size Is also quite consistent with an iono- 
spheric obstacle. Not only is the solar wind much weaker at Venus than at 
Mars, but also the Mars gravitational field is weaker. Hence, the upper 
atmosphere and ionosphere of Mars will have much greater scale heights than 
at Venus, all else being equal. On the other hand, in situ measurements of 
both the field and plasma of the Martian ionosphere will be necessary before 
the exact' nature of the ionopause is determined.' This study says only that 
any planetary field did not extend much above 400 km on 2/21/71, and certainly 
not to 1,200 km, the location of Mars-3. 

Finally, we have shown that the gasdynamic simulation provides a powerful 
tool for interpreting magnetic measurements obtained in planetary magneto- 
sheaths. In the past this important capability has not been fully exploited. 
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Table T. Mars-3 Trajectory 


Moscow Time 
;2101 
2131 
2201 
2231 
2301 


XSE 

-0.177 

0.659 

0.956 

0.240 

-0'.641 


YSE 

-2.788 

-1.313 

0.758 

2.198 

2.927 


ZSE 

-0.747 Rj^ 
-1.042 
-0.639 
0.499 
1.557 



Table 2. Median Magnetic Field Vectors 


Moscow Time 

BX • 

BY SE"^ 

BZ SE"^ 

N 

2115 

-0.6 

3.3 

-6.2 

6 

2117 

-2.4 

3.9 

-7.8 

8 

2119 

-2.3 

3.9 

-5.5 

8 

2121 

-2.0 

4.5 

-6.2 

8 

2123 

-2.6 

5.5 

-5.2 

3 

2125 

-3.9 

7.8 

-7.8 

5 

2127 

-3.9 

5.2 

-7.8 

1 

2129 

-6.2*. 

, -1.6* 

-6.5* 

0 

2131 

-8.4 

-8.4 

-5.2 

3 

2133 

.- 9-1 

-10.4 

-7.8 

2 

2135 

-11.0 

-13.0 

-9.1 

8 

2137 

-10.4 

-13.0 

-13.6 

3 

2139 

-10.4 

-12.6 

-16.9 

8 

2141 

-10.4 

-12.3 

-20.2 

5 

2143 

-10.4 

-8.4. 

-23.4 

5 

2145 

-8.1 

-5.2. 

-22.1 

4 


"*^Zero levels unknown 

* 

Interpolated value 




Fi qure 1 . Magnetic field observed by the Mars 3 soacecraft on 01/21/72. 

(a) Time series in solar ecliptic coordinates. 

(b) Vectors along the trajectory projected into the Y-Z solar ecliotic 
plane. No zero level corrections have been applied. 


Figure 2. Magnetic field calculated along Mars 3 trajectory using a 

ratio of specific heats, y» of 2 a nannetosonic Mach number 
of 7.7, and an iononause altitude of ,430km, with an IMF of 
C-I.g, -4.5, 0.90) nT. 

(a) Model A with no aberration of the solar wind. 

(b) Model B with a 9^' aberration of the solar wind in the-di rection of 
the average expected solar wind observation. 


Figure 3. Magnetic field along the Mars 3 trajectory 

(a) Panel on left shows observations corrected so that average IMF is 
same as IMF in model C. 

(b) Panel on right shows model C. All parameters are the same as in 
models B, except IMF has been changed to (-4, -4.5, -3). 


Figure 4. Ecliptic plane projection of magnetic field along Mars 3 
trajectory. 

(a) Corrected observations. 

(b) Model C. 

Figure 5 . Dawn-dusk plane projection of magnetic field along Mars 3 
trajectory. 

(a) Corrected observations. 

(b) Model C 

Figure 6 . Solar-cylindrical projection of magnetic field along Mars 3 
trajectory. 

(a) Corrected observations. 

(b) Model C. 
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